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Executive  Summary 


A  new  technique  to  delay  transition  to  turbulence  in  hypervelocity  air  flows  is  introduced 
and  investigated  in  this  work.  The  main  motivation  for  such  a  technique  is  the  pressing 
need  to  reduce  aerodynamic  heating  rates  on  hypersonic  vehicles.  Turbulent  heat  transfer 
rates  can  be  an  order  of  magnitude  higher  than  laminar  rates  at  hypersonic  Mach  numbers. 
Hence,  schemes  to  delay  transition  to  turbulence  could  provide  an  important  means  to 
reduce  the  heating  rates  on  hypersonic  vehicles. 

The  technique  presented  here  makes  use  of  the  recently  discovered  behavior  of  high- 
enthalpy  nonequilibrium  CO2.  Experimental  and  computational  data  show  that  for  the 
same  stagnation  enthalpy,  the  transition  Reynolds  number  is  larger  for  CO2  flows  than  for 
either  air  or  N2  flows.  The  explanation  for  this  phenomenon  lies  in  the  fact  that  when  CO2 
is  in  vibrational  and  chemical  nonequilibrium,  these  relaxation  processes  absorb  energy  from 
the  acoustic  disturbances  whose  growth  in  the  boundary  layer  is  responsible  for  transition  in 
hypervelocity  flows  (i.e.,  the  second  or  Mack  mode).  In  light  of  this  behaviour,  the  injection 
of  CO2  into  the  boundary  layer  of  interest  is  explored  as  a  transition-delaying  technique 
in  this  work.  The  particular  injection  geometry  investigated  consists  of  a  series  of  round 
orifices  located  at  the  upstream  end  of  the  test  model.  The  test  model  is  a  5  degree  half 
angle  cone.  This  geometry  was  chosen  due  to  the  wealth  of  previous  information  available 
with  which  the  transition  Reynolds  number  can  be  compared.  Seven  injector  geometries 
were  built  in  which  the  injection  angle,  the  number  of  orifices  and  the  diameter  of  the 
orifices  were  varied. 

Five  experimental  series  were  conducted.  Initial  experiments  used  mixtures  of  CO2/N2 
as  the  test  gas,  rather  than  introducing  CO2  directly  to  the  boundary  layer  of  the  cone.  This 
allowed  the  effect  of  the  presence  of  varying  amounts  of  CO2  on  transition  to  be  isolated 
from  any  specific  injection  scheme.  The  obtained  results  demonstrated  that  the  addition 
of  CO2,  even  as  a  minor  component  of  the  test  mixture,  significantly  delays  the  onset  of 


transition.  For  example,  in  the  case  of  60%N2/40%C02  by  mole  fraction,  the  value  of 
the  transition  Reynolds  number  was  more  than  double  that  of  the  case  with  100%  N2.  A 
similar  effect  was  noted  in  several  experiments  using  mixtures  of  air  and  CO2  as  the  test 
gas.  Several  runs  were  then  performed  in  100%  air  to  determine  a  baseline  condition  for 
injection  experiments.  The  requirement  for  this  condition  was  for  transition  to  occur  near 
the  mid-point  of  the  cone,  so  that  any  adjustment  to  the  transition  location  due  to  injection 
would  be  easily  observed.  Finally,  tests  were  carried  out  with  the  purpose  of  visualizing 
the  injection  process  to  provide  confirmation  that  the  triggering  of  the  injection  system  was 
appropriate  to  provide  adequate  CO2  by  the  time  of  the  arrival  of  the  main  test  flow.  The 
timing  of  the  CO2  injection  was  successfully  verified  and  it  was  established  that  the  flow  of 
C02  begins  approximately  100  ms  before  the  main  flow  arrives. 

As  this  project  has  progressed,  a  switch  has  been  made  from  an  old  instrumented  cone 
used  by  previous  researchers,  which  was  adequate  to  launch  the  project,  to  a  new  cone 
model  having  approximately  four  times  the  number  of  thermocouples.  The  thermocouples 
are  now  uniformly  distributed  in  the  circumferential  and  axial  directions.  This  will  allow 
us  to  pinpoint  the  onset  of  transition  more  accurately.  Also,  based  on  the  results  obtained 
thus  far,  a  new  injector  tip  has  been  designed  and  built,  having  approximately  twice  as 
many  holes  as  previous  ones.  This  will  allow  lower  injection  pressures  to  be  used  for  the 
same  mass  flow  rate. 
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Chapter  1 

Introduction 


A  novel  method  for  delaying  transition  in  hypervelocity  flows  of  air  by  injecting  CO2  into 
the  boundary  layer  of  interest  is  investigated  in  this  work.  Experimental  and  numerical 
data  show  that  when  pure  CO2  is  in  vibrational  and  chemical  non-equilibrium,  these  relax¬ 
ation  processes  absorb  energy  from  acoustic  disturbances  whose  growth  is  responsible  for 
transition  in  high  enthalpy  flows  [9;  10;  13;  14  .  This  is  an  inviscid  instability  mode  (sec¬ 
ond  or  Mack  mode  [16J )  associated  with  wave-like  disturbances  that  are  trapped  within  the 
boundary  layer.  In  contrast,  the  first  mode,  dominant  in  low-speed  flows,  is  a  viscous  mode 
associated  with  Tollmien-Schlichting  waves  [21J.  By  absorbing  energy  at  the  same  frequen¬ 
cies  as  the  acoustic  disturbances,  nonequilibrium  CO2  delays  transition  in  hypervelocity 
flows  as  compared  to,  in  particular,  N2  and  air. 

In  the  present  investigation,  several  injection  pieces  have  been  designed,  built  and  in¬ 
corporated  into  an  injection  system  for  an  existing  slender  cone  model.  Several  series  of 
experiments  have  been  performed  in  the  T5  hypervelocity  shock  tunnel  at  Caltech  using 
this  cone  model.  Exploratory  testing  was  carried  out  to  determine  the  effect  of  adding  an 
initial  component  of  CO2  to  the  test  gas  for  both  air  and  N2  flows,  with  delayed  transi¬ 
tion  observed  in  both  cases.  Further  testing  was  performed  to  identify  a  suitable  baseline 
condition  in  air  for  the  injection  experiments:  the  requirement  for  this  condition  was  that 
transition  occurred  approximately  midway  along  the  surface  of  the  cone.  A  limited  number 
of  injection  tests  were  then  performed,  but  questionable  instrumentation  readings  meant 
that  no  firm  conclusions  could  be  drawn  from  these.  Finally,  visualizations  of  the  injected 
flow  were  made  with  a  high  speed  digital  camera  to  verify  that  the  timing  of  injection  was 
appropriate  for  the  requirements  of  the  project,  i.e,  that  adequate  injection  flow  would  be 
established  by  the  time  of  arrival  of  the  main  flow  in  the  test  section  during  a  shot. 
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2.5  Schematic  of  the  injection  components  within  the  cone  model:  (top)  the  up¬ 
stream  end  of  the  cone,  which  holds  the  cone  tip  injector;  (below)  the  down¬ 
stream  end  of  the  cone,  showing  the  flange  which  holds  the  rigid  metal  tube 
for  transporting  the  CO2.  Note  that  the  test  configuration  differed  slightly 


from  that  shown .  13 

2.6  Schematic  and  drawings  of  the  cone  tip  injector  piece .  17 

2.7  Operational  amplifier .  18 

2.8  Timing  circuit  diagram .  19 

2.9  Relay  circuit  diagram .  20 

2.10  Run  tank  pressure  trace  from  shot  2448,  together  with  the  profile  obtained 

from  the  present  calculation .  29 


2.11  Schlieren  visualizations  from  injection  tests:  (top)  injection  into  ambient  air 

with  an  initial  run  tank  pressure  of  149  psi;  (bottom)  injection  into  evacuated 
T5  dump  tank  with  an  initial  run  tank  pressure  of  439  psi.  In  each  case,  the 
time  from  triggering  is,  from  left  to  right:  70,  183,  and  300ms .  30 

2.12  Run  tank  pressure  traces  from  shot  2448  and  from  an  ambient  injection  test. 

Both  the  original  test  signal  and  a  scaled/shifted  version  are  plotted:  the  latter 
is  scaled  on  the  pressure  axis  so  that  the  initial  pressures  match,  and  shifted 

on  the  time  axis  (by  182  ms)  so  that  the  initial  parts  of  the  profiles  overlap.  .  32 

2.13  The  run  tank  pressure  trace  from  an  ambient  injection  test,  with  voltage  traces 
from  two  microsensors  attached  to  the  ball  valve,  showing  the  timing  of  the 
valve  opening.  £=0  corresponds  to  the  trigger  signal  being  sent  to  the  injection 
system.  Schlieren  images  from  various  points  in  the  injection  sequence  are 
shown:  (top  left)  £=30  ms,  opening  of  the  valve;  (top  right)  £=62  ms,  startup 
of  the  injection  flow  at  the  cone  surface;  (bottom  left)  £=183  ms,  beginning  of 
test  time  if  main  flow  were  present;  (bottom  right)  £=300  ms,  injection  flow  at 

a  later  time .  33 

2.14  Schlieren  visualisations  from  low  enthalpy  shots  with  various  run  tank  pres¬ 
sures:  (top  left)  361  psi;  (top  right)  162  psi;  (middle  left)  50  psi;  (middle  right) 

16  psi;  (bottom)  no  injection. 
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1.1  Motivation 

Before  engineers  can  reliably  design  devices  to  survive  flight  through  or  into 
an  atmosphere  at  hypersonic  speeds,  they  must  somehow  provide  for,  avoid,  or 
otherwise  accommodate  the  enormous  heat-transfer  rates  to  the  vehicle  engen¬ 
dered  by  such  flight  speeds. 

-William  H.  Dorrance,  1962.[ 3j 

The  main  justification  for  this  work  is  the  pressing  need  to  reduce  aerodynamic  heating 
rates  on  hypersonic  vehicles,  including  the  configurations  being  studied  for  Rocket  Based 
Combined  Cycle  (RBCC)  propulsion.  As  a  simple  approximation,  hypersonic  aerodynamic 
heating  grows  as  the  cube  of  the  free-stream  velocity  and  linearly  with  the  free-streanr 
density  |2j.  Furthermore,  turbulent  heat  transfer  rates  can  be  an  order  of  magnitude  higher 
than  laminar  rates  for  hypersonic  Mach  numbers  [12  .  Hence,  for  high  velocities,  heating 
loads  become  an  important,  if  not  the  dominant,  consideration  in  the  design  of  hypersonic 
vehicles  [2  .  A  reduction  in  heating  loads  means  less  thermal  protection  needed  and  hence 
less  structural  mass,  or  conversely  greater  payload  for  a  given  thrust.  While  one  might 
not  have  control  over  the  free-stream  velocity  for  a  given  application,  heating  rates  could 
hypothetically  be  reduced  by  keeping  the  boundary  layer  laminar  over  larger  portions  of 
the  vehicle  using  schemes  such  as  the  one  proposed  here. 

1.2  Relevant  Previous  Work 

Hornung  et  al.  [13  summarized  the  progress  made  over  a  decade  at  the  T5  reflected  shock 
tunnel  at  Caltech  on  the  subject  of  transition  and  transition  control.  The  first  relevant 
study  was  performed  by  Germain  [llj,  who  investigated  transition  in  hypervelocity  bound¬ 
ary  layers  on  a  5°  slender  cone  instrumented  with  thermocouples.  Heat  flux  profiles  were 
obtained  and  the  transition  location  was  estimated  in  each  case  from  the  point  at  which 
a  rapid  rise  in  heating  was  observed.  The  model  was  tested  in  air  and  N2  with  total  flow 
enthalpies  ranging  from  5  to  15MJ/kg.  The  data  obtained  from  this  study  is  shown  in 
figure  1.1. 

These  results  were  extended  by  Adam  [lj  to  include  experiments  with  CO2  as  the  test 
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Figure  1.1:  Relation  between  transition  Reynolds  number  calculated  at  reference  conditions 
and  stagnation  enthalpy.  Open  symbols  indicate  most  fully  laminar  flows.  Reproduced  from 
[1. 

gas  using  the  same  cone  model.  The  range  of  flow  enthalpies  was  4-10  MJ/kg.  Figure  1.2 
compiles  the  transition  data  from  all  air,  N2  and  CO2  experiments.  The  open  symbols 
correspond  to  cases  in  which  the  flow  was  laminar  to  the  end  of  the  cone.  Figures  1.1 
and  1.2  both  include  data  from  cold  hypersonic  facilities,  obtained  from  DiCristina  |6j  and 
Demetriades  |5  .  In  these  figures  Re^r  is  the  reference  transition  Reynolds  number  based 
on  the  distance  from  the  cone  tip  to  the  transition  location,  measured  along  the  surface  of 
the  cone  (rather  than  the  axial  distance).  The  flow  properties  are  evaluated  at  the  Eckert’s 
reference  temperature  [8  .  Figure  1.2  demonstrates  that  the  value  of  for  CO2  flows  is 
typically  four  times  larger  than  those  for  flows  of  either  air  or  IV2  at  the  same  total  enthalpy. 
This  phenomenon  is  unique  and  serves  as  the  basis  for  this  study. 

These  experimental  trends  were  confirmed  by  Johnson  et  al.  [14J  who  used  a  linear 
stability  analysis  to  compute  Re*j.r  in  CO2,  air,  and  N2  flows.  These  authors  concluded  that 
the  low  dissociation  energies  and  large  number  of  vibrational  modes  of  the  CO2  molecule 
absorb  energy  from  the  acoustic  disturbances  in  the  boundary  layer  and  thus  delay  transi- 
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Figure  1.2:  Relation  between  the  transition  Reynolds  number,  calculated  at  reference  con¬ 
ditions,  and  the  stagnation  enthalpy,  for  air,  N2  and  CO2  shots  in  T5  (cold  tunnel  data  is 
included  for  reference).  Open  symbols  indicate  most  fully  laminar  flows.  Reproduced  from 
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tion.  Figure  1.3  compares  the  amplification  rates  for  reacting  and  nonreacting  disturbances 
in  reacting  mean  flows  of  air  and  CO2.  The  figure  illustrates  that  in  the  case  of  CO2  flows, 
reacting  disturbances  are  damped  considerably  in  comparison  to  nonreacting  disturbances 
for  most  frequencies,  while  for  air  the  disturbances  are  similar  whether  or  not  they  are 
reacting.  Even  though  transition  is  a  non-linear  phenomenon,  Johnson  et  al.  note  that 
the  frequency  of  the  most  amplified  disturbance  found  by  linear  analysis  is  the  same  as  in 
the  non-linear  regime.  Also,  the  linear  amplification  step  is  the  slowest  of  the  steps  leading 
to  transition  [20  .  Therefore,  linear  analysis  can  be  used  to  give  an  order  of  magnitude 
estimate  of  the  growth  rate  and  dominant  frequencies  of  the  acoustic  disturbances  in  the 
boundary  layer. 

Fujii  et  al.  [9;  10  studied  the  effects  of  relaxing  processes  on  transition  in  high  enthalpy 
flows  over  swept  cylinders  (sweep  angles:  45  and  60°).  In  this  work,  the  sound  absorption 
rate  per  wavelength  resulting  from  finite  rate  relaxation  processes  was  computed  for  N2, 
CO2,  and  air  using  an  inviscid  linear  stability  analysis.  These  absorption  rates  were  then 
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Figure  1.3:  Amplification  rates  for  reacting  and  nonreacting  disturbances  in  a  reacting  mean 
flow  of  (left)  air  (Shot  1162)  and  (right)  carbon  dioxide  (Shot  1150).  Figures  reproduced 
from  [14 

compared  with  the  growth  rate  of  wave  disturbances  in  the  boundary  layer.  Figure  1.4 
shows  that  for  CO2,  the  sound  absorption  rates  are  greater  than  the  amplification  rates  for 
acoustic  disturbances  in  the  frequency  range  (1-10  MHz)  where  the  acoustic  disturbances 
are  greatest.  From  this,  we  might  expect  the  relaxation  processes  to  lead  to  significant 
acoustic  damping  in  CO2.  In  comparison,  for  the  case  of  air,  the  sound  absorption  rates 
peak  much  before  the  acoustic  amplification  rates  become  significant  and  nonequilibrium 
processes  would  thus  not  be  expected  to  lead  to  transition  delay.  The  experimental  data  of 
these  authors  showed  good  agreement  with  the  numerical  predictions. 


Absorption  rate  Absorption  rate 
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(a)N2  (b)Air 


Figure  1.4:  Sound  absorption  rate  per  wavelength  as  a  function  of  frequency  for  N2,  air  and 
C02.  Reproduced  from  [10  . 
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Chapter  2 

Experimental  Setup 

2.1  The  GALCIT  T5  Hypervelocity  Shock  Tunnel 

The  facility  used  in  all  experiments  in  the  current  study  was  the  T5  hypervelocity  shock 
tunnel.  It  is  the  fifth  in  a  series  of  free-piston  driven,  reflected  shock  tunnels  built  by  R.J. 
Stalker,  H.G.  Hornung  and  colleagues.  The  T5  facility  consists  of  four  major  components: 
the  secondary  air  reservoir  (2R),  the  compression  tube  (CT),  the  shock  tube  (ST),  and  the 
test  section/dunrp  tank.  The  first  three  of  these  components  are  illustrated  in  figure  2.1. 


Figure  2.1:  Schematic  of  the  T5  hypervelocity  shock  tunnel. 

Two  types  of  experiments  were  performed  during  this  study:  in  the  majority  of  shots 
the  facility  was  run  under  normal  operating  conditions,  but  a  few  shots  were  performed 
at  low-enthalpy  (LE)  conditions,  in  which  T5  was  essentially  run  as  a  large  Ludwieg  tube. 
Under  normal  conditions,  a  heavy  (120  kg)  piston  is  driven  down  the  CT  by  the  release  of 
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high-pressure  air  from  the  2R.  For  the  experiments  in  the  present  series,  the  initial  pressure 
in  the  2R  was  between  800  and  1150  psi,  while  the  light  gas  mixture  in  the  CT,  typically 
85%  He  and  15%  Ar,  was  initially  at  an  absolute  pressure  of  98  to  116  kPa.  The  CT  gas 
is  compressed  adiabatically  by  the  advance  of  the  piston  until  the  pressure  at  the  primary 
diaphragm,  located  at  the  junction  of  the  CT  and  ST,  reaches  a  sufficiently  high  pressure 
that  the  diaphragm  bursts.  In  the  present  shots,  the  primary  diaphragm  consisted  of  a  0.187 
to  0.270-inch  thick  stainless  steel  plate  etched  with  a  X-shaped  groove,  between  0.042”  and 
0.073”  deep,  to  control  the  burst  pressure.  Burst  typically  occurred  at  between  75  and 
110  MPa.  The  test  gas  in  the  ST  was  either  air  or  N2,  at  an  initial  absolute  pressure  of 
between  76  and  117kPa.  The  burst  of  the  primary  diaphragm  produces  a  shock  wave  that 
travels  the  length  of  the  ST  and  reflects  from  the  end  wall,  producing  stagnation  conditions 
at  the  end  of  the  ST,  which  serves  as  reservoir.  The  incident  shock  also  bursts  the  secondary 
diaphragm,  consisting  of  a  0.002”  mylar  membrane  located  at  the  ST-nozzle  junction.  The 
test  gas  then  expands  through  the  nozzle,  flowing  into  the  test  section  and  finally  into  the 
dump  tank.  The  test  section  and  dump  tank  are  initially  evacuated,  separated  from  the  ST 
by  the  secondary  diaphragm.  Startup  of  the  flow  in  the  test  section  typically  takes  1  ms 
from  the  time  of  arrival  of  the  incident  shock  at  the  nozzle  throat;  the  test  time  is  of  the 
order  of  1-2  ms.  In  all  experiments,  a  contoured  nozzle  of  area  ratio  100  was  used. 

Unfortunately,  midway  through  these  experimental  series  the  facility  compressor  broke 
down,  meaning  that  sufficiently  high  2R  pressures  to  run  the  facility  under  normal  condi¬ 
tions  were  unattainable  for  a  period  of  approximately  six  months.  The  decision  was  made 
to  continue  testing  under  low  enthalpy  conditions,  with  the  principal  aim  of  obtaining  visu¬ 
alizations  of  flow  injection.  In  the  LE  configuration,  the  primary  diaphragm  was  removed, 
meaning  the  CT  and  ST  were  at  the  same  initial  pressure,  typically  122  kPa.  The  piston 
was  again  driven  by  a  compressed  gas  from  the  2R,  but  at  a  much  lower  initial  pressure 
than  under  normal  conditions.  For  the  present  experiments,  the  2R  gas  was  N2,  supplied 
by  bottles  (N2  was  found  to  be  more  readily  available  than  air);  three  to  four  bottles  were 
required  to  fill  to  the  desired  pressure  of  170  psi.  The  secondary  diaphragm  was  still  used  in 
this  configuration,  now  serving  as  the  main  diaphragm.  Two  mylar  membranes  were  used 
to  obtain  the  desired  burst  pressure  of  approximately  2  MPa.  Other  details  of  the  operation 
are  the  same  as  in  the  normal  configuration:  see  Kaneshige  [15  for  a  more  complete  descrip¬ 
tion  of  the  LE  configuration.  It  should  be  noted  that  Kaneshige  built  a  special  nozzle  throat 
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for  his  experiments,  whereas  in  the  present  case  the  standard  30  mm  contoured  throat  was 
used. 

The  T5  facility  is  instrumented  with  various  diagnostic  tools.  Two  LVDTs  are  used  to 
measure  the  recoils  of  the  2R  and  CT.  An  accelerometer  is  attached  to  the  CT,  and  was  in 
the  past  used  as  part  of  the  laser-timing  system;  in  the  current  experiments  the  signal  from 
this  accelerometer  was  used  to  trigger  injection.  Pressure  transducers  are  located  at  various 
points  within  the  facility.  Two  are  located  at  the  primary  diaphragm  station,  and  are  used 
to  determine  the  burst  pressure.  Pressure  transducers  placed  along  the  length  of  the  ST 
mark  the  arrival  times  of  the  incident  shock;  these  arrival  times  and  the  known  distance 
between  the  pressure  transducers  are  used  to  determine  the  shock  speed,  which  is  used  to 
calculate  the  stagnation  enthalpy  in  the  reservoir.  Two  further  pressure  transducers  are 
located  just  ahead  of  the  end  wall  of  the  ST,  and  are  used  both  to  determine  the  reservoir 
pressure  of  the  test  gas  and  to  trigger  the  T5  data  acquisition  system. 

The  T5  data  acquisition  system  (DAS)  is  in  two  parts.  The  older  DSP-based  system  is 
used  to  record  facility  data,  as  well  as  providing  approximately  30  channels  for  model  data. 
The  newer  Nl-based  system  provides  an  additional  48  channels.  Both  systems  record  at 
200  MHz. 

The  T5  optical  setup  is  a  typical  Z-arrangement  schlieren  system.  A  high-speed  digital 
camera  allows  the  recording  of  a  sequence  of  images  during  the  test  time;  this  is  used  in 
conjunction  with  either  a  continuous  white  light  source  or  a  spark-gap  source. 

2.2  The  Cone  Model 

The  model  employed  in  the  current  experiments  was  the  slender  cone  used  in  a  number 
of  previous  experimental  studies  in  T5,  including  those  of  Germain  [11,  Adam  [1  ,  and 
Rasheed  [19J.  It  is  a  slender,  5°  half-angle  cone  of  approximately  lrn  in  length  and  is 
composed  of  three  sections:  a  sharp  tip  made  of  molybdenum  (to  withstand  the  high  heat 
fluxes),  a  mid-section,  and  the  main  body.  A  schematic  of  the  model  is  shown  in  figure  2.2. 
Given  the  long,  slender  nature  of  this  body,  it  was  important  to  ensure  that  the  flow  exiting 
the  T5  nozzle  was  parallel,  hence  the  use  of  the  contoured  nozzle.  Also,  to  ensure  that  the 
expansion  wave  emanating  from  the  nozzle  edge  did  not  affect  the  flow  upstream  of  any  of 
the  thermocouples  located  on  the  surface  of  the  model,  the  cone  was  positioned  so  that  the 
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cone  tip  intruded  slightly  into  the  nozzle. 


Main  Body  -  790mm 


Figure  2.2:  The  slender  cone  model  used  in  the  current  experiments. 

For  part  of  the  present  experimental  series,  an  injection  system  was  added  to  the  cone 
model.  The  modifications  to  the  cone  model  itself  included  a  rigid  tube  running  down  the 
centerline,  as  well  as  an  injector  piece  replacing  the  solid  aluminum  mid-section  shown  in 
figure  2.2.  Further  details  of  the  injection  system  are  given  in  section  2.3.  Aside  from  these 
modifications,  the  current  model  configuration  is  that  of  Rasheed,  who  added  a  sleeve  over 
the  cone  to  divide  it  into  smooth  and  porous  half-conical  sections,  each  with  27  thermocou¬ 
ples.  A  layout  of  the  thermocouple  positions  on  the  smooth  side,  together  with  tabulated 
values  of  the  axial  coordinates,  is  shown  in  figure  2.3.  For  the  present  experiments,  only 
the  readings  of  the  thermocouples  on  the  smooth  side  were  of  interest;  unfortunately,  it  was 
found  that  several  of  the  thermocouples  had  sustained  damage,  and  only  21  were  opera¬ 
tional. 

2.3  The  Injection  System 

A  schematic  of  the  injection  system  in  shown  in  figure  2.4,  together  with  a  photograph  of 
the  components  outside  the  test  section.  The  injector  is  fed  from  a  high-pressure  run  tank, 
which  is  filled  from  a  CO2  bottle  (and  then  disconnected)  prior  to  the  shot.  The  tank  is  a 
high-pressure  model  of  approximately  5-gallon  capacity,  and  is  monitored  by  both  a  digital 
pressure  gauge  and  a  Kulite  pressure  transducer,  the  latter  connected  to  the  T5  DAS  to 
allow  recording  of  the  changing  pressure  during  the  tank  discharge.  The  run  tank  leads  to 
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Figure  2.3:  Schematic  diagram  showing  the  locations  of  the  thermocouples  on  a  developed 
view  of  the  smooth  half  of  the  cone  model.  The  azimuthal  lines  are  drawn  at  12.7  mm 
(0.5 in)  intervals  with  the  second  line  located  at  165  mm  (6.5 in).  The  first  line  shows  the 
location  of  the  interface  between  the  cone  tip  and  the  injector/mid-section.  All  distances 
given  are  along  the  cone  surface. 
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the  cone  model  via  a  series  of  flexible  metal  hoses  and  rigid  tubing,  described  in  more  detail 
in  section  2.5,  An  actuated  ball  valve,  controlled  in  turn  by  a  solenoid  valve,  lies  along  the 
flow  path  and  controls  the  timing  of  the  injection  initiation.  In  experiments,  the  system 
was  triggered  from  a  signal  generated  by  an  accelerometer  attached  to  the  CT,  allowing 
sufficient  time  for  the  injection  ffow  to  be  established  before  the  arrival  of  the  main  flow  in 
the  test  section. 


Figure  2.4:  Schematic  and  photograph  of  the  injection  system. 


Drawings  of  the  final  part  of  the  flow  path  within  the  cone  model  can  be  found  in  fig¬ 
ure  2.5.  This  was  an  earlier  design,  however,  and  the  configuration  was  modified  partway 
through  testing,  so  this  should  only  be  used  as  a  rough  guide. 1  All  injector  tips  were  de¬ 
signed  such  that  the  flow  would  choke  at  the  injector  holes.  However,  due  to  an  error  in 
manufacturing,  the  holes  in  the  tip  used  in  all  experiments  described  here  (tip  #1)  were 
of  larger  diameter  than  designed.  The  minimum  area  thus  occurred  at  the  1/4-inch  tubing 
leading  to  the  plenum  from  which  the  injector  holes  were  fed  (see  figure  2.5).  It  is  thus 
assumed  in  mass  flow  calculations  that  the  sonic  point  was  located  at  the  end  of  this  tub¬ 
ing,  immediately  before  the  plenum.  The  injector  tip,  containing  the  plenum  and  injection 
holes,  is  described  in  the  following  section. 


1A11  Solidworks  drawings  shown  in  this  chapter  were  produced  by  Bahrain  Valiferdowsi. 
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Figure  2.5:  Schematic  of  the  injection  components  within  the  cone  model:  (top)  the  up¬ 
stream  end  of  the  cone,  which  holds  the  cone  tip  injector;  (below)  the  downstream  end  of 
the  cone,  showing  the  flange  which  holds  the  rigid  metal  tube  for  transporting  the  CO2. 
Note  that  the  test  configuration  differed  slightly  from  that  shown. 


14 


2.3.1  Injector  tip  design 

The  injector  tip  took  the  place  of  the  cone  mid-section  piece  shown  in  figure  2.2,  and,  as 
such,  the  design  of  this  piece  was  limited  by  the  constraints  that  this  external  geometry 
placed  upon  it.  Even  within  these  constraints,  however,  there  were  many  factors  that  needed 
to  be  considered  within  the  design.  A  review  of  the  literature  revealed  that  relevant  work 
has  focused  principally  on  injecting  gases  onto  turbine  blades  for  cooling  or  film  cooling,  and 
no  reference  could  be  found  to  injecting  gases  into  hypervelocity  boundary  layers.  Design 
inspiration  was  thus  taken  from  previous  work  in  film  cooling  [4  . 

Thin-film  cooling  is  employed  in  gas  turbines  because,  by  cooling  the  blades,  it  is  possible 
to  run  the  turbine  at  higher  temperature  and,  therefore,  higher  efficiency.  In  thin  film 
cooling,  the  object  is  to  inject  a  thin,  uniform  layer  of  gas  that  covers  as  much  of  the  blade 
as  possible  for  uniform  cooling.  While  the  current  application  also  requires  a  thin  layer  of 
gas,  the  needs  are  somewhat  different.  Instead  of  a  complete  blanket  of  CO2  adhering  to  the 
cone  surface,  the  requirement  is  for  the  gas  to  mix  thoroughly  into  the  test  gas  boundary 
layer  as  quickly  as  possible  in  order  that  the  temperature  of  the  injected  CO2  (initially  at 
room  temperature,  «300  K)  is  raised  to  a  value  close  to  that  of  the  test  gas  in  the  boundary 
layer  (~2000K),  thereby  exciting  the  vibrational  modes  and  causing  dissociation.  Only 
under  these  circumstances  can  the  injected  gas  be  expected  to  have  a  transition-suppressing 
effect. 

The  first  design  decision  was  whether  to  employ  holes  or  slots.  Early  efforts  in  thin-film 
blade  cooling  concentrated  on  slots,  as  these  produced  2-D  sheet-like  jets  which  were  well 
suited  to  the  requirements  of  this  application.  In  the  current  application,  however,  the  use 
of  holes  seemed  to  hold  several  key  advantages,  both  in  the  ease  of  manufacturing  and  in 
that  the  injected  gas  in  this  case  would  be  mixed  more  rapidly  into  the  existing  boundary 
layer. 

Having  decided  upon  holes,  the  relevant  design  parameters  were  the  size  and  arrange¬ 
ment  of  holes  on  the  injector  surface,  and  the  hole  shaping  and  exit  angle.  The  aim  was 
to  produce  a  minimal  disturbance  to  the  boundary  layer  while  injecting  a  sufficient  mass 
flow  of  CO2,  with  adequate  coverage,  for  the  transition-delaying  effects  to  become  apparent. 
The  standard  in  thin-film  cooling  of  turbine  blades  is  to  design  the  holes  with  as  shallow  an 
exit  angle  as  possible  and  with  a  spacing  of  at  least  two  hole  diameters  between  adjacent 
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holes  [4  .  The  present  design  follows  this  convention.  The  geometry  of  the  injector  tip 
dictated  that  the  minimum  exit  angle  for  a  straight  hole  tunnel  was  11°  for  one  row  of  holes 
or  12°  for  two.  The  circumferential  spacing  was  at  least  2  hole  diameters  in  all  cases;  for 
configurations  with  two  rows  of  holes,  the  axial  spacing  between  rows  was  2.5  to  5  hole 
diameters,  and  the  rows  were  staggered  with  respect  to  one  another. 

A  further  consideration  was  the  ratio  of  the  length  of  the  hole-tunnel  to  the  hole  diam¬ 
eter,  L/D.  By  minimizing  this  value,  total  pressure  losses  due  to  friction  would  be  kept  to  a 
minimum.  In  the  literature,  the  value  of  L/D  typically  ranges  from  1  to  10  |4|.  To  minimize 
the  injection  angle  while  ensuring  adequate  strength  of  the  injector  tip,  a  somewhat  larger 
value  was  necessary  in  the  current  configuration.  The  tunnel  length  was  typically  25  mm 
while  the  hole  diameter  was  0.5  mm,  resulting  in  a  L/D  ratio  of  approximately  50. 

Recently,  a  major  advance  in  thin-film  cooling  has  been  the  incorporation  of  shaping  into 
the  exit  hole  to  improve  blade  surface  coverage  [4  .  Such  shaping  has  not  been  adopted  in 
the  current  configuration,  but  could  possibly  be  used  in  the  future  to  improve  the  tangency 
of  the  flow  (axial  shaping  is  envisaged,  rather  than  the  lateral  fanning  that  is  typically 
employed  to  improve  injection  coverage). 

The  position  of  injection  was  also  a  compromise.  Injection  near  the  cone  tip  would 
provide  increased  time  for  the  injected  gas  to  mix  with  the  test  gas  and  reach  nonequilib¬ 
rium  conditions.  Injecting  further  downstream,  where  the  boundary  layer  is  thicker,  would 
decrease  the  likelihood  of  the  injected  flow  blowing  off  the  boundary  layer.  In  the  end, 
the  former  of  these  considerations,  together  with  the  geometry  constraints  provided  by  the 
existing  model,  led  to  the  selection  of  an  injection  position  near  the  cone  tip. 

The  injector  tip  was  designed  such  that  the  injected  flow  would  choke  at  the  exit  holes,  in 
order  that  the  injection  conditions  could  be  easily  calculated.  The  injection  flow  rate  would 
be  controlled  through  adjusting  the  initial  pressure  in  the  CO2  run  tank.  Approximate 
calculations  of  the  required  and  supplied  mass  flow  rates  in  the  present  configuration  are 
provided  in  section  2.5. 

Because  of  the  uncertainty  in  both  the  optimal  design  for  the  injector  tip  and  the  optimal 
manufacturing  method,  a  number  of  tips  were  constructed,  with  the  intention  that  several 
of  these  would  be  tested  to  determine  the  best  compromise.  The  important  parameters  of 
these  tips  are  listed  in  table  2.1. 

The  original  mid-section  was  made  from  aluminum.  The  first  version  of  the  cone  tip 
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Table  2.1:  Injector  tip  parameters. 


Injector 
tip  no. 

manufac¬ 
tured  by: 

avg  hole  size 
[inch] 

number 
of  rows 

angle 

[deg] 

image  comments 

nos. 

1 

Solid  Cone 

0.0258 

2 

12 

insert 

lb 

Scicon 

blocked 

1 

12 

solid  metal 

2 

Scicon 

0.0163 

1 

12 

nickel  plated;  insert 

3 

Solid  Cone 

0.0310 

2 

12 

4 

Solid  Cone 

blocked 

2 

12 

5 

Scicon 

0.0350 

1 

12 

nickel  plated 

6 

Scicon 

0.0271 

1 

11 

nickel  plated 

7 

Scicon 

0.0171 

2 

12 

8 

Solid  Cone 

0.0301 

1 

11 

white 

injector  was  made  by  rapid  prototyping  a  high  impact  plastic.  The  second  version  was  also 
made  by  rapid  prototyping,  but  a  metal  powder  was  used.  The  second  tip  was  discarded 
since  the  holes  could  not  be  resolved  and  the  surface  was  extremely  rough.  The  added  cost 
of  drilling  holes  at  an  angle  and  machining  the  surface  to  give  it  a  smooth  finish  made 
it  economically  unviable.  The  third  version  was  made  again  by  rapid  prototyping  a  high 
impact  plastic,  but  it  was  then  nickel  plated  to  protect  it  from  the  high  heat  flux.  The  metal 
coating  was  on  the  order  of  10  gm,  but  since  the  maximum  temperature  rise  experienced 
by  the  thermocouples  is  30  K  for  a  test  time  of  4  ms,  it  was  calculated  that  the  penetration 
depth  would  be  sufficiently  small  that  this  metal  coating  would  be  adequate  to  protect  the 
plastic.  Visual  inspection  of  the  coating  after  a  high  enthalpy  shot  confirmed  this. 

The  injector  tip  (#1)  that  was  used  in  all  relevant  experiments  described  in  this  study 
(i.e.  series  4  and  5;  for  other  series  the  original  solid  midsection  was  used)  is  shown  in 
figure  2.6.  This  tip  was  expected  to  provide  injection  conditions  closest  to  those  desired, 
so,  given  the  limited  testing  time,  only  this  tip  was  tested.  The  internal  geometry  of  the 
tip  may  be  seen  in  the  drawings:  a  short  section  of  pipe  opens  into  a  hollowed  out  plenum, 
that  leads  to  the  surface  of  the  cone  by  way  of  two  rows  of  36  holes,  spaced  10°  apart  with 
the  rows  staggered  from  each  other  by  5°.  The  hole  path  is  at  a  12°  angle  (in  the  same 
sense  as  the  half-angle  of  the  cone)  and  each  hole  has  a  diameter  of  0.5  mm. 

2.3.2  Injection  timing 

Controlling  the  timing  of  the  injection  was  an  important  consideration,  as  there  was  a 
significant  startup  time  from  triggering  to  the  establishment  of  quasi-steady  injected  flow 
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Figure  2.6:  Schematic  and  drawings  of  the  cone  tip  injector  piece. 


at  the  surface  of  the  cone.  Initiation  of  the  injection  was  controlled  by  an  actuated  ball 
valve,  located  a  small  distance  down  the  flow  line  from  the  run  tank.  This  ball  valve  was 
in  turn  controlled  by  a  solenoid  valve.  Both  valves  normally  resided  in  the  closed  position: 
when  the  solenoid  valve  received  a  12-volt  trigger  signal,  it  opened  to  allow  compressed  air 
to  pass  through  and  actuate  the  ball  valve.  The  solenoid  stayed  opened  until  the  12-volt 
signal  ceased;  when  it  closed,  the  compressed  air  supply  was  cut  off  and  the  spring  action 
of  the  ball  valve  returned  it  to  its  original  position.  Because  this  spring  action  opposed  the 
opening  of  the  valve,  the  opening  time  was  somewhat  longer  than  the  closing  time,  and  there 
was  some  concern  that  the  ball  valve  would  not  open  fast  enough  for  the  CO2  to  reach  the 
surface  of  the  cone  by  the  onset  of  the  test  flow.  For  this  reason,  it  was  decided  to  trigger 
the  injection  by  the  earliest  signal  produced  during  the  shot:  that  of  the  accelerometer 
attached  to  the  CT. 

Several  electronic  components  were  designed  as  part  of  the  triggering  system:  these  are 
described  in  the  following  subsections. 
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2. 3. 2.1  Operational  amplifier 


The  signal  coming  from  the  accelerometer  was  very  weak  (around  40  mV),  so  it  was  nec¬ 
essary  to  amplify  this  signal.  Originally  it  was  input  directly  to  the  operational  amplifier 
circuit  shown  in  figure  2.7,  However,  there  was  an  impedance  mismatch  which  altered  the 
input  signal.  It  was  therefore  necessary  to  bring  the  accelerometer  signal  to  a  delay  genera¬ 
tor  and  use  the  TTL  output  from  that  to  input  to  the  op  amp.  The  delay  generator  was  set 
at  zero  delay  and  the  ~  3  volt  TTL  output  pulse  was  then  amplified  to  trigger  the  timing 
circuit  shown  in  figure  2.8.  This  op  amp  doubles  the  TTL  amplitude,  which  is  more  than 
the  voltage  necessary  to  trigger  the  timing  circuit. 

A  large  (1 MQ)  resistor,  R3,  was  placed  to  draw  down  the  voltage  once  the  pulse  ended, 
otherwise  the  op  amp  would  continue  to  output  the  amplified  TTL  pulse  until  the  power 
supply  was  switched  off.  This  resistor  can  be  seen  in  figure  2.7, 


Op  Amp 


Notes: 

TC=timing  circuit  ; 

Vs  -  +ve  of  9  V  power  supply  8  V  power  supply  goes  to  optoelectronic  interrupters  (plug-in  unit  Jameco) 

GRD  -  -ve  of  9V  power  supply  28  V  power  supply  for  kulite  -  mounted  to  DAS  (OMROM) 


Figure  2.7:  Operational  amplifier. 


2. 3. 2. 2  Timing  circuit 

The  timing  circuit  shown  in  figure  2.8  takes  the  output  from  the  op  amp,  in  this  case  the 
amplified  TTL  pulse,  generates  a  delay  and  then  outputs  a  pulse  of  amplitude  ~6  volts 
and  the  desired  duration.  The  delay  is  determined  by  the  time  constant  t±  =  R,\  x  C\. 
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In  this  case,  R\  =  4kO  and  C\  =  0.1  fi F,  so  the  delay  was  less  than  0.4  ms.  The  output 
pulse  duration  was  determined  by  the  time  constant  72  =  i?2  x  C2.  Since  R2  =  500  kO  and 
C2  =  22  fiF,  the  pulse  should  have  lasted  11s,  but  in  practice  was  found  to  last  just  over 
7  s.  This  was  more  than  sufficient,  however,  given  the  timescales  of  the  experiment. 


Timing  Circuit 


GRD  -  -ve  of  1 2V  power  supply 

R1-C1  control  delay  time 

R2-C2  control  time  circuit  stays  open 


Figure  2.8:  Timing  circuit  diagram. 


2.3. 2.3  Relay 

The  purpose  of  the  relay  in  figure  2.9  is  to  close  the  circuit  between  the  power  supply  and 
the  solenoid  valve  which  controls  the  flow  of  compressed  air  that  actuates  the  ball  valve. 
When  the  solenoid  valve  opens,  the  compressed  air  rushes  through  and  turns  the  ball  valve 
by  90  deg.  Once  the  relay  stops  receiving  an  input  voltage,  it  breaks  the  circuit,  causing 
the  solenoid  valve  to  close,  thereby  cutting  off  the  compressed  air  supply  to  the  actuated 
ball  valve.  The  ball  valve  is  controlled  by  a  quick  return  actuator,  which  uses  an  internal 
spring  to  rotate  the  ball  valve  by  —90  deg  and  close  the  supply  from  the  run  tank. 

2.4  Diffusion  calculation 

One  concern  is  that  the  injected  CO2  would  not  fully  mix  with  the  air  boundary  layer  over 
the  length  of  the  cone.  Considering  only  mass  diffusion  and  no  convective  or  vorticity-driven 
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Relay 


-ve  solenoid 

-ve  power  supply  12V 


12V 

suP^y  Solenoid  +ve 


Timing  Circuit  (TC)  +ve 
Vcc 


© 


AAVMMn  p 

R  R  1 


R=1kQ 


TC  GRD  +Op  Amp  GRD 
Relay  output  pin  3 


Notes: 

12V  power  supply  is  mounted  on  chassis  and  powers  the 

breadboard  and  the  relay/solenoid 

Relay  has  same  ground  as  the  breadboard 


Figure  2.9:  Relay  circuit  diagram. 

mixing,  it  can  be  shown  that  the  two  gases  would  fully  mix  before  the  end  of  the  cone.  The 
effect  of  these  other  mechanisms  would  further  aid  the  mixing  of  the  two  gases. 

The  following  diffusion  calculation  follows  Anderson|3j  pages  696-699.  Based  on  kinetic 
theory,  the  general  diffusion  equation  is 


.  Vf* 

Dab  =  K'd  — 

where  Dab  is  the  binary  diffusion  coefficient  (or  diffusivity)  for  species  A  into  B,  cr  is 
the  collision  cross  section,  and  K'D  is  a  constant|3j. 

More  specifically,  from  Fick’s  law  (equation  2.1  3;  7  ) ,  Dab  is  given  by 


JA 


Dab 


-p  Dab^ca 


0.0018583 


Pd2AB£lc>,AB 


(2.1) 

(2.2) 


where  j a  is  the  mass  flux  of  species  A  in  mass  per  second  per  unit  area,  ca  is  the 
mass  fraction  of  A  and  Mi  is  the  molecular  weight  of  species  i|3  .  Dab  is  given  in  units  of 
[■ cm2 / s\ ,  T  in  [K],  P  in  [atm],  and  cLab  in  [A] [3].  Refer  to  Table  16.1  [3  for  tabulated  values 
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of  Qd,ab  as  a  function  of  k\T/EAB ■  Here,  e/k\  and  d  are  the  Lennard-Jones  parameters 
which  represent  the  characteristic  energy  of  interaction  and  the  characteristic  molecular 
diameter  respectively [3  .  eab  and  cIab  can  be  approximated  by 

£AB  =  V£a£b 
dAB  ~  ^  A  +  ^ B ) 


From  Table  16.2 [3j , 


Table  2.2:  Parameters  used  in  diffusion  computations,  see  Anderson  |3j,  Table  16.2 
Gas  Molecular  Weight  Lennard-Jones  Parameters 
M  a  e/ki 

[g/rnol]  [A]  [K] 


air 

28.97 

3.617 

97.0 

n2 

28.02 

3.681 

91.5 

co2 

44.01 

3.996 

190 

eab  =  j£_A£B  =  ^97.0  X  190  =  135.8 
k\  y  k\  k\ 

dAB  ~  ^  (3.617  +  3.996)  =  3.807 

Assume  that  the  temperature  is  0.7  of  the  stagnation  temperature,  which  is  typically 
around  3000  K.  A  pressure  of  0.25  atm  is  typical  at  the  edge  of  the  boundary  layer.  Using 
T=2000K, 


T  2000 

- —  = - =  14.73 

£AB/ki  135.8 

This  is  used  to  find  the  value  for  Qd,ab  in  Table  2.3  3  .  Values  are  given  for  Tki/s  = 
10  and  20.  Interpolating  between  these  values, 

Substituting  these  values  into  Equation  2.2,  the  diffusivity  is  then 


Dab 


0.0018583- 


20003 


28.97 


+ 


1 


44.01 


) 


0.25  x  (3.8072)  x  0.7056 
15.55  cm2 js  =  0.001555  m2 / s 
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Table  2.3:  Additional  parameters  used  in  diffusion  calculations  [3 


ki  T/e 
or 

kiT/eAB 

Qfj,  —  k 
for  viscosity  and 
thermal  conductivity 

^ D,AB 
for  mass 
diffusivity 

10 

0.8242 

0.7424 

14.73 

0.7859 

0.7056 

20 

0.7432 

0.6640 

Taking  a  reference  speed  of  Uref  =  4000  m/s  at  half  the  expected  boundary  layer  thick¬ 
ness,  the  distance  of  interest,  tBL,  is  0.0005  m.  The  diffusion  distance,  Xdiff,  is  then 


xdiff 


t2BL  X  Uref  _  0.00052  X  4000 


=  0.643  m 


Dab  0.001555 

Thus,  through  diffusion  alone,  the  injected  CO2  should  be  fully  mixed  by  64  cm  along 
the  surface  of  the  cone.  Since  there  will  be  convective  mixing  in  addition  to  mass  diffusion, 
the  two  gases  should  be  fully  mixed  around  halfway  along  the  surface  of  the  cone. 


2.5  Injection  mass  flux  calculations 

In  this  section,  estimates  are  made  of  both  the  injection  flow  rate  necessary  to  supply  an 
adequate  amount  of  CO2  to  the  cone  boundary  layer  and  the  mass  flow  rate  provided  by  the 
current  injection  system.  The  latter  is  based  on  the  geometry  of  the  injection  components 
and  the  known  pressure  in  the  run  tank. 

2.5.1  Required  injection  mass  flux  in  the  cone  boundary  layer 

The  following  simple  calculation  gives  an  indication  of  the  injected  CO2  mass  flow  rate 
necessary  to  produce  a  significant  effect  on  the  transition  Reynolds  number. 

The  target  condition  for  the  injection  experiments  is  that  of  Shot  675,  for  which  it  has 
been  calculated  [lj  that  the  flow  velocity  and  density  at  the  edge  of  the  boundary  layer, 
ue  and  pe ,  are  4208 m/s  and  0.0786 kg/rn3,  respectively.  Within  the  boundary  layer,  the 
velocity  will  decrease  towards  the  cone  surface  while  the  density  will  increase  (as  the  cone 
temperature  is  lower  than  that  of  the  flow),  but  to  a  first  approximation,  we  assume  that 
their  product  remains  constant.  According  to  the  boundary  layer  calculations  of  Adam  [1, 
at  a  cone  diameter  of  dcone= 5  cm  (roughly  10  cm  downstream  of  the  injection  point),  the 
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boundary  layer  thickness,  5,  will  be  approximately  1  mm. According  to  these  assumptions 
then,  the  mass  flux  within  the  boundary  layer  at  this  position  will  be  given  by: 

hife/  =  uepe  —  [(dcone  +  2(5)  dcone]  (2-3) 

=  0.053  kg/s.  (2.4) 

Since  the  amount  of  injected  CO2  that  would  be  necessary  to  lead  to  a  significant  delay  in 
transition  is  not  certain,  let  us  simply  assume  that  an  mass  flow  rate  equal  to  the  boundary 
layer  mass  flux  at  this  point  would  be  required.  An  injection  mass  flow  rate  of  0.053  kg/s  is 
certainly  feasible  with  the  current  injection  system.  If  we  further  assume  that  the  transition 
Reynolds  number  scales  linearly  with  the  fraction  of  CO2  in  the  boundary  layer,  it  is  possible 
to  estimate  the  expected  Re^r  in  this  case.  From  previous  experiments  [13J,  the  transition 
Reynolds  number  for  /io~10MJ/kg  increases  from  ~lxl06  for  100%  air  to  ~5xl06  for 
100%  CO2.  For  equal  amounts  of  air  and  CO2  then,  a  value  of  Rej!r~3xl06  is  obtained. 

2.5.2  Description  of  injection  system  components 

To  estimate  the  mass  flow  rate  in  the  injection  system,  a  complete  description  of  the  physical 
components  of  the  system  is  first  necessary.  The  components  were  as  follows.  Leading  from 
the  run  tank  to  the  ball  valve  is  ~5”  long,  1/2-inch  (0.4”  ID)  rigid  pipe.  From  the  valve  to 
the  feed-through  plate  is  a  90°  bend  fitting,  of  0.4”  ID,  with  a  total  length  of  ~4”,  followed 
by  a  10”  long  1/2-inch  rigid  pipe  (0.4”  ID)  connected  to  a  36”  long,  1/4-inch  flexible  metal 
hose.  According  to  the  manufacturer  (Swagelok),  this  hose  has  a  minimum  internal  diameter 
of  0.19”.  From  the  feed-through  plate  to  the  model,  there  is  first  a  90°  bend  fitting,  of  0.4” 
ID,  with  a  total  length  of  ~5”,  followed  by  a  48”  long  1/2-inch  flexible  metal  hose,  with  a 
quoted  minimum  ID  of  0.4”.  Attached  to  either  end  of  this  hose  is  a  3  1/2”  long,  0.4”  ID 
fitting.  The  hose  leads  to  a  23”  long,  0.4”  ID  rigid  tube,  situated  inside  the  model,  which 
contracts  to  a  6”  long,  1/4-inch  (0.17”  ID)  rigid  section  leading  to  the  injector  plenum. 

For  the  calculation,  the  flow  path  is  modelled  as  a  series  of  straight  pipes  and  area 
changes.  The  flow  through  each  pipe  section  is  assumed  to  follow  the  Fanno  relations,  while 
the  flow  through  the  area  changes  is  assumed  to  take  place  ander  isentropic  conditions. 
Such  a  treatment  ignores  the  total  pressure  losses  induced  at  both  the  ball  valve  and  the 
two  90°  bends  in  the  flow  path,  but  these  are  expected  to  be  small  relative  to  the  losses 
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caused  by  friction  in  the  large  L/D  pipe  sections  in  the  current  configuration. 


2.5.3  The  Fanno  flow  relations 


A  discussion  of  Fanno  flow  may  be  found  in  almost  any  textbook  on  compressible  flow:  the 
summary  in  this  section  follows  Shapiro  [22  .  Fanno  flow  refers  to  the  steady  adiabatic  flow 
of  fluid  through  a  constant-area  pipe  with  friction.  The  frictional  effects  are  characterized 
through  a  coefficient  of  friction,  /,  defined  as 
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7~w 

pu2  / 2 ’ 


(2.5) 


where  tw  is  the  wall  shearing  stress.  It  may  be  shown  that  in  Fanno  flow,  whether  the  flow 
is  subsonic  or  supersonic,  the  Mach  number  always  tends  to  unity.  If  the  length  of  the  pipe 
is  larger  than  a  critical  length,  Lmax/D  (here  normalized  by  the  pipe  diameter),  at  which 
M= 1,  choking  will  occur.  For  a  perfect  gas,  this  critical  length  can  be  related  to  the  initial 
Mach  number,  M,  as 


jLmax  _  1  -  Af 2  7+1  / 7  +  1  M2  \ 

'  D  7A42  +  27  °g  ^  2  l  +  IxlM2  J  ’ 


(2.6) 


where  /  is  the  mean  friction  coefficient  with  respect  to  pipe  length.  Other  flow  quantities 
may  also  be  obtained,  and  are  usually  presented  in  terms  of  their  values  at  the  point  at 
which  the  Mach  number  reaches  unity,  denote  by  a  star.  Of  interest  here  are  the  density, 
velocity,  temperature,  and  sound  speed: 
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The  value  of  the  friction  coefficient  depends  on  a  number  of  factors,  in  particular  whether 
the  flow  is  laminar  or  turbulent,  the  Reynolds  number  of  the  flow,  and  the  roughness  of 
the  pipe.  A  representative  value  for  high  Reynolds  number  flows  is  0.0025.  A  complete 
discussion  of  the  behaviour  of  the  friction  coefficient  under  different  ffow  regimes  may  be 
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found  in  Schlicting  [21  ,  who  uses  a  slightly  different  parameter,  A,  defined  as 
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7~w 

pu2/8' 


(2.9) 


It  follows  immediately  that  A=4/.  For  smooth  pipes,  as  will  be  assumed  the  case  here,  a 
relation  has  been  derived  from  experimental  data  with  Reynolds  numbers  of  up  to  3.4xl06, 
as  follows: 
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2.01og(Re\/A)  -0.8. 


(2.10) 


this  relation  is  known  as  Prandtl’s  universal  law  of  friction  for  smooth  pipes. 


2.5.4  Area  changes 


Flow  between  area  changes  is  assumed  to  take  place  under  isentropic  conditions,  so  the 
standard  isentropic  flow  relations  hold: 


as  well  as: 
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The  respective  quantities  here  are  referenced  either  to  the  sonic  values,  denoted  by  an 
asterisk,  or  the  stagnation  conditions,  denoted  by  a  zero  subscript. 


2.5.5  Initial  calculation 

An  initial  calculation  was  made  using  a  constant  value  of  /=0.0025  in  equation  2.6.  The 
flow  was  assumed  to  choke  at  the  end  of  the  1/4-inch  rigid  pipe  leading  into  the  injector 
plenum,  meaning  there  were  8  locations  at  which  the  flow  conditions  had  to  be  calculated. 
Following  the  flow  path  from  the  run  tank  to  the  injector  plenum,  with  the  internal  geometry 
between  each  consecutive  pair  indicated,  these  were  as  follows: 

1.  Downstream  junction  of  the  run  tank  and  1/2-inch  pipe 
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— +  0.4”ID  pipe  of  length  19”,  L/D= 47.5 

2.  Upstream  junction  of  the  1/2-inch  pipe  and  1/4-inch  hosing 
— >  Area  change,  ^3/^2=0.2256 

3.  Downstream  junction  of  the  1 /2-inch  pipe  and  1/4-inch  hosing 
— »  0.19”ID  hose  of  length  36”,  L/D= 189 

4.  Upstream  junction  of  the  1/4-inch  hosing  and  1/2-inch  pipe 
— >  Area  change,  A^/ A^=AA?>2 

5.  Downstream  junction  of  the  1 /4-inch  hosing  and  1 /2-inch  pipe 
— >■  0.4”ID  pipe/hose  of  length  83”,  L/D= 208 

6.  Upstream  junction  of  the  1/2-inch  pipe  and  1/4-inch  pipe 
— >  Area  change,  A7/Aq=0.214 

7.  Downstream  junction  of  the  1/2-inch  pipe  and  1/4-inch  pipe 
— »  0.185”  pipe  of  length  6”,  L/D= 32 

8.  Exit  of  1/4- inch  pipe  into  plenum 

Using  the  assumed  condition  of  M= 1  at  location  8,  equations  2.6  and  2.11  were  used  to 
calculate  the  Mach  number  at  every  location  back  through  the  system.  The  expansion  in  the 
run  tank  takes  place  under  adiabatic  conditions,  with  no  work  done  either  on  or  by  the  gas, 
so  this  expansion  will  to  a  good  approximation  be  isothermal.  Using  the  known  reservoir 
conditions  and  equations  2.7,  2.8,  2.12,  and  2.13,  the  other  flow  quantities  of  interest  can 
then  traced  back  forward  through  the  system  to  the  location  of  choking. 

Table  2.4  lists  the  flow  quantities  of  interest  calculated  by  this  method.  The  density 
depends  on  the  run  tank  pressure  and  so  is  expressed  at  each  station  as  a  fraction  of  the 
run  tank  density  (note  that  po  in  the  relevant  column  refers  to  this  run  tank  density  rather 
than  the  stagnation  value  at  that  point  in  the  system);  the  other  flow  quantities  may  be 
considered  independent  of  the  run  tank  pressure.  Values  of  7=1.3  and  77=189  J/kgK  were 
assumed  for  CO2. 

Using  the  conditions  at  location  1  (at  which  the  area  is  8.11x10“ 5m2),  the  mass  flow 
rate,  m=pV A,  can  be  easily  calculated.  Given  that  the  quantities  Tq  and  po  are  known 
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Location 

M 

T(K) 

a  (nr/s) 

V  (nr/s) 

p/po 

0 

0 

298 

271 

0 

1 

1 

0.083 

298 

270 

22 

0.997 

2 

0.083 

298 

270 

22 

0.995 

3 

0.402 

291 

267 

108 

0.922 

4 

0.588 

283 

264 

155 

0.639 

5 

0.110 

298 

270 

30 

0.751 

6 

0.112 

298 

270 

30 

0.738 

7 

0.663 

280 

262 

174 

0.600 

8 

1 

259 

252 

252 

0.410 

Table  2.4:  Flow  conditions  at  various  points  along  the  injection  system 


directly,  and  that  T\k,Tq  and  pi~po5  it  is  useful  to  write  the  mass  flux  as 


m 


PaMlAl{ik- 


(2.14) 


Then,  if  po=300psi  (=2.07  x  106Pa),  a  mass  flux  of  m~0. 067  kg/m3  is  obtained. 

Using  table  2.4  and  the  known  system  geometry,  we  may  also  estimate  the  time  for  the 
fluid  to  pass  from  the  run  tank  through  to  the  injector  plenum  as  0.098  s,  meaning  that 
the  total  flow  time  from  run  tank  to  injection  tips,  once  the  flow  is  fully  established,  is  of 
the  order  of  0.1s.  This  compares  with  the  time  taken  for  a  sound  wave  to  propagate  this 
distance  of  approximately  0.014  s  (under  stationary  flow  conditions). 


2.5.6  Refined  calculation 

In  the  previous  calculation,  a  friction  coefficient  of  0.0025  was  assumed  in  all  pipe  sections. 
Using  the  results  of  this  calculation,  together  with  equation  2.10,  a  more  accurate  deter¬ 
mination  of  this  friction  coefficient  is  possible.  Assuming  a  run  tank  pressure  of  ~300psi, 
and  thus  a  density  of  po~36kg/m3,  together  with  a  viscosity  of  /i=1.38xl0-5  (suitable  for 
CO2  at  273  K),  an  average  Reynolds  number  between  locations  7  and  8  of  approximately 
1.2xl06  is  obtained.  Solving  equation  2.10  graphically  gives  A=2.25xl0-3,  from  which 
/=0.56xl0~3.  Repeating  the  calculation  from  the  previous  subsection  with  this  value  of 
/  throughout,  a  Mach  number  at  location  1  of  Mi=0.1085  results.  Substituting  this  into 
equation  2.14  with  the  other  values  as  before  (po=300psi)  now  yields  m~0.087kg/m3. 

This  calculation  may  be  tested  in  the  following  manner.  According  to  the  above  as¬ 
sumptions,  the  run  tank  temperature,  To,  is  constant,  and  so  differentiating  the  ideal  gas 
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equation,  po= PqRTo,  the  rate  of  change  of  the  run  tank  pressure  is 

dPO  nrr,  dPo  ,0  , 

IT  =  RT°W  (2'15) 

As  dpo / dt=—m/Vo ,  where  Vo  is  the  run  tank  volume,  we  thus  have 

dpo  m 

i  =  <2-i6> 

~  -Mi  x/rRToirrPo,  (2.17) 

ho 

where  we  have  made  use  of  equation  2.14.  This  differential  equation  may  be  easily  solved 
to  yield: 

Po(t)  =  po(t0)  exp  {-Mi  spyKTo-±(t  -  t0)).  (2.18) 

ho 

We  thus  see  that  the  characteristic  time,  r,  for  the  run  tank  pressure  drop,  i.e.,  that  time 
required  to  drop  to  1/e  of  its  original  value,  is  (Mi yj' jRTq  - ) ~ 1 .  For  Vb=0.0189m3  (5 
gallons),  Mi=0.1085,  and  the  other  quantities  as  before,  r=7.95s. 

In  figure  2.10,  the  measured  run  tank  pressure  from  shot  2448  is  plotted  together  with 
the  theoretical  profile  given  by  equation  2.18.  The  run  tank  pressure  transducer  signal 
was  recorded  on  a  channel  of  each  of  the  two  T5  data  acquisition  systems:  on  the  new  NI 
system,  the  signal  was  recorded  over  a  period  of  5  s,  and  thus  shows  an  appreciable  drop 
in  jjq:  on  the  old  DSP  system,  only  a  brief  20  ms  window  was  recorded,  and  this  reading 
is  thus  only  useful  for  determing  the  tank  pressure  during  the  test  time.  Actually,  as  may 
be  seen,  there  is  some  difference  in  the  two  readings,  with  the  DSP  signal  indicating  a 
pressure  approximately  lOpsi  higher.  In  fact,  the  initial  pressure  reading  from  the  digital 
pressure  gage  was  also  different  from  the  initial  Kulite  transducer  reading,  at  229  psi.  These 
discrepancies  were  not  investigated  further,  but  do  add  a  further  degree  of  uncertainty  to 
the  mass  flow  rate  calculation. 

The  longer-term  measured  run  tank  profile  shows  a  sharp  initial  drop  in  pressure  as  the 
injection  flow  starts  up.  This  startup  period  appears  to  extend  slightly  past  t=0  and  then 
a  relatively  steady  drop  in  pressure  is  established.  A  further  examination  of  the  injection 
startup  period  is  made  in  section  2.6.  The  reference  point,  to,  in  equation  2.18  is  chosen 
as  fo=50ms,  and  the  resulting  theoretical  curve  is  shown.  The  theoretical  and  measured 
profiles  are  almost  indistinguishable  for  a  period  of  approximately  2  s,  sugggesting  that  the 
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t(s) 


Figure  2.10:  Run  tank  pressure  trace  from  shot  2448,  together  with  the  profile  obtained 
from  the  present  calculation 


present  calculation  gives  a  good  indication  of  the  mass  flow  rate  during  this  time.  However, 
the  theoretical  curve  diverges  from  the  measured  trace  at  earlier  times.  This  divergence  is 
only  slight  at  t= 0,  but  still  suggests  that  the  injection  flow  has  yet  to  be  fully  established 
at  this  time.  The  actual  mass  flux  during  the  test  time  may  thus  be  slightly  higher  than 
the  above  calculation  indicates. 


2.6  Injection  timing  experiments 


One  of  the  most  pressing  concerns  regarding  the  injection  system  was  whether  the  startup 
would  be  sufficiently  quick  that  there  would  be  adequate  injection  at  the  cone  surface  by 
the  time  of  arrival  of  the  main  test  flow.  A  number  of  experiments  involving  visualization 
of  the  injected  flow  were  thus  performed  in  order  to  gain  information  about  the  injection 
timing. 
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2.6.1  Injection  without  main  flow 

The  first  series  of  tests  involved  visualizing  the  injection  process  without  flow  in  the  test 
section  while  recording  the  unsteady  run  tank  pressure  over  a  period  of  several  seconds, 
with  the  data  acquisition  and  injection  triggered  by  the  same  signal  (generated  by  tapping 
the  accelerometer  on  the  CT).  As  such  a  run  tank  profile  had  been  obtained  during  shot 
2448  (see  figure  2.10),  test  and  shot  profiles  could  be  correlated  to  estimate  the  extent  of  the 
injected  flow  during  the  test  time  of  the  shot.  Two  sets  of  tests  were  performed:  in  the  first 
set,  the  injection  was  into  ambient  air;  in  the  second,  the  injection  was  into  the  evacuated 
dump  tank.  For  the  first  series,  as  the  cone  tip  protruded  from  the  front  of  the  test  section, 
the  dump  tank  was  simply  moved  into  a  position  so  that  the  injector  tip  was  within  the 
line-of-sight  of  the  optical  system;  for  the  second,  the  cone  was  moved  back  in  the  test 
section  so  that  the  injector  tip  was  in  line  with  the  windows.  Figure  2.11  shows  schlieren 
images  obtained  during  two  test  runs:  the  first  sequence  shows  injection  into  ambient  air, 
the  second  injection  into  the  dump  tank  vacuum.  The  initial  run  tank  pressure  and  time 
elapsed  after  triggering  in  each  image  are  indicated  in  the  caption. 


Figure  2.11:  Schlieren  visualizations  from  injection  tests:  (top)  injection  into  ambient  air 
with  an  initial  run  tank  pressure  of  149  psi;  (bottom)  injection  into  evacuated  T5  dump 
tank  with  an  initial  run  tank  pressure  of  439  psi.  In  each  case,  the  time  from  triggering  is, 
from  left  to  right:  70,  183,  and  300ms. 
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While  the  injected  gas  is  very  much  apparent  in  the  ambient  injection,  in  the  vacuum 
injection,  even  with  the  higher  run  tank  pressure,  the  flow  features  are  very  difficult  to 
discern.  Expansion  fans  are  barely  visible  at  the  injection  tips  in  the  last  image  of  the 
sequence,  but  the  startup  of  the  flow  is  impossible  to  estimate  from  these  images.  The  images 
obtained  from  the  ambient  injection  tests  were  thus  used  as  a  guide.  The  first  appearance  of 
injected  flow  in  these  tests  was  typically  visible  at  around  60  ms  after  triggering.  Visually, 
the  flow  is  fairly  well  established  within  100  ms  of  this  first  appearance,  although  a  slight 
broadening  of  the  injected  layer  continues  after  this  time. 

The  run  tank  pressure  was  also  recorded  during  these  experiments,  and  the  obtained 
traces  could  be  compared  with  that  recorded  during  shot  2448.  In  figure  2.12  the  shot  2448 
trace  is  shown  together  with  the  test  profile  at  the  initial  run  tank  pressure  closest  to  that 
of  the  shot.  Two  versions  of  the  latter  signal  are  shown:  the  first  is  the  original  signal, 
while  the  second  has  been  scaled  on  the  pressure  axis  so  that  the  initial  pressure  matches 
that  of  the  shot,  and  the  time  axis  has  been  shifted  so  that  the  initial  parts  of  the  profiles 
overlap.  The  best  match  was  found  for  a  shift  of  182  ms,  indicating  that  the  delay  between 
the  accelerometer  signal  trigger  and  the  triggering  of  the  data-acquisition  system  by  the 
reservoir  transducers  in  the  shot  was  approximately  182  ms  (with  the  test  time  beginning 
approximately  1ms  after  the  triggering  of  the  DAS).  This  is  important  information,  as  it 
enables  us  to  return  to  the  injection  images  and  estimate  to  what  extent  the  injected  flow 
had  been  established  at  the  time  of  arrival  of  the  main  flow  in  the  test  section  in  the  shot. 

In  figure  2.13,  the  run  tank  pressure  trace  from  the  test  shown  in  figure  2.12  is  again 
plotted,  this  time  with  schlieren  images  showing  the  state  of  the  injected  flow  at  several 
relevant  times.  Also  included  on  this  plot  are  the  voltage  signals  from  two  photomicrosensors 
attached  to  the  ball  valve,  showing  the  timing  of  the  valve  opening.  The  injected  flow  at  the 
estimated  beginning  of  the  test  time,  at  t=183ms,  is  shown.  The  flow  at  this  point  is  fairly 
well  established,  with  only  a  slight  broadening  of  the  injected  layer  occurring  afterwards.  If 
anything,  the  startup  under  vacuum  conditions  would  be  quicker  than  under  the  ambient 
conditions  considered  here,  due  to  the  stronger  shocks  propagating  down  the  system  in 
the  former  case.  These  tests  thus  show  that  the  injection  system  developed  here  is  indeed 
suitable  for  the  requirements  of  the  experiments. 
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Figure  2.12:  Run  tank  pressure  traces  from  shot  2448  and  from  an  ambient  injection  test. 
Both  the  original  test  signal  and  a  scaled/shifted  version  are  plotted:  the  latter  is  scaled 
on  the  pressure  axis  so  that  the  initial  pressures  match,  and  shifted  on  the  time  axis  (by 
182  ms)  so  that  the  initial  parts  of  the  profiles  overlap. 

2.6.2  Low-enthalpy  tests 

Finally,  a  series  of  low-enthalpy  experiments  were  run  in  the  facility,  in  the  configuration 
outlined  in  section  2.1  of  this  chapter.  As  the  time  from  piston  launch  to  flow  arrival  in  the 
test  section  is  different  in  low-enthalpy  experiments  to  operation  under  normal  conditions, 
these  experiments  could  not  be  used  as  a  realistic  test  for  the  injection  timing.  The  intent 
was  rather  to  visualize  injection  under  conditions  in  which  the  main  test  section  flow  was 
present.  In  figure  2.14,  visualizations  from  five  low  enthalpy  shots  with  different  run  tank 
pressures  are  shown:  run  tank  pressures  ranged  from  361  to  16psi,  and  one  shot  was 
performed  without  injection,  As  may  be  seen,  for  these  low-enthalpy  cases,  the  injected  gas 
quickly  extends  out  past  the  conical  shock.  The  stagnation  pressure  in  these  low-enthalpy 
shots  is  only  of  the  order  of  2  MPa,  whereas  in  a  normal  shot,  this  number  would  be  closer 
to  50  MPa.  Thus,  a  significantly  higher  mass  flux  would  be  present  in  the  main  test  stream 
under  normal  conditions,  which  would  cause  the  injected  flow  to  “hug”  the  cone  more  closely. 
Also,  to  further  decrease  the  exit  angle  of  the  injected  flow,  a  new  tip  has  been  designed 
and  built.  This  tip,  discussed  further  in  section  4.2,  has  four  rows  of  holes,  allowing  the 
same  injected  mass  flux  for  a  lower  run  tank,  and  thus  injector  tip,  pressure. 
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Figure  2.13:  The  run  tank  pressure  trace  from  an  ambient  injection  test,  with  voltage  traces 
from  two  microsensors  attached  to  the  ball  valve,  showing  the  timing  of  the  valve  opening. 
t=0  corresponds  to  the  trigger  signal  being  sent  to  the  injection  system.  Schlieren  images 
from  various  points  in  the  injection  sequence  are  shown:  (top  left)  t=30ms,  opening  of 
the  valve;  (top  right)  t=62ms,  startup  of  the  injection  flow  at  the  cone  surface;  (bottom 
left)  f=183ms,  beginning  of  test  time  if  main  flow  were  present;  (bottom  right)  t=300ms, 
injection  flow  at  a  later  time. 
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Figure  2.14:  Schlieren  visualisations  from  low  enthalpy  shots  with  various  run  tank  pres¬ 
sures:  (top  left)  361  psi;  (top  right)  162  psi;  (middle  left)  50psi;  (middle  right)  16  psi;  (bot¬ 
tom)  no  injection. 
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Chapter  3 

Analysis 


This  chapter  follows  the  analysis  procedure  from  the  experiment  through  to  calculating  the 
transition  Reyolds  number. 


3.1  T5  measurements 

During  the  experiment,  the  shock  tube  initial  pressure  (Pi),  the  ambient  temperature  ( Tamb ), 
and  the  run  tank  pressure  are  recorded  as  part  of  the  checklist  procedure.  After  the  experi¬ 
ment,  the  run  tank  pressure  is  noted  again.  These  quantities  are  used  in  postprocessing  the 
recorded  data.  Data  files  from  the  old  DAS  (Data  Acquisition  System)  can  be  found  on  the 
GALCIT  server  in  /home/hyper/T5/das/shot/T5/xxxx/,  where  xxxx  is  the  4  digit  shot 
number.  Data  files  from  the  new  DAS  are  transferred  to  the  GALCIT  server  and  copied  to 
this  folder. 

Two  PCB  piezo-electric  pressure  transducers  mounted  at  points  along  the  shock  tube 
and  two  more  mounted  at  the  reservoir  near  the  end  allow  measurements  of  the  shock  speed 
and  the  stagnation  pressure  to  be  made.  The  passing  shock  is  recorded  by  these  transducers 
as  a  sharp  rise  in  pressure.  From  the  time  interval  between  the  passage  of  the  shock  past 
the  first  and  second  transducer,  and  from  the  time  interval  between  the  passage  of  the  shock 
past  the  second  transducer  and  the  two  reservoir  pressure  transducers  near  the  end  wall  of 
the  tube,  the  shock  speed  is  calculated.  The  two  transducers  near  the  shock  tube  end  give 
the  reservoir  pressure.  Together  with  knowledge  of  the  pre-shock  state  of  the  gas,  this  is 
sufficient  information  to  calculate  the  reservoir  state,  as  described  below. 
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3.2  Post-processing  of  Shock  Tube  Measurements 

Using  the  pressure  transducer  data,  the  program  T5_data_proc.m  calculates  the  shock  speed 
(Us)  in  the  shock  tube  and  measures  the  stagnation  pressure  (Pq)  at  the  end  of  the  shock 
tube.  The  version  of  T5_data_proc.m  used  for  the  CO_2  injection  experiments1  can  be 
found  on  the  GALCIT  server  in  /home/hyper/T5/das/data_proc_amy.  From  this  folder, 
open  Matlab  Version  6 2  and  run  T5_data_proc. 

In  most  of  the  range  of  T5  operation,  the  incident  shock  causes  gases  such  as  nitrogen, 
oxygen,  and  carbon  dioxide  to  become  vibrationally  excited  and  partially  dissociated.  In 
almost  all  parts  of  the  range  the  vibrational  and  dissociational  relaxation  distances  are 
negligibly  small,  so  the  assumption  that  the  gas  outside  of  a  thin  shock  is  everywhere 
in  equilibrium  holds.  We  first  calculate  the  conditions  after  the  incident  shock  using  the 
initial  state  of  the  gas,  the  shock  speed,  and  the  equilibrium  properties  of  the  gas  using 
ESTC  (Equilibrium  Shock  Tube  Calculation),  a  program  written  by  M.  McIntosh  (Ph.  D. 
thesis,  Australian  National  University,  1971).  The  same  program  is  used  to  calculate  the 
conditions  after  the  reflected  shock.  The  pressure  after  the  reflected  shock  differs  from 
the  measured  reservoir  pressure,  unless  the  shock  tube  is  operated  in  perfectly  tailored- 
interface  conditions.  To  determine  the  reservoir  conditions,  ESTC  assumes  that  the  process 
transforming  the  gas  from  the  reflected-shock  state  to  the  measured  reservoir  pressure  is 
isentropic.  The  ESTC  calculation  thus  gives  the  complete  thermodynamic  state  of  the 
reservoir  gas. 

The  quantities  Tamb,  Pi,  Us,  and  Pq  are  then  input  into  ESTC.  At  the  end  of  the  T5 
shock  tube,  the  gases  are  shocked  with  both  the  incident  and  reflected  waves,  so  the  program 
calculates  the  thermodynamic  and  chemical  properties  for  twice  shocked  gases.  The  input 
file  is  fort. 8.  The  parameters  of  interest  are: 


Table  3.1:  ESTC  Parameters  of  Interest. 


parameter 

meaning 

n2 

C02 

air 

mixtures 

IRUN 

shot# 

2331 

2333,35,36;2444-45 

2433-40,46-49 

2332,34,37;2441-43 

ISC 

#  elements 

1 

2 

2 

3 

ISS 

#  species 

2 

5 

5 

8 

1The  more  general  version  can  be  found  on  the  GALCIT  server  in  /home/hyper/T5/das/data_processing 

2  Some  modifications  need  to  be  made  before  it  will  run  in  more  recent  versions  of  Matlab 


37 


Table  3.2:  Cont’d  ESTC  Parameters  of  Interest. 


parameter 

comments 

IC 

0 

#  ions,  excluding  E- 

NFIT 

1 

thermofit  cards:  1  used,  0  none  used 

KHO 

1 

harmonic  oscillator  cards:  1  used,  0  none  used 

CTMXX 

5000 

temp  [K]  where  switch  is  made  from  harmonic  oscillator  to  thermofit  data 

BZERO 

0 

constant  used  in  imperfect  gas  correction 

THETA 

0 

wedge  angle;  0  for  vertical  wall 

IGJ 

1 

thermofit  data  available  for  species  J 

IGJ 

2 

only  thermofit  (no  harmonic  oscillator)  data  for  species  J;  use  for  CO2 

From  the  output,  the  stagnation  temperature  (To),  stagnation  enthalpy  (ho)3 ,  and  stag¬ 
nation  composition  are  used  in  the  next  step. 


3.3  Nozzle  expansion 

There  are  two  programs  that  calculate  the  flow  conditions  after  nozzle  expansion.  Either 
NENZF  or  Xscriptnozflow.m,  the  nozzle  flow  expansion  script  which  is  written  in  Matlab 
and  uses  Cantera,  can  be  used. 


3.3.1  Stagnation  Conditions  to  Throat 

Starting  from  the  stagnation  conditions  calculated  in  ESTC,  an  isentrope  on  the  P  —  v 
diagram  is  followed  to  reach  the  sonic  condition  at  the  nozzle  throat. 

Given  two  thermodynamic  properties  of  the  initial  reservoir  condition,  stagnation  tem¬ 
perature  (To)  and  stagnation  pressure  (To),  and  the  molar  concentrations  of  each  species  in 
the  reservoir,  all  the  thermodynamic  properties  are  determined  in  Cantera.  The  program 
makes  note  of  the  specific  volume  (vo),  stagnation  enthalpy  (ho),  stagnation  entropy  (so), 
and  the  soundspeed  of  the  gas.  Since  we  are  considering  the  reservoir  conditions,  the  fluid 
is  taken  as  quiescent,  i.e.  u  =  0. 

'^Please  note  that  for  CO2  gas,  8.933;  MJ/kg  should  be  added  to  the  value  given  by  ESTC  to  compensate 
for  the  enthalpy  of  formation  at  STP.  The  enthalpy  in  this  output  is  measured  relative  to  molecular  oxygen, 
which  is  assigned  a  heat  of  formation  at  STP  of  zero.  For  fluid  mechanical  purposes,  however,  it  is  convenient 
to  think  of  the  enthalpy  relative  to  carbon  dioxide. 
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3.3.2  Expansion  Through  the  Nozzle 

The  following  derivation  obtains  the  differential  equations  for  nozzle  flow  solved  in  Mat- 
lab  code.  Starting  with  the  mass  (3.1),  momentum  (3.2),  energy  (3.3)  and  species  (3.4) 
equations, 


1  dp  1  dU  1  dT 
p  dx  JJ  dx  A  dx 

=  0  = 

1  dp  1  d  U  1  dT 
p  dx  U  dx  T  dx 

(3.1) 

pUdXJ 

=  -dP 

dP  TTdU 

=*■  *  =  ~pVd^ 

(3.2) 

,  u2 
h+T  = 

ht  =* 

dh  ^d  U 
dx  dx 

(3.3) 

up  = 

dx 

-WiUi 

P 

(3.4) 

where  yt  is  the  mass  fraction  of  species  i  and  u )j  is  the  production  rate  of  species  i. 
For  a  mixture  of  species,  specific  heat  and  enthalpy  are  written  as: 


Cp  —  ^  ^  ViCpi 


h  = 


d  h 
dx 


i= 1 
k 

^Vih 

1=1 

dh 


dh 

dT 


Vi 


dh  , 

Tj  - 

oyi 


dh  =  — 

dT 


Vi 


i=l 


„  dT  ^ 
CPdx  +  I> 
1=1 


d  Vi 
dx 


(3.5) 


Solve  (3.5)  for  the  temperature  gradient  and  then  use  (3.3)  to  simplify. 


dT 

dx 


J_dh  _  hj  d jjj 
Cv  dx  Cv  dx 

^  i= 1  v 

_U_dU  % 

Cv  dx  Cv  dx 

y  i=l  y 


(3.6) 


Prescribing  an  equation  of  state 


1  dP  1  dp  1  d  R  1  dT 
P  dx  p  dx  R  dx  T  dx 


P  =  pRT 


(3.7) 
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Using  the  mass  dependent  gas  constant  expression 


R~w 


Rv  =  RW 


where 


k 

2=1 


and  Xi  is  the  mole  fraction  of  species  i.  Then  the  gradient  of  the  gas  constant  becomes: 


d  R 
dx 

d  W 
dx 


d  R 
dx 


Ru  d  W 
~Wi  "dx" 

Vi  \  1  dyi 

^  Wi  dx 
2=1 

y-  Ru  d  m  _  y^  IT  d  m 

Wi  dx  Wi  dx 

i= 1  i=  1 


Substitute  this  result  into  (3.7)  and  we  have  the  following,  which  we  subsequently  solve 
for  the  temperature  gradient  and  equate  to  (3.6). 


1  d  P 

Id  p  A  W  dyi  1  dT 

p  dx  ^  Wi  dx  T  dx 

2=1 

P  dx 

dT 

dx 

=  T 

’  1  dP  1  dp  A  W  dVi 

P  dx  p dx  ^ Wi  dx 

2=1 

U  d  U 
Cp  dx 

y^  hi  d i/i 

~  ^C~p  d^ 
2=  1  y 

=  T 

’  1  dP  1  dp  A  W  dVi 

P  dx  p  dx  ^ Wi  dx 

L  '  i= 1 

Divide  this  result  by  T  and  substitute  (3.1)  and  (3.2)  to  the  right  hand  side.  Then 
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collect  like  terms  and  multiply  by  U. 


U  dU  y,  hi  d  Vi 
CPT  dx  CpT  dx 


dU 

dx 

dU 

dx 


U  PU  1 

CpT  P  U 

U 2  pU2 
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pUdU  Id  U  1  dA  W  dyi 
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UdA 
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where  p  is  defined  as  the  sonic  parameter.  Substituting  p  back  into  (3.8),  we  obtain 


dU 

dx 


UdA 
A  dx 


E 

2=1 


W 

wi 


ho 


CpT 


C-A 

dx 


which  gives  us  the  following  set  of  differential  equations  to  solve. 


dU 
dx 
dP 
dx 
dp 
dx 
d  Vi 
dx 


1 

V 


UdA 
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pU 
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where  the  thermicity,  a,  is  defined  as: 


3.3.3  NENZF 

Although  the  shock  tube  flow  remains  in  equilibrium,  the  flow  through  the  nozzle  does 
not.  While  the  gas  is  in  equilibrium  in  the  reservoir  and  through  the  throat,  and  even  up 
to  an  area  ratio  of  3  or  more,  the  chemical  reaction  rates  become  important  as  the  gas 
density  decreases  in  the  nozzle  expansion.  Nozzle  expansion  is  calculated  with  NENZF 
(Non-Equilibrium  NoZzle  Flow),  a  program  developed  by  Lordi,  Mates  and  Moselle  (NASA 
CR-472,  1966).  This  is  an  inviscid  quasi-one-dimensional  flow  code.  The  input  information 
for  this  program  consists  of  the  temperature  and  pressure  in  the  reservoir,  the  nozzle  shape, 
the  equilibrium  gas  properties,  the  chemical  reactions  that  are  active,  and  the  reaction  rate 
constants.  Vibrational  excitation  may  be  set  to  be  in  equilibrium  or  frozen. 

The  most  recent  version  of  the  NENZF  executable  is  nenzfh.  The  original  nenzf.f  com¬ 
piles  using  f77  in  the  unix  environment  only.  Modifications  would  have  to  made  to  the 
Fortran  script  in  order  to  run  it  in  the  linux  environment.  The  other  two  existing  Fortran 
scripts,  nenzfg.f  and  nenzfh. f,  compile  using  f77  in  linux.  The  executable  nenzfg  is  similar 
to  nenzfh  with  the  exception  that  To  and  Pq  are  entered  when  prompted  instead  of  entered 
in  the  input  file.  The  input  file  in  all  cases  (nenzf,  nenzfg,  and  nenzfh)  is  fort. 8,  while  the 
output  file  is  fort. 9.  The  input  file  must  be  in  the  same  folder  as  the  NENZF  executable. 
The  parameters  of  interest  are  summarized  in  the  following  table.  IRUN,  ISC,  ISS,  IC, 
NFIT,  KHO,  CTMXX,  BZERO,  and  IGJ  are  the  same  as  for  ESTC.  See  Table  3.1  above. 

The  output  quantities  are  given  in  the  same  units  as  the  output  of  ESTC.  The  enthalpy 
is  again  referenced  to  an  assigned  value  of  zero  for  molecular  oxygen,  and  the  species  con¬ 
centrations  are  given  in  moles  per  gram  of  mixture.  However  these  are  also  given  as  mole 
fractions  at  the  very  end  for  the  free-streanr  gas  composition.  The  two  numbers  called  St 
and  REYNOLDS  N  are  the  normalizing  factor  to  produce  the  Stanton  number  from  heat 
flux  value,  and  the  product  of  free-streanr  density  and  free-streanr  velocity.  These  are  merely 
convenient  numbers  for  determining  Stanton  and  Reynolds  numbers. 

The  NENZF  outputs  for  all  the  shots  are  stored  in  a  directory  in  the  T5  computer  and 
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Table  3.3:  NENZF  Parameters  of  Interest. 


parameter 

comments 

ISR 

NQS 

Number  of  reactions  -  depends  on  mechanism  model  used. 

IUPD 

0 

Starting  downstream  of  throat. 

IUPD 

1 

Starting  upstream  of  throat. 

KKUR 

1 

Third  body  reactions  used. 

KKUR 

0 

No  third  body  rxns  used. 

NAFIT 

1 

Fitted  (e.g.  contour)  nozzle  used. 

NAFIT 

0 

Standard  (e.g.  conical)  nozzle  used. 

CTAP 

xxxx 

Stagnation  temperature  in  K. 

PRESA 

CXMAX 

x.xxE+07 

Stagnation  pressure  in  Pa. 

Maximum  (non-dimensional)  distance  x. 

CXMAX 

100.0 

Length  of  contour  nozzle. 

CXMAX 

109.0 

Length  of  conical  nozzle. 

SL 

CECHII 

ALPIJ 

1.0 

Characteristic  length  in  cm,  used  to  non-dimensionalize  x. 

Check  that  there  are  ISR  number  of  entries  for  CECHII. 

The  first  ISC  number  of  rows  of  this  matrix  must  represent  the  elements. 

are  available  on  demand. 

NENZF  output  seemed  to  be  incorrect  for  the  contour  nozzle  geometry  used  to  accom¬ 
modate  the  slender  cone  in  the  test  section,  so  an  alternate  program  was  used  for  these 
calculations.  However,  NENZF  seemed  to  work  for  conical  geometry. 


3.3.4  Nozzle  Flow  Matlab/Cantera  Script 

For  the  Matlab  script,  we  input  To,  Po>  the  test  gas  stagnation  composition,  the  nozzle 
geometry  (contour,  conical  or  constant  area),  the  gas  mechanism  (thermodata  and  reaction 
rates),  and  the  length  of  the  cone  (i.e.  the  distance  at  which  to  stop  the  integration). 

Xscriptnozflow.nr  is  a  short  Matlab  script  containing  a  list  of  shot  numbers  which  calls 
nozflowfun.m  recursively  until  all  the  desired  shots  have  been  processed.  This  in  turn 
calls  the  following  non-Cantera  subroutines:  nozflowfun.m,  oneDflow.m,  areafun.m,  non- 
ideal_eq_soundspeed.m,  nonideal_soundspeed.nl,  and  isenfun.m.  These  subroutines  make 
use  of  subroutines  that  are  included  in  the  Cantera  package:  importPhase.m,  set. nr,  pres¬ 
sure. nr,  density.m,  temperature. nr,  entropy _mass. m,  enthalpy _mass. m,  moleFraction.m,  cp_mass.m, 
cv_mass.ni,  nSpecies.m,  massFractions.m,  moleFractions.m,  netProdRates.m,  molecular- 
Weights.m,  and  meanMolecular Weight. m. 

The  output  file,  which  is  in  the  form  T5_xxxx_output.plt  gives  the  initial  stagnation 


43 


conditions,  the  sonic  conditions  at  the  throat,  and  the  flow  properties  at  the  end  of  the 
nozzle. 


3.4  Taylor-MacColl 

After  exiting  the  nozzle,  the  flow  passes  through  a  conical  shock  which  sits  at  some  angle 
from  the  surface  of  the  cone.  To  find  the  flow  properties  at  the  edge  of  the  boundary  layer, 
the  Taylor-MacColl  solution  is  used  to  solve  for  flow  conditions  through  conical  shocks. 
First  the  ratio  of  specific  heats,  the  freestream  Mach  number,  the  cone  half-angle,  and  an 
initial  guess  of  the  conical  shock  angle  are  entered.  For  the  cases  in  the  present  experiments, 
an  initial  guess  of  12  or  13  degrees  was  usually  sufficient.  Once  the  program  converges  on  a 
solution,  the  program  prompts  the  user  for  the  specific  gas  constant,  the  stagnation  enthalpy 
(which  was  obtained  from  ESTC),  the  freestream  temperature,  and  the  freestream  pressure. 
Here,  freestream  refers  to  the  flow  at  the  nozzle  exit.  The  Taylor-MacColl  program  gives 
the  flow  properties  at  the  edge.  In  order  to  calculate  the  Stanton  and  Reynolds  numbers 
at  each  thermocouple,  the  only  remaining  properties  needed  are  the  boundary  and  edge 
viscosities.  The  assumption  that  the  conical  shock  is  a  weak  one  allows  for  the  use  of  the 
freestream  composition  as  the  edge  composition. 


3.5  Viscosity 

Using  the  edge  temperature  and  the  freestream  composition,  the  viscosity  can  be  calculated. 
Wilke’s  Rule  |27_  is  used  for  the  viscosity  of  the  mixed  gases.  A  Matlab  script  was  written 
by  Eric  Marineau  for  the  mixing  of  Ar ,  O2,  and  IV2  for  air.  This  was  modified  by  Ivett 
Leyva  to  include  CO2  and  its  subspecies.  The  edge  temperature  is  then  used  to  calculate 
the  reference  temperature.  The  reference  temperature  for  air,  T* .  given  by  Eckert  [8  is 

*  j  ^  t  j  * 

—  =  0.5  +  0.039 Ml  +  0.5  ,  (3.9) 

J-e  J-e 

where  Tw  is  the  wall  temperature  (taken  to  be  Tamb,  the  ambient  temperature)  and  the 
subscript  e  indicates  the  edge  conditions. 

Dorrance  |7_  showed  that  a  more  general  expression  could  be  derived  from  similarity 
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considerations  as 

=  0.5+  Ml  +  0.5^  (3.10) 

Te  2  6  Te  v  ' 

where  Pr/am  is  given  by  yfy  and  Prtur6  is  given  by  ^7. 

The  viscosity  is  recalculated  using  the  reference  temperature  so  that  the  Stanton  number 
vs  the  reference  Reynolds  number  can  be  plotted.  When  both  axes  are  on  the  log  scale,  a 
slope  of  -1/2  indicates  a  laminar  boundary  layer.  A  sharp  jump  indicates  transition. 


3.6  St  vs  Re* 

The  discussion  in  this  section  follows  that  of  Adam  [lj.  Similar  material  may  also  be  found 
in  the  thesis  of  Germain  [llj,  although  a  mistake  appears  in  this  work  in  equation  2.11b;  the 
correct  form  is  given  by  Adam  in  equation  4.14  (here  equation  3.21).  Refer  to  White  [25 
for  a  fuller  treatment  of  this  material. 

It  is  desired  to  reduce  the  experimentally  obtained  heat  transfer  profiles  to  suitable 
nondimensional  variables:  the  heat  transfer  rate  is  typically  expressed  as  the  Stanton  num¬ 
ber,  whereas  the  distance  along  the  cone  surface  is  normalized  into  the  Reynolds  number. 
The  Reynolds  number  is  defined  as 


D  PeUeX 

Re  =  - 

Pe 


(3.11) 


where  x  is  the  distance  along  the  centerline  from  the  cone  tip,  and  the  subscript  e  indicates 
that  the  fluid  properties  are  to  be  evaluated  at  the  edge  of  the  boundary  layer.  The  Stanton 
number  is  given  bythe  following  expression 


peue  [h0  -\u2e{l-r)  -  CpTw\  ’ 

where  q(x)  is  the  heat  transfer  rate  and  r  is  the  recovery  factor.  For  laminar  flows,  riam  = 
\J Pr.  and  for  turbulent  flows,  rturb  ~  Pr1^ ,  where  Pr  is  the  Prandtl  number.  Pr  is  assumed 
to  be  constant,  which  is  a  reasonable  approximation  under  the  conditions  of  interest. 

Theoretical  or  empirical  expressions  relating  the  Stanton  number  to  the  Reynolds  num¬ 
ber  have  been  obtained  for  both  laminar  and  turbulent  boundary  layers.  These  are  discussed 
separately  in  the  following. 


(3.12) 
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3.6.1  The  laminar  case 


For  self-similar  boundary  layers  at  sufficiently  high  speeds,  regardless  of  whether  the  flow 
is  laminar  or  turbulent,  the  Renolds  analogy  postulates  that  the  Stanton  number  may  be 
related  to  the  skin  friction  coefficient,  Cf(x )  =  Tw(x)/^peU2,  where  tw(x)  is  the  shear  stress 
at  the  wall,  as  |25j 

St  «  ^CfPr~2/3.  (3.13) 

Thus,  if  an  expression  relating  the  skin  friction  coefficient  to  the  Reynolds  number  can  be 
found,  this  would  give  the  desired  Re-St  relationship. 

For  an  incompressible  two  dimensional  boundary  layer  in  a  zero  pressure  gradient,  the 
classical  Blasius  flat  plate  solution  may  be  solved  numerically,  allowing  the  incompressible 
skin  friction  coefficient  to  be  obtained  as 


0.664 

CfMc{x)  =  -=.  (3.14) 

yJRex 

If  compressible  effects  are  to  be  included,  a  correction  factor  must  be  introduced.  The 
now-approximate  expression  then  becomes 


C 


f,comp 


(x)  = 


0.664 


y/Re 


(3.15) 


where  Cw  =  pwpw/ pepe  is  the  Chapman- Rubesin  parameter  evaluated  at  the  wall.  While 
this  expression  is  quite  accurate  for  an  adiabatic  wall,  for  a  hot  or  cold  wall  unfortunately 
it  is  not  as  good.  This  problem  may  be  remedied,  however,  by  employing  the  reference 
temperature  of  either  Eckert  or  Dorrance,  as  discussed  in  the  previous  section  and  given  in 
equations  3.9  and  3.10.  Then,  using  C*  =  p*p*/pepe , 

Cf}Comp(x)  ~  (3.16) 

yJRex 


A  final  consideration  is  that  the  conical  geometry  of  the  case  under  consideration  must 
be  accounted  for.  This  is  done  by  making  use  of  the  Lees-Illingworth  transformation,  which 
results  in 

C/,cone  =  '/sc f  plate-  (3-17) 
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The  desired  St-Re  relationship  for  the  cone  is  thus 

„  0.664  \/3C* 

2  Pr2/3y/R, Tx' 


(3.18) 


3.6.2  The  turbulent  case 


As  in  the  laminar  case,  we  begin  with  equation  3.13  as  the  basis  for  our  desired  St-Re 
relationship,  again  seeking  to  express  the  skin  friction  coefficient  in  terms  of  the  Reynolds 
number.  While  no  exact  expression  exists  in  the  turbulent  case,  a  good  approximation  for 
the  incompressible  turbulent  skin  friction  on  a  flat  plate  is  given  by  (see  White  [25  ): 


Cftinc  — 


0.455 


log  (0.06i?ex) 


(3.19) 


It  is  now  desired  to  introduce  a  compressibility  compression  factor,  equivalent  to  the 
Chapman-Rubesin  parameter  in  the  laminar  case.  This  can  only  be  done  if  it  is  assumed 
that  the  compressible  turbulent  boundary  layer  does  not  differ  significantly  in  a  qualitative 
manner  from  an  incompressible  one.  According  to  the  hypothesis  of  Morkovin  [17_ ,  this  holds 
if  the  density  fluctuations  are  small,  which  holds  for  boundary  layers  at  Mach  numbers  of 
up  to  around  5.  It  is  typical  to  introduce  to  correction  factors  as: 


Cf:comp  —  , .  C fine  (  R^x  S' Re  ) 

r, 


1  0.455 

i?clog2(0.06RexFRe)’ 


(3.20) 


where  Ffte  is  the  ’’stretching”  factor,  modifying  the  Reynolds  number  within  the  incom¬ 
pressible  equation  3.19,  and  Fc  is  the  correction  factor  for  the  resulting  skin  friction. 


Two  formulations  for  Fne  and  Fc  are  considered:  the  first  by  van  Driest  [24  ,  the  second 
by  White  and  Christoph  [26J.  The  former  is  the  more  widely  used.  In  both  models,  Fc  is 
given  by 


F  =  Taw/Te  -  1 
c  (arcsin  A  +  arcsin  B)2  ’ 

where  Taw  is  the  adiabatic  wall  temperature  defined  as 


(3.21) 
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The  parameters  A  and  B  are  given  by 


2a2  —  b 
Vb2  -  4 a2  ’ 


\Jb2  -  4 a2  ’ 


where  a  and  b  are  defined  as 
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(3.23) 
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(3.25) 

(3.26) 


The  model  differ  in  the  specification  of  the  stretching  factor  Fr£.  In  the  van  Driest 
model,  commonly  referred  to  as  ”  van  Driest  II” ,  this  is  given  by 


FRe 


1  He 
Fc  Hie 


whereas  in  White  and  Christoph’s  model,  it  takes  the  form 


FRe 


1  He 
\/ Fc  Hw 


(3.27) 


(3.28) 


Allowance  must  finally  be  made  for  the  fact  that  the  flow  is  conical  rather  than  over 
a  flat  plate.  For  turbulent  conical  flows,  the  transformations  used  in  the  laminar  case  are 
approximate,  and  yield  a  factor  of  approximately  1.10,  in  comparison  to  \/3  in  laminar 
conical  flow  (van  Driest  [23 J. 
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Chapter  4 

Results  and  discussion 

4.1  Results  of  present  experiments 

Over  the  last  two  years,  five  series  of  experiments  relevant  to  the  current  study  have  been 
conducted.  Tables  of  run  conditions  for  each  of  these  series,  together  with  heat  transfer 
profiles,  where  applicable,  may  be  found  in  Appendix  A. 

4.1.1  N2/C02  mixtures 

The  first  series  consisted  of  testing  with  mixtures  of  N2  and  C02  in  various  ratios  as  the 
test  gas,  with  no  injection.  These  runs  serve  as  a  baseline  to  assess  the  effect  of  adding  C02 
to  delay  boundary  layer  transition  because  C02  is  not  being  introduced  at  the  boundary 
layer  and  so  the  results  are  independent  from  any  given  injection  scheme.  The  average  po 
for  these  runs  was  51  MPa  and  the  average  mixture  ho  was  10.7MJ/kg.  The  transition 
location  for  each  shot  was  deduced  from  the  heat  transfer  profile  obtained  along  the  length 
of  the  cone.  Figure  4.1  shows  a  typical  heat  flux  profile,  and  in  this  case  the  beginning  of 
transition  from  laminar  to  turbulent  flow  in  the  boundary  layer  is  seen.  That  the  flow  is 
initially  laminar  is  shown  by  the  good  agreement  between  the  experimental  data  and  the 
theoretical  prediction  for  laminar  flow.  Then,  as  the  Re  increases,  a  sharp  increase  in  the 
heat  flux  is  seen.  The  intersection  between  the  laminar  prediction  line  and  a  least-squares 
fit  to  these  data  can  be  used  as  a  good  approximation  to  the  onset  of  transition.  From  this 
intersection,  Retr  =  peUextr/pe  (where  subscript  e  denotes  the  conditions  at  the  edge  of 
the  boundary  layer)  can  be  determined. 

The  results  of  the  experiments  in  Series  1  are  shown  in  figure  4.2,  in  which  the  transition 
Reynolds  number,  Re*,  is  plotted  against  the  N2  mole  fraction  of  the  test  gas.  The  test 
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Figure  4.1:  Typical  heat  flux  trace  for  a  transitional  flow.  Shot  2332,  Pq=51.2  MPa, 
/io=10.9MJ/kg,  N2  mole  fraction=0.95,  CO2  mole  fraction=0.05 

gas  consists  of  mixtures  of  N2  and  CO2.  Note  that  for  this  data,  a  switch  has  been  made 
from  the  Reynolds  number  based  on  the  edge  conditions  to  Re*r  =  p*Uextr/ p* ,  where  the 
density  and  viscosity  are  evaluated  at  the  reference  conditions  denoted  by  the  *  superscript 
(see  section  3.5  and  3.6).  The  use  of  such  conditions  was  shown  by  Eckert  |8j  to  be  more 
appropriate  for  discerning  trends  between  different  gases  at  different  stagnation  enthalpies. 
The  results  from  these  data  are  best  read  from  right  to  left.  For  100%  N2,  the  minimum 
value  of  Re 1r  is  obtained,  but,  as  CO2  is  added,  Re^r  increases  significantly,  more  than 
doubling  by  the  time  the  N2  fraction  is  reduced  to  60%.  The  data  pertaining  to  0%  mole 
fraction  N2  was  inconclusive:  although  transition  appeared  to  begin  toward  the  end  of 
the  cone,  there  was  not  a  clear  trend.  The  data  point  plotted  here  signifies  the  minimum 
value  that  Re^r  could  take,  corresponding  to  the  position  of  the  last  thermocouple.  Previous 
experiments  have  been  performed  with  either  pure  N2  or  pure  CO2  in  which  i?e^r  was  shown 
to  be  significantly  larger  for  CO2  flows  at  a  given  h^.  However,  this  is  the  first  time  that 
experiments  have  been  performed  in  the  T5  facility  with  mixtures  of  CO2/N2  in  the  free 
stream  for  the  purpose  of  measuring  transition  locations.  The  data  from  these  experiments 
shows  that  CO2  is  effective  in  delaying  the  transition  of  the  boundary  layer  from  laminar 
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to  turbulent  flow  when  present  as  a  component  of  the  free-stream  gas. 

The  results  shown  here  employed  the  Wilke  mixing  rule  [27J  to  compute  mixture  vis¬ 
cosities  and  the  data  for  the  viscosity  of  CO2  obtained  from  Olynick  et  al.  [18J.  This  result 
lends  optimism  to  the  effectiveness  of  the  transition  delay  technique  being  investigated  here. 


x  10® 


Figure  4.2:  Results  from  series  1  experiments:  the  effect  of  CO2  mole  fraction  on  Re^r  on  a 
freestream  otherwise  consisting  of  N2  for  a  5°  slender  cone. 


4.1.2  Air  and  air/C02  mixtures 

While  series  1  focused  on  mixtures  of  the  test  gas,  in  series  2  the  aim  was  to  find  a  pure 
air  baseline  condition  in  which  the  boundary  layer  transitioned  somewhere  near  the  middle 
of  the  cone,  so  that  when  CO2  was  added  to  the  boundary  layer  via  the  injection  system, 
the  effect  on  the  transition  location  could  be  easily  seen.  For  this  series  the  stagnation 
pressure  was  varied  from  39  to  56  MPa,  with  the  stagnation  enthalpy  falling  between  8  and 
HMJ/kg.  These  ranges  were  chosen  based  on  a  previous  condition  run  in  T5  (shot  675 
/io=10.5  MJ/kg,  po=58.5MPa),  where  transition  occurred  near  the  desired  location.  Based 
on  the  results  from  this  series  and  the  previous  results  from  shot  675,  it  was  decided  to  use 
Po~54MPa  and  /io~8  MJ/kg  as  the  target  conditions  for  the  injection  experiments. 

The  series  3  experiments  were  similar  to  those  of  series  1,  except  air/C02  mixtures  were 
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used  in  place  of  N2/CO2  mixtures.  In  addition  to  allowing  us  to  assess  the  effect  of  the 
presence  of  CO2  on  boundary  layer  transition  in  air,  this  also  enabled  an  assessment  of  any 
effects  from  the  other  elements  in  air  besides  N2.  Furthermore,  such  mixtures  represent  more 
practically  realistic  conditions,  since  if  this  technique  were  to  be  employed  in  a  hypersonic 
vehicle,  the  freestream  flow  would  be  air.  The  percentage  of  CO2  in  the  free-stream  was  0%, 
20%,  40%,  70%  and  100%  by  mass  fraction,  respectively,  for  the  various  experiments:  St-Re 
plots  from  the  first  three  of  these  are  shown  in  figure  4.3.  Comparing  the  first  two  cases, 
the  presence  of  CO2  appears  to  delay  the  onset  of  transition,  which  agrees  with  the  results 
obtained  in  Series  1  with  N2/CO2  mixtures.  However,  for  the  last  two  cases  (70%  and 
100%  CO2),  the  heat  transfer  plots  were  radically  different:  the  heat  fluxes  were  typically 
higher  by  an  order  of  magnitude  than  expected  from  the  theoretical  laminar  predictions  but 
the  profiles  had  slopes  close  to  these  predictions  (see  figure  4.3).  Therefore,  it  is  not  clear 
whether  these  two  cases  are  transitional  or  turbulent,  or,  in  fact,  whether  this  result  is  due 
to  the  thermocouples  failing  or  malfunctioning  due  to  excessive  wear.  Comparing  the  100% 
CO2  profiles  from  series  1  and  3  suggests  the  latter  possibility. 

4.1.3  Injection  experiments 

In  series  4,  for  the  first  time,  CO2  was  introduced  to  the  cone  boundary  layer  through  an 
injector  tip.  The  baseline  air  condition  determined  from  series  2  was  run  and  CO2  was 
injected  at  different  mass  flow  rates.  The  injector  employed  (injector  tip  #1)  had  72  holes 
with  diameters  of  0.0258”,  and  is  shown  in  figure  4.4.  Four  runs  were  completed,  all  with 
/io~6-7MJ/kg  and  po~46  —  56  MPa.  The  reason  these  numbers  differ  slightly  from  the 
target  condition  is  that  the  primary  diaphragms  used  previously  were  no  longer  available, 
so  this  condition  could  not  be  exactly  reproduced. 

One  of  the  obtained  heat  transfer  profiles  is  shown  in  figure  4.5;  the  remainder  may 
be  found  in  Appendix  A.  As  may  be  seen,  the  St-Re  profiles  for  these  experiments  are 
all  similar  to  one  another,  and  also  resemble  the  profiles  from  series  3  in  which  the  CO2 
fraction  was  60%  or  above.  The  experimental  Stanton  numbers  are  significantly  higher  than 
the  corresponding  theoretical  values,  and  the  profiles  typically  show  a  slightly  negative  slope 
(according  to  the  least-squares  fits  shown),  although  the  large  scatter  in  the  data  makes 
meaningful  conclusions  difficult  to  draw.  The  similarity  of  these  profiles  to  the  high-C02- 
fr action  experiments  of  series  3  may  indicate  one  of  two  things:  either  that,  despite  the 
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Figure  4.3:  Heat  transfer  profiles  with  different  mixtures  of  C02/air  as  the  test  gas:  (top 
left)  0%  mole  fraction  CO2;  (top  right)  20%  mole  fraction  CO2;  (bottom)  70%  mole  fraction 

C02. 


fact  that  the  main  test  gas  is  air  and  CO2  is  being  injected  only  over  a  small  area  close 
to  the  tip  of  the  cone,  the  injected  CO2  is  dictating  the  behavior  of  the  heat  transfer  over 
the  surface  of  the  cone,  or,  that  the  thermocouples  were  continuing  to  malfunction.  If  the 
former  is  true,  and  these  are  valid  data,  it  may  be  possible  to  discern  a  further  trend  in  the 
data,  namely  a  trend  towards  an  increasingly  positive  slope  in  the  St-Re  profile  from  shot 
2446  to  2449.  This  corresponds  to  both  decreasing  stagnation  pressure  and  increasing  run 
tank  pressure.  Note,  for  example,  that  the  least-squares  fit  gives  a  positive  value  for  the 
profile  slope  for  shot  2449,  which  corresponds  to  both  the  lowest  po  and  the  highest  run  tank 
pressure  among  the  shots.  This  may  indicate  a  trend  towards  transition  at  such  conditions, 
though  this  would  be  counterintuitive  since  this  shot  had  both  the  lowest  overall  Re  and 
the  largest  flux  of  injected  CO2.  Another  possibility  is  that  at  higher  injection  rates,  the 
injected  gas  is  disrupting  the  boundary  layer,  causing  an  earlier  transition  to  turbulence. 
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Figure  4.4:  The  configuration  for  the  injection  experiments  of  Series  4:  (left)  photograph  of 
injector  tip  #1  installed  in  the  cone  model;  (right)  schematic  of  the  flow  path  of  the  injector 
tip. 


Again,  however,  the  scatter  in  the  data  and  the  possibility  of  malfunctioning  thermocouples 
prevent  any  firm  conclusions  being  drawn. 


Figure  4.5:  Normalized  heat  flux  vs  Re  for  Shot  2446,  in  which  CO2  was  injected  with 
injector  tip  #1. 


With  the  possibility  of  faulty  instrumentation  in  mind,  at  this  point  it  was  decided 
to  switch  to  visualization  of  the  CO2  injection  with  the  existing  model  while  a  new  cone 
model  was  being  designed  and  manufactured.  The  new  model  was  to  have  more  than  four 
times  the  number  of  thermocouples  of  the  existing  model,  allowing  for  much  more  detailed 


measurements . 
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4.1.4  Injection  visualization 

In  the  last  of  these  series  of  experiments,  series  5,  it  was  decided  to  focus  on  visualizing 
the  injection  flow  and  validating  that  the  flow  of  CO2  would  be  established  before  the 
main  test  gas  arrives  at  the  cone.  Uncertainty  about  the  synchronization  of  the  injection 
with  the  tunnel  flow  prompted  us  to  investigate  the  timing  by  using  a  Phantom  high¬ 
speed  movie  camera  to  determine  the  delay  from  the  time  of  sending  the  trigger  signal 
to  the  injector  valve  to  the  establishment  of  injected  CO2  at  the  surface  of  the  cone.  By 
comparing  run  tank  pressure  traces  from  these  tests  and  from  shot  2448,  it  was  found  that 
this  delay  needed  to  be  less  than  ~180ms.  Further  details  of  these  experiments  are  given 
in  section  2.6.  Figure  4.6  shows  the  results  from  one  test  in  which  CO2  was  injected  into 
ambient  air.  These  experiments  demonstrated  that  the  onset  of  CO2  injection  at  the  surface 
of  the  cone  occurred  approximately  70  ms  after  the  trigger  signal.  The  injection  timing  is 
thus  compatible  with  the  operation  of  the  tunnel  and  this  result  confirms  that  the  injection 
mechanism  (valving  time,  trigger  signal,  etc.)  that  has  been  designed  is  suitable  for  this 
experiment. 


Figure  4.6:  Flow  visualization  of  CO2  injection  from  injector  tip  1  into  ambient  air  with  an 
initial  run  tank  pressure  of  398  psi:  (left)  flow  at  73  ms  after  trigger;  (right)  flow  at  78  ms 
after  trigger. 
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4.2  Future  Work 

With  experience  gained  from  the  testing  performed  thus  far,  a  new  series  of  experiments 
is  planned  with  important  improvements  to  both  the  instrumentation  and  the  injection 
method.  The  first  improvement  is  the  design  and  manufacture  of  a  new  cone,  carrying 
significantly  more  instrumentation  than  the  cone  model  described  in  this  report.  This 
older  cone  was  used  for  the  present  experiments  because,  having  been  already  built  and 
successfully  used  in  a  past  project,  it  represented  no  added  cost  to  our  limited  budget. 
As  discussed,  this  cone  performed  quite  well  in  series  1  and  2  but  the  thermocouples  may 
have  started  to  show  signs  of  aging  in  series  3  and  4.  The  new  cone  has,  in  fact,  already 
been  designed  and  built,  and  is  currently  being  instrumented.  Photographs  of  the  main 
body  are  shown  in  figure  4.7,  In  addition  to  having  80  thermocouples  in  comparison  to  the 
21  of  the  old  cone,  the  new  cone  also  has  the  thermocouples  distributed  uniformly  in  the 
circumferential  direction  as  well  as  the  longitudinal  axis,  whereas  the  old  cone  had  usable 
thermocouples  only  on  one  half. 


Figure  4.7:  New  cone  with  80  thermocouple  locations  (as  yet  uninstalled):  (left)  overall 
view;  (right)  blow-up  view  detailing  several  of  the  thermocouple  orifices. 

The  second  improvement  for  the  future  series  is  a  new  injector  tip,  which  has  been 
designed  and  manufactured.  This  tip  is  longer  than  those  previously  used,  allowing  the 
accommodation  of  a  larger  number  of  holes,  and  thus  an  increased  mass  flow  rate  for  the 
same  run  tank  pressure.  A  schematic  of  the  new  design  is  shown  in  figure  4.8, 
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Figure  4.8:  Schematic  of  the  new  injector  tip  (Injector  5)  to  be  tested  in  the  next  set  of 
experiments. 

In  the  next  series  of  experiments  it  is  intended  to  start  with  the  previously-determined 
baseline  condition  of  Fb=54MPa,  /io=8MJ/kg  with  100%  air  in  the  test  gas  and  CO2  only 
being  introduced  via  the  injection  tips.  For  this  baseline  condition,  for  which  transition 
occurs  on  the  downstream  half  of  the  cone  when  no  injection  is  present,  several  CO2  mass 
flow  rates  (high,  medium,  low)  will  be  tested  with  at  least  one  old  injector  tip  (fitted  with 
an  adapter  section  for  the  new  cone)  and  the  new  injector  tip.  It  is  expected  that  a  few 
calibration  shots  with  no  CO2  will  be  performed  to  check  the  thermocouples  in  the  new 
cone,  and  to  make  sure  previous  conditions  can  be  repeated.  We  expect  to  demonstrate 
delay  of  the  boundary  layer  transition  through  varying  of  the  CO2  mass  flow  rate  introduced 
to  the  baseline  condition. 

4.3  Conclusions 

In  this  project,  a  new  technique  to  delay  the  transition  in  hypervelocity  air  flows  by  injecting 
CO2  into  the  boundary  layer  of  interest  has  been  proposed  and  investigated.  The  motivation 
for  this  idea  comes  from  experimental  and  computational  data  that  shows  that  for  the  same 
stagnation  enthalpy,  the  transition  Reynolds  number,  Rejjj.,  is  larger  for  CO2  flows  than  for 
either  air  or  N2  flows.  The  explanation  for  this  phenomenon  lies  in  the  fact  that  when  CO2 
is  in  vibrational  and  chemical  nonequilibrium,  these  relaxation  processes  absorb  energy  from 
the  acoustic  disturbances  whose  growth  in  the  boundary  layer  is  responsible  for  transition 
in  hypervelocity  flows  (second  or  Mack  mode). 
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Five  experimental  series  have  been  conducted  in  this  project  to  test  this  technique  in 
high  enthalpy  flows  over  a  5°  slender  cone.  This  geometry  was  chosen  due  to  the  wealth 
of  previous  information  available  with  which  obtained  Re Jj,  data  can  be  compared.  A  side 
benefit  of  using  this  geometry  was  that  an  instrumented  cone  already  existed  and  therefore 
allowed  us  to  start  this  project  with  a  limited  budget.  Initial  experiments  used  mixtures 
of  CO2/N2  as  the  test  gas,  rather  than  introducing  CO2  directly  to  the  boundary  layer  of 
the  cone.  The  obtained  results  demonstrate  that  the  addition  of  CO2  delays  the  onset  of 
transition.  For  example,  in  the  case  of  60%N2/40%C02  by  mole  fraction,  the  value  of  Re^r 
was  more  than  double  in  comparison  to  the  case  of  100%  N2.  A  similar  effect  was  noted 
in  several  experiments  using  mixtures  of  air  and  CO2  as  the  test  gas.  Several  runs  were 
then  performed  in  100%  air  to  determine  a  baseline  condition  for  injection  experiments:  the 
requirement  for  this  condition  was  that  transition  occurred  near  the  mid-point  of  the  cone, 
so  that  any  adjustment  of  transition  location  due  to  injection  would  be  easily  observed. 
A  small  number  of  experiments  were  subsequently  carried  out  with  CO2  injection  into  an 
air  freestream,  but  questionable  thermocouple  readings  meant  that  no  firm  conclusions 
could  be  drawn  from  these.  Finally,  tests  were  carried  out  in  which  the  injected  flow  was 
visualized,  providing  confirmation  that  the  triggering  of  the  injection  system  is  appropriately 
timed  to  provide  adequate  injected  CO2  during  the  test  time.  In  fact,  the  CO2  flow  starts 
approximately  100  ms  before  the  main  flow  arrives,  providing  enough  time  for  the  injection 
flow  to  be  fully  established. 

A  new  cone  model  has  been  constructed  with  80  thermocouples  uniformly  distributed 
in  the  circumferential  and  axial  directions.  This  compares  with  the  older  model  that  was 
restricted  to  21  thermocouples  located  only  on  one  side  of  the  cone.  A  new  injector  tip  has 
also  been  designed  and  built,  having  twice  as  many  holes  as,  and  approximately  double  the 
flow  area  of  the  previously  used  tip.  This  will  allow  lower  injection  pressures  to  be  used 
for  the  same  mass  flow  rate.  The  next  test  campaign  is  scheduled  to  be  completed  by  the 
end  of  this  year.  The  first  focus  of  this  campaign  will  be  to  quantify  the  effect  of  different 
injection  conditions,  for  a  given  free-stream  condition,  on  transition  delay. 
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Appendix  A 

Run  Conditions  and  Data  Plots 


This  appendix  contains  the  relevant  run  condition  data  from  each  of  the  experimental  series, 
together  with  plots  of  the  Stanton  number  (normalized  heat  flux)  versus  the  Reynolds  num¬ 
ber  for  each  relevant  shot.  In  these  plots,  theoretical  curves  corresponding  to  laminar  and 
turbulent  heat  transfer  profiles  (both  the  Van  Driest  II  and  White  &  Christoph  turbulent 
rnoels  are  shown)  are  plotted  together  with  the  experimental  data,  and  a  best-fit  transition 
line  is  also  indicated  where  relevant. 


A.l  Series  1  Experiments  -  Seeded  N2  Flows 

In  the  Series  1  experiments,  the  N2  test  gas  was  seeded  with  varying  fractions  of  CO2. 


Table  A.l:  Series  1  run  conditions.  Baseline  and  seeded  flow. 


Shot  No. 

Po 

[MPa] 

To 

[K] 

h0 

[■ MJ/kg ] 

m 

[mol /mol] 

[co2] 

[mol  /  mol] 

2331 

52.2 

7375 

10.75 

1.0000 

0.0000 

2332 

51.2 

7095 

10.90 

0.9500 

0.0500 

2333 

50.0 

4434 

9.43 

0.0000 

1.0000 

2334 

50.1 

5380 

10.89 

0.6111 

0.3889 

2335 

37.6 

3459 

5.58 

0.0000 

1.0000 

2336 

37.2 

3606 

6.12 

0.0000 

1.0000 

2337 

51.2 

7379 

11.37 

0.9676 

0.0324 

Stanton  number  (St)  [-]  Stanton  number  (St)  [-] 
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Plot  of  St  vs  Rex  for  T5-2331 ;  P0  =  52.2  MPa,  hQ  =  10.8  MJ/kg,  TQ  =  7375  K 
10'° 


io-3 


io~4 

io6  io7 

Reynolds  number  (Re  )  [-] 


Reynolds  number  (Re  )  [-] 


Reynolds  number  based  on  reference  conditions  (Re  )  [-] 


Figure  A.l:  Heat  flux  profiles  from  series  1  experiments. 
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A. 2  Series  2  Experiments  -  Calibration  in  Air 

In  Series  2,  the  test  gas  was  air  in  all  cases.  A  baseline  condition  was  sought  in  which 
transition  was  observed  to  occur  at  approximately  the  midpoint  of  the  cone. 

Table  A. 2:  Series  2  run  conditions. 


Shot  No. 

Po 

[MPa] 

To 

[K] 

ho 

[M  J /  kg] 

2433 

39.1 

5305 

8.34 

2434 

42.5 

5538 

8.87 

2435 

48.9 

5843 

9.53 

2436 

49.5 

6069 

10.07 

2437 

46.8 

6567 

11.33 

2438 

48.9 

6069 

10.08 

2439 

55.5 

5715 

9.15 

Stanton  number  (St)  [-]  Stanton  number  (St)  [-]  Stanton  number  (St)  [-] 
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Reynolds  number  (Re  )  [-] 


Reynolds  number  (Re  )  [-] 


Plot  of  St  vs  Rex  for  T5-2438;  P0  =  48.9  MPa,  hQ  =  10.08  MJ/kg,  TQ  =  6069  K 


O  experimental 

- laminar 

Van  Driest  II 

:-r~  -  --  ------- 

White  &  Christoph 

io'4L 

10 


Reynolds  number  (Re  )  [-] 


Figure  A. 2:  Heat  flux  profiles  from  series  2  experiments. 
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A. 3  Series  3  Experiments  -  Seeded  Air  Flows 


In  series  3,  the  air  test  gas  was  seeded  with  varying  fractions  of  CO2.  Shot  2440  was  used 
as  a  baseline  and  CO2  was  gradually  introduced  until  the  test  gas  was  100%  CO2. 


Table  A. 3:  Series  3  run  conditions. 


shot  no. 

Po 

[MPa] 

T0 

[K] 

ho 

[MJ/kg] 

[air] 

[9/9] 

[co2] 

[9/9] 

2440 

53.8 

5216 

7.99 

1.0000 

0.0000 

2441 

52.6 

4689 

7.33 

0.8000 

0.2000 

2442 

53.3 

4250 

7.18 

0.6000 

0.4000 

2443 

53.9 

3702 

4.38 

0.3000 

0.7000 

2444 

53.0 

3616 

5.99 

0.0000 

1.0000 

2445 

53.3 

3782 

6.61 

0.0000 

1.0000 

Stanton  number  (St)  [-]  Stanton  number  (St)  [-]  Stanton  number  (St)  [-] 
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Reynolds  number  (Re  )  [-] 


Plot  of  St  vs  Rex  for  T5-2443;  PQ  =  53.9  MPa,  hQ  =  4.38  MJ/kg,  TQ  =  3702  K 

10”2r - - - * - * - * - * - * - * - ^ ■ 


Van  Driest  II 


White  &  Christoph 


Reynolds  number  (Re  )  [-] 


Figure  A. 3:  Heat  flux  profiles  from  series  3  experiments. 
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A. 4  Series  4  Experiments  -  High  Enthalpy  Injection 

In  Series  4,  CO2  was  injected  into  the  boundary  layer  of  the  cone  in  a  test  flow  consisting 
of  100%  air.  The  runtank  pressure  reading,  Prt,  indicated  below  is  the  fill  pressure  as  read 
from  the  digital  gauge,  and  the  (approximate)  mass  flow  rate  is  calculated  from  this.  This 
gauge  was  found  to  read  slightly  lower  than  the  Kulite  pressure  transducer,  and  the  initial 
value  thus  may  be  a  good  approximation  to  the  actual  tank  pressure  during  the  test  time. 


Table  A. 4:  Series  4  run  conditions.  CO2  injection  in  high  enthalpy  air  flow. 


shot  no. 


2446 

2447 

2448 

2449 


Po 

[MPa] 


To 

[K\ 


m 


ho  Prt 

[MJ/kg]  [MPa]  [kg/s] 


55.8 

50.5 

48.4 

45.7 


4742 

4357 

4358 
4286 


6.85 

6.06 

6.07 

5.94 


1.35  0.057 

1.58  0.067 

2.43  0.103 


Reynolds  number  (Re  )  [-] 


Reynolds  number  (Re  )  [-] 


Figure  A. 4:  Heat  flux  profiles  from  series  4  experiments. 
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Appendix  B 

Mat  lab  Scripts 


This  appendix  contains  the  Matlab  scripts  and  function  files  used  to  calculate  the  flow  from 
the  nozzle  through  to  the  conditions  at  the  edge  of  the  cone  boundary  layer.  Comments 
have  been  added  throughout  where  appropriate. 

B.0.1  Xscriptnozflow.m 

7. 

7.  SCRIPT:  Xscriptnozflow.m 

7.  This  script  runs  nozflowfun.m  recursively  for  all  the  shot  numbers  indicated 
clear  all; 

index  =  [2446  2447  2448  2449] ; 

[a  b]  =  size(index); 

for  i  =  l:b 

input_file  =  [‘Xshot’,  num2str (index (i) ) ,  ‘.dat’]; 
nozf 1 owf un( input _f ile) ; 

end 


B.0.2  nozflowfun.m 

7. 

7.  SCRIPT:  nozflowfun.m 

/  This  script  creates  arrays  for  P-v  Curves  to  the  sonic  point,  then  solves 
7.  for  quasi  1-D  flow  down  a  nozzle  with  area  change 

7. 

7.  Calls  non-Cantera  functions  -  nonideal_soundspeed.m,  isenfun.m, 

7.  nonideal_eq_soundspeed.m,  oneDflow.m 

function  nozf lowfun(inp) 

7.  read  in  #  of  species 

AA  =  textread(inp, ‘7.d’ ,  1 ,  ‘ commentstyle ’, ‘matlab’ )  ; 
n  =  AA(1) ; 

N  =  n  +  4; 

7.  read  in  species  concentrations 

BB  =  textread(inp, ‘7»s’ ,N, ‘commentstyle’ , ‘matlab’)  ; 


68 


mech  =  BB{N-2};  "/,  N.B.  curly  braces  to  read  inside  cell 

fname  =  BB{N-1};  "/  (to  read  them  in  as  string) 

nozzle  =  BB{N};  "/,  parantheses  for  entire  cell 

7.  read  in  initial  temperature  and  pressure 

CC  =  textread(inp, ‘ ‘headerlines ’ ,N> ‘ commentstyle ‘matlab’ ) ; 

TO  =  CC(1) ; 

PO  =  CC(2) ; 
outfile  =  CC(3); 
plotfns  =  CC(4); 
xfinal  =  CC(5); 

7.  reads  the  species  concentrations  into  a  string 

for  i  =  l:n 

Q{i>  =  BB{i+l>; 

end 

q  =  Q{1>; 

if  n  >=  2 

for  i  =  2:n 

q  =  [q, ‘  ’ ,Q{i>] ; 

end 

end 

7.  this  gets  rid  of  the  concentrations  and  leaves  a  string 
7.  with  just  the  species 
for  i  =  l:n 

S  =  sscanf  (Q{i} , ‘‘/,c  ’  )  ; 
j  =  i; 

while  (S(j)~=‘ : 3 ) && ( j  <  8) 

SS(j)  =  S ( j )  ; 

j=j+i; 

end 

spec{i}  =  SS; 
clear  SS  S 

end 


species  =  spec .  ’  ; 

7o7.7o7.7.7o7.7.7.7.7o7X/.7.7o7.7.7o7.7.7.7.7o7.7.7.7.7o7.7.7o7.7o7.7.7o7X/.7.7o7.7.7.7.7o7.7.7o7.7.7.7.7o7.7.7.7.7.7.7.7o7.7o7.7.7o7.7.7.7.7o7.7.7.7. 

7.  ISENTROPE  CURVE  ON  P-V  DIAGRAM 

7.y.y.7.y.7.y.y.y.y.7.y.y.7.y.7.7.y.7.y.y.7.y.7.y.7.7.y.7.y.y.7.y.y.7.y.7.y.y.y.y.7.7.y.7.y.7.7//.7.y.y.7.77.y.y.7.y.y.y.y.7.y.y.7.7.7.y.7.7.y.7.7.y.y.y. 


I  GET  FIRST  POINT  ON  ISENTROPE 


gas  =  importPhase (mech) ; 

set(gas, ‘T’ ,T0, ‘P’ ,P0, ‘X’ ,q) ;  "/,  sets  the  state  of  the  gas 

rO  =  density  (gas)  ;  "/,  initial  density 

st  =  entropy_mass  (gas)  ;  "/,  total  entropy 

ht  =  enthalpy _mass  (gas)  ;  "/,  total  enthalpy 


i  =  1; 
va  =  1/rO; 


"/,  specific  volume 
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P(i)  =  pressure (gas) ; 
v(i)  =  va; 

Tisen(i)  =  temperature (gas) ; 

Sisen(i)  =  entropy_mass (gas) ; 

c(i)  =  nonideal_soundspeed(gas ,mech) ; 

u(i)  =  0; 

while(u(i)  <  c(i)) 
i  =  i  +  1; 
va  =  va  +  0.001; 

"/,  Always  starts  at  new  volume  and  previous  temperature  to  find  next  state 
array  =  [va;st]; 

options  =  optimset ( 'Display ’ , ‘ of f ’) ; 

[Ta,fval]  =  f solve (Sisenfun, TO .options , gas , array) ; 

set (gas , ‘T’ ,Ta, ‘Density’ , 1/va) ; 

P(i)  =  pressure (gas) ; 
v(i)  =  va; 

Tisen(i)  =  temperature (gas) ; 

Sisen(i)  =  entropy_mass (gas) ; 

c(i)  =  nonideal_eq_soundspeed(gas ,mech) ; 

u(i)  =  sqrt(2*(ht  -  enthalpy_mass (gas) ) ) ; 

end 


7.0/.7o0/.7o7.0/.7o0/o0/.7.0/.7o0/o7.7.7o7.%0/.7o0/o7.7.7.7o0/o%%7.7o7.0/.7o0/.7o7.0/.7o0/o7.7.7o7o0/o7.7.0/o0/o0/.0/.“/o0/o7o7.7.7o0/o%%7.7o7.0/.7o7.7o7.0/.7o0/o0/.7.7o7o0/o0/. 

I  MAKE  P-v  AND  T-S  DIAGRAMS 

%7.7.7.7o7.7.7.7.7.7.7.7.7.7.7.7o7.7.7.7o7o7.7o7.7o7o7.%7.7o7.7.%7.7o7.7.7.7o7.7.7o7.7.7.7.7o7.7.7.7o7o7.7o7.7o7o7.%7.7o7.7.7.7.7o7.7.7.7.7o7.7o7o7.7. 


if  (plotfns  ==  0)  "/,  0  for  no  plots,  anything  else  to  create  plots 

status  =  ‘DONE  ISENTR0PE’ 

else 

plot(v,P) 
figure ; 

plot (Tisen, Sisen) 


end 


7.7. %7.7o7.7.7.7.7.7.7.7.7.7.7.7o7.7.7.7o7o7.7o7.7o7o7.7o7.7o7.7.7.7.7o7.7.7.7.7.7.7o7o7.7.7.7o7.7.7.7o7o7.7o7.7o7o7.7o7.7o7.7.7.7.7o7.7.7.7.7o7.7o7o7.7. 

"/.  CREATE  OUTPUT  TEXT  FILE 

7.7. %7.7o7.7.7o7.7.7.7.7o7o7«7.7o7.7.7.7o7o7.7.7.7o7o7.%7.7o7.7.%7.7o7.7.7o7.7«7.7o7o7o7o7.7o7o7o7.7o7o7o7o7.7o7o7.%7.7o7.7.%7.7o7.7.7o7.7.7.7o7o7o7o 


if(outfile  ==  0) 

disp(‘No  output  files  created  from  nozf lowfun .m’ ) 

else 

fn  =  [fname,  ‘ _output .pit ’ ] ; 
d  =  date ; 

fid  =  fopen(fn,  ‘w’); 
fprintf (fid,  ...);... 

y  =  [v ;  P ;  c ;  u]  ; 

fprintf  (fid,  ‘"/0.4f  \t  '/,  .Of  \t  '/, .  2f  \t  '/, .  2f  \n\n\n’ ,  y)  ; 
f close(f id) ; 

end 


"/.  CALCULATE  SONIC  CONDITIONS 
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7.%%yoyo7.%%yo%7.%%7o%7.%%%%7o0/o0/.7.7.“/o7o0/.%0/o0/o7.%7.0/o7o7.%7o0/o7.7.7o0/o7o0/.7.0/o7o7o7.7o0/o0/o7o7.“/o7o0/.%0/o7o7.0/.7.0/o7o7.%7o0/o7.7.7o0/o°/o0/. 


mFs  =  moleFractions  (gas)  ;  7  mole  fractions  at  sonic  condition 

Ts  =  temperature  (gas)  ;  "/,  temperature  at  sonic  condition 

Ps  =  pressure  (gas)  ;  7  pressure  at  sonic  condition  [Pa] 

rs  =  density  (gas)  ;  7  density  at  sonic  condition 

cs  =  nonideal_soundspeed(gas ,  mech) ; 

Cp  =  cp_mass (gas) ; 

Cv  =  cv_mass (gas) ; 
gam  =  Cp/Cv; 

xl  =  0.001;  7  put  xl  slightly  downstream  of  throat 

cl  =  cs  +  1;  */.  make  the  velocity  slightly  supersonic 

7  another  way  of  expressing  mole  fractions  at  sonic  condition 
qs  =  0; 
for  i  =  l:n 

qson-(i}  =  sprintf  ( ‘’/,s"/„s'/,f  ’  ,  spec{i} ,  ‘  ’  ,mFs  (i)  )  ; 
qs  =  [qs,  ‘  ’  ,qson-[i]-]  ; 

end 


7.7. %7.7o7.7.7o7o7.7.7.7o7.7.7.7o7o7.7.7o7o7o7.7.7o7o7.%7.7o7.7.%7.7o7.7.7.7.7.7.7o7o7o7.7.7.7.7.7.7o7o7o7.7.7o7o7.%7.7o7.7.7.7.7o7.7.7.7o7.7.7o7o7o7. 

7.  RUN  l-D  FLOW  CALCULATION  w/AREA  CHANGE 

7.7. %7.7o7.7.7o7o7.7.7.7o7.7.7.7o7o7.7.7o7o7o7.7.7o7o7.%7o7o7.7.%7.7o7.7.7o7o7.7.7o7.7.7.7.7o7o7.7.7o7o7.7o7.7o7o7.7o7o7o7.7.7.7.7o7.7.7o7o7.7.7o7o7o7o 


nsp  =  nSpecies (gas) ; 

yO  =  [1  7  initial  pressure/Ps 

rs  "/,  initial  density 

cl  "/,  initial  velocity  -  slightly  supersonic 

massFractions(gas)] ; 

xel  =  [xl  xfinal]  ; 

options  =  odeset ( ‘RelTol 1 . e-5 , ‘AbsTol’ , 1 . e-8 , ‘Stats ‘ of f ’) ; 
tO  =  cputime; 

warning  off  MATLAB : odel5s : IntegrationTolNotMet ; 

7,  out.x  =  position,  out.y  =  temperature  and  species  profiles 
out  =  odel5s (OoneDf low, xel ,y0 , options , gas ,Ps , nozzle) ; 


7.7.7.7.7o7.7.7.7.7.7.7.7.7o7.7.7o7.7o7.7o7o7.7o7.7o7o7o%7.7o7.7.7.7.7o7.7.7.7.7.7.7o7.7.7.7.7o7.7o7.7o7o7.7o7.7o7o7.%7.7o7.7.7.7.7o7.7.7.7.7.7.%7.7.7. 


7.  EXTRACT  OUTPUT  INFORMATION 


%7.7.7.7o7.7.7.7.7.7.7.7.7.7o7.7o7.7o7.7o7o7.7o7.7o7o7.7.7.7o7.7.7.7.7o7.7.7.7.7.7.7o7.7.7o7.7o7.7o7.7o7o7.7o7.7o7o7o%7.7o7.7.7.7.7o7.7.7.7.7.7.7o7.7.7. 


[a,b]  =  size(out.y); 

7o  a  =  3  +  species;  b  =  number  of  iterations 

7.  Species  profile  for  current  iteration 
yy  =  zeros (nsp, 1) ; 


7.  GET  POSITION,  PRESSURE,  DENSITY,  AND  TEMPERATURE  ARRAYS 
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throat  =  areafun(0 .nozzle) ; 

for  i  =  l:b 

for  m  =  4:a 

yy(m-3,l)  =  abs (out .y (m. i) ) ; 

end 

rho(i)  =  out.y(2,i); 

set (gas , ‘Rho ’ ,rho(i) , ‘ Y’ .transpose (yy) ) ; 
wt(i)  =  meanMolecularWeight (gas) ; 

R(i)  =  gasconstant/wt (i) ; 

T(i)  =  out .y (1 , i) *Ps/ (out .y (2 , i) *R(i) ) ; 

pos(i)  =  out.x(i); 

area  =  areafun(pos (i) .nozzle) ; 

AR(i)  =  area/throat; 

Pr(i)  =  out .  y  (1 ,  i)  *Ps ;  7  pressure  [Pa] 

set (gas , ‘T’ ,T(i) , ‘Rho ’ , rho (i) , ‘ Y’ .transpose (yy) ) ; 

"/.  Vectors 

wdot  =  netProdRates (gas) ; 
mw  =  molecularWeights (gas) ; 

"/.  Scalars 
cp  =  cp_mass (gas) ; 
cv  =  cv_mass (gas) ; 
g(i)  =  cp/cv; 

c(i)  =  nonideal_soundspeed(gas ,mech) ; 

U(i)  =  out ,y(3,i) ; 

M(i)  =  U(i)/c(i);  7  Mach  Number 

mF(i,l:nsp)  =  moleFractions (gas) ; 

end 

Pexit  =  Pr(b); 
rexit  =  rho(b); 
xexit  =  pos (b) ; 

Texit  =  T  (b) ; 

Rexit  =  R(b) ; 
wtexit  =  wt(b); 

Uexit  =  U(b) ; 

Mexit  =  M(b) ; 
cexit  =  c (b)  ; 
gexit  =  g(b) ; 
mFexit  =  mF(b,l:nsp); 


7777777777777777777777777777777777777777777777777777777777777777777777777777.7 

7  APPEND  TO  OUTPUT  TEXT  FILE 

77777777777777777777777777777777777777777777777777777777777777777777777777777 


if(outfile  ==  0) 

disp(‘No  output  files  created  (script_nozzle) ’ ) 

else 

fid  =  fopen(fn,  ‘a’); 
fprintf (fid, 
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y  =  [pos;  rho;  Pr;  T;  U] ; 

fprintf(fid,  ‘°/01.6e  \t  "/,1.6e  \t  7,1. 6e  \t  7,1. 6e  \t  7,1. 6e  \n’,  y)  ; 

fprintf (fid, 
f close(f id) ; 

end 


%0/.%7.7o7.0/.7o0/o7.7.%0/.7o0/.7.7o7o0/o7.7o0/o0/o7.7.7o0/oy.7o7.7o7.0/.%7.7o7.0/.7o0/o7.7.0/o7.7.0/.7.7o7.0/o7.7o0/o0/.7.7.7o7oy.7o0/o7o7.0/.%7.7o7.%7o0/o0/.7.7o7o°/o0/. 

7.  MAKE  PLOTS  OF  P,  T,  U,  AND  CONCENTRATIONS  VS  AREA  RATIO 

7.7.%7.7o7.7.7.7.7.7.7.7.7.7.7.7o7.7o7.7o7o7.7.7.7o7o7.7o7.7o7.7.%7.7o7.7.7.7o7.7.7o7.7.7.7.7o7.7o7.7o7o7.7.7.7o7o7.%7.7o7.7.7.7.7o7.7.7.7.7.7.7o7.7.7. 


if  (plotfns  ==  0)  %  0  for  no  plots,  anything  else  to  create  plots 

status  =  ‘DONE’ 

else 

loglog(AR,Pr , ‘r’) 
figure ; 

plot(AR,T, ‘b’) 
figure ; 

plot (AR,U, ‘r’) 
figure ; 

loglog(AR,mF(l :b, 1) , ‘ or ’ ) 
figure ; 

loglog(AR,mF(l :b,2) , ‘b’) 


end 


B.0.3  isenfun.m 

function  F  =  isenfun(Tl , gas , array) 

vl  =  array(l) ; 
si  =  array(2) ; 

set (gas , ‘Density’ , 1/vl , ‘Temperature ’ ,T1) ; 
equilibrate (gas , ‘TV’ ) ; 
sgl  =  entropy .mass (gas) ; 
sg2  =  si; 

F  =  sg2-sgl ; 

B.0.4  nonideal_eq_soundspeed.m 

function  b  =  nonideal_soundspeed_eq(a,mech) 
7.  SOUNDSPEED  -  Speed  of  sound  (m/s)  . 

rhoO  =  density (a); 

VO  =  1/rhoO; 
tO  =  temperature (a) ; 
pO  =  pressure (a) ; 
sO  =  entropy_mass (a) ; 
qO  =  moleFractions (a) ; 

temp_gas  =  importPhase (mech) ; 
set (temp_gas , ‘T’ ,t0, ‘P’ ,p0, ‘X’ ,q0) ; 
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pi  =  1.001*p0; 

set (temp_gas , ‘P ’ ,pl , ‘S’ , sO , ‘X’ ,qO) ; 
equilibrate (temp_gas , ‘ SP ’ ) ; 
rhol  =  density (temp_gas) ; 

dpdrho_s  =  (pi  -  pO)/(rhol  -  rhoO) ; 
b  =  sqrt (dpdrho_s) ; 

B.0.5  nonideal  soundspeed.m 

function  b  =  nonideal_soundspeed(a,mech) 

7.  SOUNDSPEED  -  Speed  of  sound  (m/s)  . 

rhoO  =  density (a); 
tO  =  temperature (a) ; 
pO  =  pressure (a) ; 
sO  =  entropy_mass (a) ; 
qO  =  moleFractions (a) ; 

temp_gas  =  importPhase (mech) ; 

set (temp_gas , ‘T’ ,t0, ‘P’ ,p0, ‘X’ ,q0) ; 

rhol  =  1.001*rho0; 

set (temp_gas , ‘Density’ ,rhol , ‘Entropy’ , sO, ‘X’ ,q0)  ; 
pi  =  pressure (temp_gas) ; 

dpdrho_s  =  (pi  -  pO)/(rhol  -  rhoO); 
b  =  sqrt (dpdrho_s) ; 

B.0.6  oneDflow.m 

7.  oneDflow.m 

7.  Set  of  ODEs  to  solve  quasi  1-D  flow  with  area  change 

7. 

7,  FUNCTION  SYNTAX 

7  =============== 

7,  dydx  =  oneDf low(x,y ,gas ,P1  ,nf ile) 

7.  ALWAYS  called  by  an  ODE  solver,  in  this  case 

7.  out  =  odel5s  (@oneDf  low, xel  ,y0  .options  , gas  ,Ps  .nozzle) 

7.  where : 

7.  xel  =  distance  span 

7.  yO  =  initial  conditions  array  (pressure,  density,  velocity) 

7.  options  =  specified  in  program  that  calls  oneDflow.m 

7.  gas  =  Cantera  gas  object 

7.  Ps  =  Pressure  at  throat 

7.  nozzle  =  conical,  contour  or  constant 

7. 

7.  INPUT 
7.  ===== 

7.  pos  =  position 
7.  y  =  solution  array 
7.  gas  =  Cantera  gas  object 
7.  PI  =  pressure  at  throat 

7.  nfile  =  filename  of  contoured  nozzle  data 

7. 

7  OUTPUT 
7.  ====== 

7.  dydx  =  Array  of  ODEs  to  be  solved  by  ode  15s 
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7 

7  SUBFUNCTION  CALLS 
7  ================= 

7  Cantera  Functions:  set.m,  meanMolecularWeight  .m,  gas  const  ant  .m, 

"/.  density,  m,  nSpecies.m,  netProdRates  .m,  enthalpies_RT.m, 

7  molecularWeights .m,  cp_mass.m,  soundspeed .m, 

function  dydx  =  oneDf low(pos ,y ,gas ,P1 ,nf ile) 

x  =  pos;  7  initial  position 

Pinit  =  y(l);  7  initial  presure/pressure  at  sonic  condition  (PI) 

rinit  =  y(2);  7  initial  density 

Uinit  =  y(3);  7  initial  velocity 

set (gas , ‘Rho ’ , rinit , ‘ Y’ ,y (4 : end) ) ; 
avgW  =  meanMolecularWeight (gas) ; 

Rspec  =  gasconstant/avgW;  7  specific  gas  constant  [J/(kg*K)] 
nsp  =  nSpecies (gas) ; 

Tinit  =  Pinit*Pl/ (rinit*Rspec) ; 

set (gas , ‘ T’ , Tinit , 'Rho ’ , rinit , ‘ Y’ ,y (4 : end) ) ; 

7  Vectors 

wdot  =  netProdRates (gas) ; 

7  wdot  =  zeros (nsp);  7  test  -  compare  to  no  rxns 
Hi  =  enthalpies_RT(gas) *gasconstant*Tinit ; 
molwt  =  molecularWeights (gas) ; 

7  Scalars 
cp  =  cp_mass (gas) ; 
c  =  soundspeed(gas) ; 
r  =  density(gas) ; 

U  =  Uinit*rinit/r ;  7  flow  velocity  from  continuity  r*U  =  rinit*Uinit 

M  =  U/c ;  7  Mach  Number 

eta  =  1  -  M~2;  7  Sonic  Parameter 

sigma  =  0; 
for  z  =  l:nsp 

wi  =  molwt (z) ; 
hi  =  Hi(z)/wi; 
wd  =  wdot (z) ; 
dyidt  =  wi*wd/r; 

sigma  =  sigma  +  (avgW/wi  -  hi/ (cp*Tinit) ) *dyidt ;  7  Thermicity 

end 

[Area,  dAdx]  =  areafun(x,nf ile) ; 

Pprime  =  -(r*U/Pl)*(sigma  -  U*dAdx/Area)/eta;  7  Pressure  Derivative 

rprime  =  -(r/U)*(sigma  -  (M~2)*U*dAdx/Area)/eta;  7  Density  Derivative 

uprime  =  (sigma  -  U*dAdx/Area)/eta;  7  Velocity  Derivative 


77777777777777777777777777777777777777777777777777777777777777777777777777777 

7  SET  UP  COLUMN  VECTOR  FOR  DYDX 

77777777777777777777777777777777777777777777777777777777777777777777777777777 

dydx  =  [  Pprime 


rprime 
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uprime 

zeros  (lisp,  1)  ]  ; 

7,  species  equations 
rrho  =  1 .0/density(gas) ; 

for  i  =  l:nsp 

dydx(i+3)  =  rrho*molwt(i)*wdot(i)/U; 
end 


B.0.7  areafun.m 

This  function  takes  the  distance  downstream,  ‘x’,  and  the  name  of  the  file  containing  the 
nozzle  area  function  data,  ‘nozfile’,  as  inputs  and  returns  the  appropriate  Areci(x)  and 
d"dx£a (x) •  The  inPut  fifo  format  can  be  found  in  subsection  B.1.1,  in  section  B.l, 

7. 

7.  x  =  position  in  the  x  direction 

7.  nozfile  =  name  of  the  file  containing  the  nozzle  area  function  data 

7 

7  INPUT  FILE  : nozfile . dat 
7. 

7  To  use  this  function,  x  must  be  non-dimensionalized  by  1  cm. 

7  i.e.  x*  =  x/0.01 

7  The  area  coming  out  has  been  non-dimensionalized  by  the  throat  area. 

7,  i.e.  A*  =  A/A_throat  =>  A  =  A*  x  A_throat 

7,  The  change  in  area  coming  out  has  been  non-dimensionalized  by  both. 

"/,  i.e.  dA*/dx*  =  (0.01  x  dA/dx)/A_throat 
7  =>  dA/dx  =  100  x  A_throat  x  dA*/dx* 

function  [area,dAdx]  =  areafun(xpos , nozfile) 

filnam  =  [nozfile,  ‘  .dat ’]  ; 

fid  =  fopen(f ilnam, ‘r ’ ) ; 
f  scanf  (fid,  ‘  ’/,s  ’  ,  2)  ; 

ATP  =  fscanf  (fid,  <y„14e’ ,5)  ; 
f  scanf  (fid,  ‘  ’/,s  ’  ,  2)  ; 

A  =  fscanf  (fid,  7,14e’  ,10)  ; 
f  scanf  (fid,  ‘  ’/,s  ’  ,  2)  ; 

B  =  fscanf  (fid,  7,14e’  ,31)  ; 
f close (fid) ; 

xstar  =  100*xpos; 

delta  =  0; 
delta2  =  0; 

if  xstar  <  ATP(l) 

coef  =  f lipud(A(l : 3) ) ; 
astar  =  polyval (coef , xstar) ; 
for  j=2:3 

delta  =  delta  +  A( j ) * ( j-1) *xstar" ( j-2) ; 

end 

elseif  xstar  <  0 

coef  =  flipud(A(4: 10) ) ; 


7,  only  two  strings  allowed  in  line  1 

7,  there  must  be  5  ATP  coeffs 

"/,  only  two  strings  allowed  in  line  3 

"/,  there  must  be  10  A  coeffs 

"/,  only  two  strings  allowed  in  line  6 

7,  there  must  be  31  B  coeffs 
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astar  =  polyval (coef ,xstar) ; 
for  j=5 : 10 

delta  =  delta  +  A( j ) * ( j-4) *xstar" ( j-5) ; 

end 

elseif  xstar  <  ATP(2) 

coef  =  f lipud(B(l : 3) ) ; 
astar  =  polyval (coef , xstar) ; 
for  j=2:3 

delta  =  delta  +  B ( j ) * ( j-1) *xstar" ( j-2) ; 

end 

elseif  xstar  <  ATP(3) 

coef  =  flipud(B(4: 10) ) ; 
astar  =  polyval (coef , xstar) ; 
for  j=5 : 10 

delta  =  delta  +  B ( j ) * ( j-4) *xstar" ( j-5) ; 

end 

elseif  xstar  <  ATP(4) 

coef  =  flipud(B(ll : 17) ) ; 
astar  =  polyval (coef , xstar) ; 
for  j  =  12 : 17 

delta  =  delta  +  B ( j ) * ( j  —  1 1 ) *xstar~ ( j-12) ; 

end 

elseif  xstar  <  ATP(5) 

coef  =  f lipud(B(18 : 24) ) ; 
astar  =  polyval (coef , xstar) ; 
for  j  =  19 : 24 

delta  =  delta  +  B ( j ) * ( j-18) *xstar~ ( j-19) ; 

end 

elseif  (ATP (2)  ==  0) 

coef  =  f lipud(B(l : 3) ) ; 
astar  =  polyval (coef , xstar) ; 
for  j=2:3 

delta  =  delta  +  B ( j ) * ( j-1) *xstar" ( j-2) ; 

end 

else 

coef  =  flipud(B(25 : 31) ) ; 
astar  =  polyval (coef , xstar) ; 
for  j=26 : 31 

delta  =  delta  +  B ( j ) * ( j-25) *xstar~ ( j-26) ; 

end 

end 

throat  =  pi* (0 . 015) ~2; 

area  =  astar*throat ; 
dAdx  =  delta* 100*throat ; 


B.0.8  NASA2Chemkin.m 


"/.  FORMAT 
"/.  Rec  1: 

7.  : 

7.  Rec  2: 


7 

7. 


7. 


7 i 
7 


species  name  A24,  col  1-24 

comments  (data  source)  A56,  col  25-80 

number  of  T  intervals  12,  col  2 

optional  identification  code  A6,  col  4-9 

chem  formulas,  sym,  nos  5(A2,F6.2),  col  11-50 

zero  for  gas  and  nonzero  for  condensed  phases  II,  col  52 

molecular  weight  F13.5  col  53-65 

heat  of  formation  at  298.15  K  [J/mol]  F13.5  col  66-80 
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7.  Rec  3:  temp  range  2F10.3  col  2-21 

7.  :  no  of  coeff  for  Cp~o/R  II  col  23 

7  :  T  exponents  in  empirical  eqn  for  Cp~o/R  8F5.1  col  24-63 

7.  :  H~o  (298 . 15)  -  H~o(0)  [J/mol]  F15.3  col  66-80 

7.  Rec  4:  first  5  coeff s  for  Cp~o/R  5D16.8  col  1-80 
7.  Rec  5:  last  three  coeff s  for  Cp~o/R  3D16.8  col  1-48 
7.  :  integration  consts  B1  &  B2  2D16.8  col  49-80 

7.  repeat  3,  4,  and  5  for  each  interval 

7. 

7. 

7.  C  Hf  :  C JP  v33  1955  pl25.  NSRDS-NBS  3  sec3  1970. 

7.  3  g  7/97  C  1.00  0.00  0.00  0.00  0 . 00  0  12.0107000  716680.000 

7.  200.000  1000.0007  -2.0  -1.0  0.0  1.0  2.0  3.0  4.0  0.0  6535.895 

7.  6 . 495031470e+02-9 . 649010860e-01  2 . 504675479e+00-l .  281448025e-05  1 . 980133654e-08 
7.-1.606144025e-ll  5 . 314483411e-15  8 . 545763110e+04  4 . 747924288e+00 

7.  1000.000  6000.0007  -2.0  -1.0  0.0  1.0  2.0  3.0  4.0  0.0  6535.895 

7.-1.289136472e+05  1 . 719528572e+02  2 . 646044387e+00-3 . 353068950e-04  1 . 742092740e-07 
7-2 . 902817829e-l 1  1 . 642182385e-15  8 . 410597850e+04  4 . 130047418e+00 

7  6000.000  20000.0007  -2.0  -1.0  0.0  1.0  2.0  3.0  4.0  0.0  6535.895 

7  4 . 432528010e+08-2 . 886018412e+05  7 . 737108320e+01-9 . 715281890e-03  6 . 649595330e-07 
7-2 . 230078776e-l 1  2 . 899388702e-16  2 . 355273444e+06-6 . 405123160e+02 

7 

7  need  to  change  fname  &  coefficients  Axx  &  Bxx 
clear  all; 
hp  =  fN2’  ; 

NASAfile  =  [hp , ‘ _highT . inp ’ ] ; 
writetxt  =  1; 

7 

7  READ  IN  NASA  HIGH  TEMP  C0EFFS 

7 

f id2  =  f open(NASAf ile , ‘r ’ ) ; 

for  m  =  1:3 

tline  =  fgets(fid2); 

temp2(m)  =  sscanf (tline (2 : 11) , ‘7d’ ) ; 

k  =  sscanf (tline (23) , ‘7d’ ) ; 

for  n  =  l:k 

b  =  23  +  5*n; 
a  =  b  -  4; 

EX  =  sscanf (tline (a:b) , f7e ’)  ; 
ex(m,n)  =  EX; 

end 

hf  =  sscanf (tline (66 : 80) , ‘7e ’) ; 
for  z  =  1:2 

txtline  =  fgets(fid2); 
for  x  =  1:5 
d  =  16*x; 
a  =  d  -  15; 
b  =  a  +  11; 
c  =  d  -  2 ; 
n  =  5*(z-l)  +  x; 
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num  =  sscanf  (txtline  (a:b)  ,  ‘‘/,e  ’  )  ; 
alpha  =  class (num) ; 
if  (alpha(l : 4)==‘ char ’ ) 
coeff(m,n)  =  0; 

else 

EX2  =  sscanf  (txtline(c  :  d)  ,  ‘  "/,d’ )  ; 
coeff(m,n)  =  num*10~EX2; 

end 

end 

end 

end 

temp2(4)  =  sscanf  (tline(12: 22)  ,  ‘7,d’  )  ;  %  end  of  range  3  [K] 

f close (f id2) ; 

plotfns  =  1;  "/,  0  for  no  plots,  anything  else  for  plots 

T1  =  temp2(l);  "/,  start  of  range  1  [K] 

T2  =  temp2(2);  "/,  start  of  range  2  [K] 

T3  =  temp2(3)  ;  "/,  start  of  range  3  [K] 

T4  =  temp2(4)  ;  "/,  end  of  range  3  [K] 

Tstart  =  Tl; 

Tbreak  =  5000;  "/,  mid  temperature  between  two  temp  ranges 

Tend  =  10000; 

Tlplot  =  Tl;  "/,  start  temperature  for  plots 

T2plot  =  Tend;  "/,  end  temperature  for  plots 

Nstart  =  1; 

Nbreak  =  Tbreak  -  Tstart  +  1 ;  7,  corresponding  index  number 

Nend  =  Tend  -  Tstart  +  1; 

Nlplot  =  Tlplot  -  Tstart  +  1;  70  corresponding  index  number  for  plots 

N2plot  =  T2plot  -  Tstart  +  1 ; 


exl  =  ex(l) 
ex2  =  ex(2) 
ex3  =  ex(3) 
ex4  =  ex(4) 
ex5  =  ex(5) 
ex6  =  ex(6) 
ex7  =  ex(7) 


l  exponents  of  T 


All  =  coeff  (1,1)  ; 
A12  =  coeff (1,2)  ; 
A13  =  coeff (1,3) ; 
A14  =  coeff (1,4) ; 
A15  =  coeff (1,5) ; 
A16  =  coeff (1,6) ; 
A17  =  coeff  (1,7)  ; 
Bll  =  coeff (1,9) ; 
B12  =  coeff (1,10) ; 


A21  =  coeff (2 , 1) 
A22  =  coeff (2,2) 
A23  =  coeff (2,3) 
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A24  =  coeff (2,4) ; 

A25  =  coeff (2,5) ; 

A26  =  coeff (2,6) ; 

A27  =  coeff (2,7) ; 

B21  =  coeff (2,9) ; 

B22  =  coeff (2,10) ; 

A31  =  coeff (3,1)  ; 

A32  =  coeff (3,2) ; 

A33  =  coeff (3,3) ; 

A34  =  coeff (3,4) ; 

A35  =  coeff (3,5) ; 

A36  =  coeff (3,6) ; 

A37  =  coeff (3,7) ; 

B31  =  coeff (3,9)  ; 

B32  =  coeff (3,10) ; 

i  =  1; 

T  =  Tstart; 

while  T  <  T2 

CpR(i)  =  All* (T~exl)  +  A12*(T~ex2)  +  A13*(T~ex3)  +  A14*(T~ex4)  + 

A15* (T~ex5)  +  A16*(T~ex6)  +  A17*(T~ex7); 

HRT(i)  =  -All*(T“exl)  +  A12* (T~ex2) *log(T)  +  A13*(T~ex3)  +  A14* (T~ex4) /2  + 
A15*(T“ex5)/3  +  A16*(T~ex6)/4  +  A17* (T~ex7) /5  +  Bll/T; 

SR(i)  =  -All* (T~exl) /2  -  A12*(T~ex2)  +  A13* (T" ex3) *log(T)  +  A14*(T~ex4)  + 
A15*(T~ex5)/2  +  A16*(T~ex6)/3  +  A17* (T~ex7) /4  +  B12; 

X(i)  =  T; 
i  =  i  +  1; 

T  =  T  +  1; 

end 


while  T  <  T3 

CpR(i)  =  A21* (T~exl)  +  A22*(T~ex2)  +  A23*(T~ex3)  +  A24*(T~ex4)  + 

A25* (T“ex5)  +  A26*(T~ex6)  +  A27*(T~ex7); 

HRT(i)  =  -A21* (T“exl)  +  A22* (T~ex2) *log(T)  +  A23*(T~ex3)  +  A24* (T~ex4) /2  + 
A25*(T“ex5)/3  +  A26*(T~ex6)/4  +  A27* (T~ex7) /5  +  B21/T; 

SR(i)  =  -A21*(T~exl)/2  -  A22*(T~ex2)  +  A23* (T~ex3) *log(T)  +  A24*(T~ex4)  + 
A25*(T~ex5)/2  +  A26*(T~ex6)/3  +  A27* (T~ex7) /4  +  B22; 

X(i)  =  T; 
i  =  i  +  1 ; 

T  =  T  +  1; 

end 


while  T  <  T4 

CpR(i)  =  A31* (T~exl)  +  A32*(T~ex2)  +  A33*(T~ex3)  +  A34*(T~ex4)  + 

A35* (T“ex5)  +  A36*(T~ex6)  +  A37*(T~ex7); 

HRT(i)  =  -A31* (T“exl)  +  A32* (T~ex2) *log(T)  +  A33*(T~ex3)  +  A34* (T~ex4) /2  + 
A35*(T“ex5)/3  +  A36*(T~ex6)/4  +  A37* (T~ex7) /5  +  B31/T; 

SR(i)  =  -A31* (T~exl) /2  -  A32*(T~ex2)  +  A33* (T~ ex3) *log(T)  +  A34*(T~ex4)  + 
A35*(T~ex5)/2  +  A36*(T~ex6)/3  +  A37* (T“ex7) /4  +  B32; 

X(i)  =  T; 
i  =  i  +  1; 

T  =  T  +  1; 

end 


X3  =  X(Nstart : Nbreak) . “2; 
X4  =  X(Nstart : Nbreak) ; 


80 


X5  =  ones (size (X3) ) ; 

AA  =  [  X3’  X4’  X5  ’  ]  ; 

[il ,  j  1]  =  size (AA) ; 

Y2  =  X(Nbreak+l :Nend) . “3; 

Y3  =  X(Nbreak+l :Nend) . “2; 

Y4  =  X (Nbreak+1 : Nend) ; 

Y5  =  ones (size (Y3) ) ; 

DD  =  [  Y2’  Y3’  Y4’  Y5’  ] ; 

[i2,  j2]  =  size(DD)  ; 

BB  =  zeros(  [il ,  j2]  )  ; 

CC  =  zeros(  [i2, j 1] )  ; 

alphal  =  [AA,  BB;  CC,  DD]  ;  "/,  solving  alphal*P  =  xx 

xxl  =  CpR(Nstart :Nend) ; 

P  =  (alphal  ’  *alphal)\ (alphal  ’  *xxl  ’)  ;  ’/,  P(l:3)  for  CpRl 

7,  P  (4 : 7)  for  CpR2 


coef 2  =  -P(4) ; 
coef 3  =  P(l)-P(5) ; 
coef4  =  P(2)-P(6) ; 
coef 5  =  P(3)-P(7) ; 

alpha2  =  [Tbreak~4  -Tbreak~4  Tbreak~3; 

4*Tbreak~3  -4*Tbreak~3  3*Tbreak~2; 

12*Tbreak~2  -12*Tbreak~2  6*Tbreak] ; 
xx2  =  [-coef 2*Tbreak~3-coef 3*Tbreak~2-coef 4*Tbreak-coef 5 ; 
-3*coef 2*Tbreak~2-2*coef 3*Tbreak-coef4; 

-6*coef 2*Tbreak-2*coef 3] ; 
c_hat  =  (alpha2 ’ *alpha2) \ (alpha2! *xx2) ; 

Pl(l)  =  c_hat (1)  ;  "/.  PI  for  CpRl,  HRT1 ,  SRI 

PI (2)  =  c_hat (3) ; 

P2(l)  =  c_hat (2)  ;  "/.  P2  for  CpR2,  HRT2 ,  SR2 

P2(2)  =  P (4)  ; 

for  n  =  1:3 

PI (n+2)=P(n) ; 

P2 (n+2) =P (n+4) ; 

end 


Tatm  =  298; 

Natm  =  Tatm  -  Tstart  +  1; 

HHatm  =  0 ; 

HH1  =  0; 

HH2  =  0; 

for  j=l:5 

HHatm  =  HHatm  +  (PI ( j )/ (6-j ) ) *Tatm~ (5- j ) ; 
HH1  =  HH1  +  (PI ( j ) / (6-j ) ) *Tbreak~ (5- j ) ; 
HH2  =  HH2  +  (P2( j ) / (6-j ) ) *Tbreak~ (5-j ) ; 
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end 

PI (6)  =  (HRT (Natm) -HHatm) *Tatm; 

P2(6)  =  (HH1  +  PI (6)/Tbreak  -  HH2)*Tbreak; 


SSatm  =  0 ; 

551  =  0; 

552  =  0; 

for  j=l:4 

SSatm  =  SSatm  +  (PI ( j )/ (5-j ) ) *Tatm~ (5- j ) ; 

551  =  SSI  +  (PI ( j ) / (5-j ) ) *Tbreak~ (5-j) ; 

552  =  SS2  +  (P2( j ) / (5-j ) ) *Tbreak~ (5-j ) ; 

end 


SSatm  =  SSatm  +  Pl(5)*log(Tatm) ; 

551  =  SSI  +  PI (5) *log(Tbreak) ; 

552  =  SS2  +  P2 (5) *log(Tbreak) ; 

PI (7)  =  SR(Natm) -SSatm; 

P2(7)  =  SSI  +  PI (7)  -  SS2 ; 


Ql(j)  =  PI  (6-j  )  ; 
Q2(j)  =  P2 (6-j )  ; 

end 

Q1 (6 : 7)  =  PI (6 : 7) ; 
Q2(6 : 7)  =  P2 (6 : 7) ; 


for  n=l:Nbreak 
T  =  X(n)  ; 

chemCpR(n)  =  Ql(l)  +  Q1(2)*T  +  Q1(3)*T~2  +  Q1(4)*T~3  +  Q1(5)*T~4; 
chemHRT (n)  =  Ql(l)  +  Ql(2)*T/2  +  (Q1 (3) *T~2) /3  +  (Q1 (4) *T~3)/4  + 

(Q1  (5)  *rT'4)/5  +  Ql(6)/T; 

chemSR(n)  =  Ql(l)*log(T)  +  Q1(2)*T  +  (Q1 (3) *T~2) /2  +  (Q1 (4) *T~3) /3  + 
(Q1 (5) *T"4)/4  +  Q1  (7) ; 
end 

for  n=Nbreak+l :Nend 
T  =  X(n)  ; 

chemCpR(n)  =  Q2(l)  +  Q2(2)*T  +  Q2(3)*T~2  +  Q2(4)*T~3  +  Q2(5)*T~4; 
chemHRT (n)  =  Q2(l)  +  Q2(2)*T/2  +  (Q2 (3) *T~2) /3  +  (Q2(4)*T~3)/4  + 
(Q2(5)*T"4)/5  +  Q2(6)/T; 

chemSR(n)  =  Q2(l)*log(T)  +  Q2(2)*T  +  (Q2(3)*T“2)/2  +  (Q2(4)*T~3)/3  + 
(Q2(5)*T"4)/4  +  Q2(7) ; 
end 


start  =  sprintf  (  ‘  "/,d’ ,Tstart)  ; 
endtemp  =  sprintf  ( '"/.d’ ,  Tend)  ; 
brtemp  =  sprintf  (  ‘"/.d’  .Tbreak)  ; 

fn  =  [hp, ‘ start , ‘ -’ , brtemp , ‘ endtemp, ‘ .txt ’] ; 
d  =  date; 

fid  =  fopen(fn, fw’ ) ; 

fprintf(fid,  [hp,‘\n’]); 

fprintf  (fid,  ‘CURVE  FITTED  ON  7.s\n\n’ ,d)  ; 

fprintf(fid,  [start,  ‘  ’ , brtemp, ‘  ’ .endtemp, ‘\n’]) ; 

fprintf (fid,  ‘  Low  range  A1-A7 :  \n’); 
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fprintf  (fid,  7.1. 8E  7.1. 8E  7.1. 8E  7.1.8E\n  7.1. 8E  7.1. 8E  7.1 . 8E\n\n’  ,  Ql)  ; 
fprintf (fid,  ‘  High  range  A1-A7:  \n’); 

fprintf  (fid,  ‘7.1. 8E  7.1. 8E  7.1. 8E  7.1.8E\n  7.1. 8E  7.1. 8E  7.1 . 8E\n’  ,Q2)  ; 
f close (fid) ; 

midtemp  =  num2str (Tbreak) ; 

if  (plotfns  ==  0)  7.  0  for  no  plots,  anything  else  to  create  plots 

status  =  ‘DONE’ 

else 

plot (X(Nlplot : N2plot) ,CpR(Nlplot :N2plot)  ,  ‘r’) 
hold  on 

plot (X(Nlplot : N2plot) , chemCpR(Nlplot :N2plot) , ‘b’ ) 
hold  on 

figure ; 

plot (X(Nlplot : N2plot) ,HRT(Nlplot :N2plot) , ‘r’) 
hold  on 

plot (X(Nlplot :N2plot) ,chemHRT(Nlplot :N2plot) , ‘b’) 
hold  on 

figure ; 

plot (X(Nlplot : N2plot) , SR(Nlplot : N2plot) , ‘r’) 
hold  on 

plot (X(Nlplot : N2plot) , chemSR(Nlplot : N2plot) , ‘b’) 
hold  on 


end 


B.0.9  Viscosity  Files 

viscmixC02.m 


function  mu  =  viscmixC02(x,m,mui) 

7.  viscmix  Compute  the  viscosity  of  a  mixture  of  N  species  based  on 
7.  C.R.  Wilke,  J.Chem.  Phys.,vol.l8  pp. 517-522,  1950 

7. 

7. 

7.  Input : 

7.  x  =  column  vector  of  length  N  containing  the  mole  fraction  of  each 
7.  specie 

7.  m  =  column  vector  of  length  N  containing  the  molecular  weight  of 
7.  each  specie 

7.  mui  =  column  vector  of  length  N  containing  the  viscosity  of  each 
7.  specie 

7. 

7.  Output : 

7.  mu  =  viscosity  of  the  mixture 

7.  Author : 

7.  Eric  Marineau 

7.  Caltech,  T5  Hypervelocity  Shock  Tunnel  Lab,  Galcit 

7.  26  July  2007 


7.  computation  of  G  matrix 
N  =  length (x) ; 

G  =  zeros(N) ; 
for  i=l:N 
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for  k=l:N 
if  i~=k 

miner  =  (  l+(mui  (i)  /mui  (k)  )  “0 . 5*  (m(k)  /m(i)  )  “0 . 25  )~2; 
denom  =  (  2“ (3/2) * (l+m(i) /m(k) ) “0 . 5  ); 

G(i,k)  =  immer /denom; 

end 

end 

end 

"/,  multiply  by  0.1  to  get  units  of  kg/m-s 
mu  =  sum (mui . /(1+G*x. /x) ) ; 

viscspec.m 

function  mui  =  viscspecC02(chem,T) 

"/„  viscspec  Compute  the  viscosity  of  each  specie  for  the  8  species  air 
"/„  model  using  coefficient  form  NASA  Ref  1232  and  data  from  The  viscosity 
’/,  of  argon  is  taken  from  M.  N.  Macrossan  and  C.  R.  Lilley, 

’/,  Physics  of  Fluids  v  15  (11),  Nov  2003 
/ 
l 
1 

' I,  Input : 

’/,  chem  =  structure  defined  in  gas_mixture  .m  containing  the  curve  fits 
‘/,  T  =  temperature  of  the  mixture 

’/,  Output : 

‘/,  mu  =  viscosity  of  the  mixture 

‘/,  Author : 

‘/,  Eric  Marineau 

‘/,  Caltech,  T5  Hypervelocity  Shock  Tunnel  Lab,  Galcit 

°/,  26  July  2007 

"/,  modified  to  add  C02  by  I  Leyva,  Jan  ’08 

"/,  Original  8  elements:  e+,  N2,  02,  Ar,  N,  0,  N0,N0+ 

"/,  Additional  elements:  C02,  CO,  C2,  C3,  CN,  C,  C+ 

‘/,  Ref:  Olynick.  ,  D.  et  al .  ,  J  of  Spacecraft  and  Rockets, 

°/,  Vol .  36,  No.  3,  May- June  1999 

"/,  Viscosity  using  curve  fit 
for  i  =  1:8 

mui(i)  =  0 . l*exp(chem. spec .MUFC(i) ) *T~ (chem . spec .MUFA(i) *log(T)  + 
chem. spec .MUFB(i) ) ; 
end; 

for  i  =  9:15 

mui(i)  =  0 . l*exp( (chem. spec .MUFA(i) *log(T)  + 
chem. spec .MUFB(i) ) *log(T)+chem. spec .MUFC(i) ) ; 
end; 

"/,  Viscosity  of  argon  T  >  1500 

Ti  =  [1500  2000  2230  3000  4000  4360  5000  6000] ; 
muAri  =  [7.3  8.9  9.6  11.8  14.5  15.4  17.1  19.5]*lE-5; 
mui(4)  =  spline (Ti, muAr i, T) ; 


B.l  Input/Output  Data  Files 

B.1.1  Nozzle  Geometry  Files 

The  nozzle  geometry  files  conform  to  the  following  guidelines: 
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The  first,  third  and  sixth  lines  can  ONLY  contain  the  coefficient  names,  immediately 
followed  by  the  format  of  the  numbers. 

There  are  exactly  5  ATP  coefficients,  10  A  coefficients,  and  31  B  coefficients. 

This  is  similar  to  the  old  NENZF  format,  except  that  the  area  transfer  points  (ATPs) 
are  given  first  instead  of  last. 

The  following  files  are  written  specifically  for  the  nozzle  geometries  used  in  T5. 

30  mm  Conical  Nozzle 

"/.ATP  1-5  (7.14e) 

0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
"/.Al-10  C/.14e) 

1.0000000e+00  0 . 0000000e+00  1 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
"/.Bl-31  (”/.14e) 

1 . 0000000e+00  1 . 6371275E-01  6 . 7004660E-03  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00 

30  mm  Contour  Nozzle 

"/.ATP  1-5  c/.14e) 

-1.0000000e+00  1 . 2000000e+00  1 . 6000000e+01  2 . 0000000e+02  1 . 0000000e+03 
"/.Al-10  C/.14e) 

1.0000000e+00  0 . 0000000e+00  1 . 0000000e+00  1 . 0000000e+00  0 . 0000000e+00 
1.0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
"/.Bl-31  (”/.14e) 

1 . 0000000e+00-6 . 1236300e-03  0 . 1620970e+00  0 . 7871880e+00  0 . 3527470e+00 
7 . 1633500e-03  3 . 3462400e-03-2 . 3991400e-04  6 . 1195700e-06  0 . 0000000e+00 
-1.1377500e+00  0 . 3739650e+00  3 . 9124500e-02-5 . 8480000e-04  3 . 0294200e-06 
6 . 1342700e-10-4 . 2136200e-ll  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00 

Constant  Area  Nozzle 

"/.ATP  1-5  c/.14e) 

0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
"/.Al-10  C/.14e) 

1.0000000e+02  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
"/.Bl-31  C/.14e) 

1.0000000e+02  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00  0 . 0000000e+00 
0 . 0000000e+00 


B.1.2  Nozzle  Flow  I/O  Files 
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Sample  Input  File 


8  •/. 

N: 0.001935  1 

N2 : 0 . 680602  % 

NO: 0.127611 
C: 0.000000 
0:0.111041 
C02: 0.000000 
CO: 0.000000 
02:0.078812 
highT_mix_ JO . cti 


shot 2440  •/. 

contour30mm  '/, 
5216  •/. 

53 . 8e6  •/. 

1  •/. 

0  •/. 

1  7. 


n  =  number  of  species 
species  q{l} 

species  q{2}  (n  lines,  one  for  each  species) 


"/„  mechanism  ‘mech’  in  CTI  format 
filename  for  output  ‘fname’,  usually  shot  number 
filename  of  nozzle  area  data  ‘nozzle’ 

Stagnation  Temperature  TO  [K] 

Stagnation  Pressure  P0  [Pa] 

outfile:  1  to  print  file,  0  to  bypass  printing 
plotfns:  1  to  make  plots,  0  to  bypass 
xfinal  [m]  -  position  to  stop  integration 


Some  comments  on  the  previous  input  file,  as  taken  from  the  notes  in  nozflowfun.m: 
q  =  Initial  Composition  MUST  use  capital  letters  -  STRING 

mech  =  Mechanism  File  name  in  CTI  format  -  STRING  (Generally  all  mechanism 
files  are  stored  in  ‘Program  Files/Common  Files/Cantera/data’  on  a  Windows 
Machine) 

fname  =  shot  number  used  output  file  name 

nozzle  =  nozzle  geometry  -  either  conical30mm.dat,  contour30mm.dat  or  con- 
stantArea.dat 

outfile  =  Output  File  Name  -  Enter  a  non-zero  number  to  generate  a  .pit  file 
OR  enter  ‘O’  (the  number  zero)  to  bypass  generating  an  output  file 

plotfns  =  Plot  Functions  -  Enter  a  non-zero  number  to  generate  plots  of  Pressure, 
Temperature,  Velocity  and  Species  Concentrations  vs  Area  Ratio  OR  enter  ‘0’ 
(the  number  zero)  to  bypass  generating  plots 
Sample  Output  File 

###################################### 

# 

#  ISENTRDPE  CALCULATION 

#  CALCULATION  RUN  ON  29-Mar-2008 

# 

#  INITIAL  CONDITIONS 

#  ================== 

#  STAGNATION  TEMPERATURE  (K)  5216 
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#  STAGNATION  PRESSURE  (Pa)  53800000 

#  DENSITY  (KG/M“3)  3.3768e+001 

#  SPECIES  MOLE  FRACTIONS:  N: 0.001935  N2: 0.680602  NO: 0.127611  C:0 

#  0:0.111041  C02:0. 000000  C0:0. 000000  02 

# 

###################################### 

THE  OUTPUT  DATA  COLUMNS  ARE: 

Volume  Pres  Sound  Vel 
[m~3/kg]  [Pa]  [m/s]  [m/s] 


0 . 0296 
0.0306 
0.0316 
0.0326 


53800000 

51697603 

49706735 

47852200 


1434.04 

1390.84 

1385.44 

1380.24 


0.00 

356.00 

500.60 

608.02 


###################################### 

# 

#  NOZZLE  CALCULATION 

#  THROAT  CONDITIONS 

#  ================= 

#  TEMPERATURE  (K)  4758 

#  PRESSURE  (Pa)  30205328 

#  DENSITY  (KG/M~3)  2.1002e+001 

#  SOUND  SPEED  [m/s]  1.3589e+003 

#  SPECIES  MOLE  FRACTIONS:  N: 0.000845  N2: 0.695317  NO: 0.114866  C:0 

#  0:0.092391  C02: 0.000000  CO: 0.000000  02 

# 

###################################### 

THE  OUTPUT  DATA  COLUMNS  ARE: 


Position  Density  Pressure 

[m]  [kg/m~3]  [Pa] 


1 . 000000e-003 
1 . 000016e-003 
1 . 000033e-003 
1 . 000049e-003 


2 . 100242e+001 
2 . 100143e+001 
2 . 100057e+001 
2 . 099983e+001 


Temperature 

[K] 

3 . 020533e+007 
3 . 020351e+007 
3 . 020192e+007 
3 . 020054e+007 


Velocity 

[m/s] 

4 . 757985e+003 
4 . 757921e+003 
4 . 757866e+003 
4 . 757818e+003 


###################################### 

# 

#  NOZZLE  CALCULATION 

#  EXIT  CONDITIONS 

#  ================= 

#  PRESSURE  [Pa]  24012 

#  TEMPERATURE  [K]  1165 

#  DENSITY  [kg/m “ 3]  7.1515e-002 

#  VELOCITY  [m/s]  3.6464e+003 

#  SOUND  SPEED  [m/s]  6.6856e+002 

#  MACH  [-]  5 . 4541e+000 

#  SPECIFIC  GAS  CONSTANT  [J/(kg*K)]  2.8818e+002 

#  GAMMA  [-]  1 . 3310e+000 

#  SPECIES  MOLE  FRACTIONS:  N: 0.000000  N2: 0.753972  NO: 0.072055  C:0 

#  0:0.000000  C02: 0.000000  CO: 0.000000  02 

# 

###################################### 


000000 

0.078812 


000000 

0.096582 


. 359930e+003 
.  359993e+003 
. 360049e+003 
. 360097e+003 


000000 

0.173973 


Sample  *.cti  File 
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# 

#  Pieced  together  from  highT_co2 . cti  &  highT_n2 . cti 

#  rxn  data  from  Joe  Olejniczak 

#  ref:  Mitcheltree  94  C02-N2  mechanism:  8  species,  12  rxns 

# 

#  originally  created  Thurs,  Nov  2,  2006 

#  corrected  on  Fri,  Mar  21,  2008 

# 

units(length  =  "cm",  time  =  "s",  quantity  =  "mol",  act_energy  =  "cal/mol") 

ideal_gas(name  =  "gas", 

elements  =  "  N  C  0  " , 

species  =  . .  N  N2  NO  C  0  C02  CO  02  """, 

reactions  =  "all", 

initial_state  =  state (temperature  =  300.0, 
pressure  =  OneAtm)  ) 


# - 

#  Species  data 

#  - 

species (name  =  "N", 
atoms  =  "  N : 1  ", 
thermo  =  ( 

NASA(  [  200.00,  5000.00],  [  2 . 551359510E+000 ,  -8 . 707774780E-005 , 

2 . 702513760E-008,  -3 . 855552400E-012 ,  8 . 158480060E-016 , 

5 . 609296790E+004 ,  3 . 926090980E+000]  ), 

NASA(  [5000.00,  10000.00],  [  3 . 467587680E-001 ,  7 . 155048950E-004 , 

-4 . 985209040E-008 ,  1 . 132633460E-012 ,  0 . 000000000E+000 , 

6 . 001740740E+004,  1 . 957078870E+001]  ) 

) 

) 

species (name  =  "N2" , 
atoms  =  "  N:2  ", 
thermo  =  ( 

NASA(  [  200.00,  5000.00],  [  3 . 367776150E+000 ,  6 . 124317000E-004 , 

-7 . 8903827 10E-008 ,  5 . 618762820E-012 ,  -8 . 265727480E-016 , 

-1 . 030630570E+003,  3 . 677876810E+000]  ), 

NASA(  [5000.00,  10000.00],  [  8 . 748225300E+000 ,  -1 . 609736510E-003 , 
1 . 849739840E-007 ,  -5 . 446528000E-012 ,  0 . 000000000E+000 , 

-9 . 938338770E+003 ,  -3 . 400417780E+001]  ) 

) 

) 

species (name  =  "NO", 

atoms  =  "  N:1  0:1  ", 

thermo  =  ( 

NASA(  [  200.00,  5000.00],  [  3 . 52970451E+000 ,  5.63745766E-004, 

-7 . 28596656E-008 ,  -4 . 51046388E-012 ,  1 . 10056467E-015 , 

9 . 90058181E+003 ,  5 . 07142356E+000]  ), 

NASA(  [5000.00,  10000.00],  [  6 . 11917281E+000 ,  -7.30414827E-004, 

1.06554271E-007,  -3 . 83975042E-012 ,  0.000000000E+000, 

6 . 23772099E+003 ,  -1 . 26114326E+001]  ) 

) 


) 


species (name  =  "C", 
atoms  =  "  C : 1  ", 
thermo  =  ( 

NASA(  [  200.00,  5000.00],  [  2 . 492008160E+000 ,  -3 . 176429740E-006 , 

1 . 230402030E-008 ,  1 . 653762660E-012 ,  -3 . 848031920E-016 , 

8 . 545333520E+004 ,  4 . 816987870E+000]  ), 

NASA(  [5000.00,  10000.00],  [  1 . 897753210E+000 ,  3 . 364812680E-004 , 

-4 . 262838170E-008 ,  1 . 883958820E-012 ,  0 . 000000000E+000 , 

8 . 619126870E+004,  8 . 797021620E+000]  ) 

) 

) 

species (name  =  "0", 
atoms  =  "  0:1  ", 
thermo  =  ( 

NASA(  [  200.00,  5000.00],  [  2 . 606288620E+000 ,  -8 . 921956950E-005 , 

1 . 919984270E-008,  8 . 465095090E-013 ,  -2 . 006709720E-016 , 

2 . 919541490E+004 ,  4 . 547141860E+000]  ), 

NASA(  [5000.00,  10000.00],  [  2 . 080292310E+000 ,  1 . 544325120E-004 , 

-1 . 012466240E-008,  1 . 699429180E-013 ,  0 . 000000000E+000 , 

2 . 998189400E+004 ,  8 . 172285430E+000]  ) 

) 

) 

species(name  =  "C02", 

atoms  =  "  C : 1  0:2  ", 

thermo  =  ( 

NASA(  [  200.00,  5000.00],  [  4.781580110E+000,  1 . 630974560E-003 , 

-2 . 232868160E-007 ,  -1 . 652436650E-011 ,  3 . 971671770E-015 , 

-4 . 882410360E+004 ,  -2 . 006776990E+000]  ), 

NASA(  [  5000.00,  10000.00],  [  1 . 158910560E+001 ,  -1 . 975370230E-003 , 
3 . 030562420E-007 ,  -1 . 214103190E-011 ,  0 . 000000000E+000 , 

-5 . 791598320E+004,  -4 . 809741630E+001]  ) 

) 

) 

species(name  =  "CO", 

atoms  =  "  C:1  0:1  ", 

thermo  =  ( 

NASA(  [  200.00,  5000.00],  [  3. 4051401 10E+000,  6 . 133646300E-004 , 

-8 . 026241880E-008,  -6 . 763063860E-012 ,  1 . 525103350E-015 , 

-1.433539910E+004,  4 . 192231220E+000]  ), 

NASA(  [5000.00,  10000.00],  [  2 . 641034490E+000 ,  9 .484758230E-004 , 

-1 . 607418360E-007,  9 . 666733490E-012 ,  0 . 000000000E+000 , 

-1 . 296441810E+004,  9 . 584425660E+000]  ) 

) 

) 

species(name  =  "02", 
atoms  =  "  0:2  ", 
thermo  =  ( 

NASA(  [  200.00,  5000.00],  [  3 . 523048030E+000 ,  6 . 331510000E-004 , 

-6 . 365018570E-008,  -4 . 741986880E-012 ,  1 . 052000240E-015 , 

-1 . 077941160E+003,  4 . 414943860E+000]  ), 

NASA(  [5000.00,  10000.00],  [  3 . 243591070E+000 ,  6 . 070521740E-004 , 

-4 . 59772 1320E-008,  2 . 632995420E-013 ,  -5.421010860E-020, 
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-2. 153368840E+002,  6 . 660544630E+000]  ) 

) 

) 


# - 

#  Reaction  data 

#  - 

#  Reaction  1 

reaction (  "C02  +  C02  <=>  CO  +  0  +  C02",  [6.90000E+018,  -1.5,  63275]) 

#  Reaction  1 

reaction (  "C02  +  CO  <=>  CO  +  0  +  CO",  [6.90000E+018,  -1.5,  63275]) 

#  Reaction  1 

reaction (  "C02  +  N2  <=>  CO  +  0  +  N2",  [6 . 90000E+018 ,  -1.5,  63275]) 

#  Reaction  1 

reaction (  "C02  +  02  <=>  CO  +  0  +  02",  [6.90000E+018,  -1.5,  63275]) 

#  Reaction  1 

reaction (  "C02  +  NO  <=>  CO  +  0  +  NO",  [6 . 90000E+018 ,  -1.5,  63275]) 

#  Reaction  1 

reaction (  "C02  +  C  <=>  CO  +  0  +  C",  [1 .40000E+019,  -1.5,  63275]) 

#  Reaction  1 

reaction (  "C02  +  N  <=>  CO  +  0  +  N",  [1 .40000E+019,  -1.5,  63275]) 

#  Reaction  1 

reaction (  "C02  +  0  <=>  CO  +  0  +  0",  [1 .40000E+019,  -1.5,  63275]) 

#  Reaction  2 

reaction (  "CO  +  C02  <=>  C  +  0  +  C02",  [2.30000E+017,  -1.0,  129000]) 

#  Reaction  2 

reaction (  "CO  +  CO  <=>  C  +  0  +  CO",  [2 . 30000E+017 ,  -1.0,  129000]) 

#  Reaction  2 

reaction (  "CO  +  N2  <=>  C  +  0  +  N2",  [2 . 30000E+017 ,  -1.0,  129000]) 

#  Reaction  2 

reaction (  "CO  +  02  <=>  C  +  0  +  02",  [2 . 30000E+017 ,  -1.0,  129000]) 

#  Reaction  2 

reaction (  "CO  +  NO  <=>  C  +  0  +  NO",  [2 . 30000E+017 ,  -1.0,  129000]) 

#  Reaction  2 

reaction (  "CO  +  C  <=>  C  +  0  +  C" ,  [3 . 40000E+017 ,  -1.0,  129000]) 

#  Reaction  2 

reaction (  "CO  +  N  <=>  C  +  0  +  N" ,  [3 . 40000E+017 ,  -1.0,  129000]) 

#  Reaction  2 

reaction (  "CO  +  0  <=>  C  +  0  +  0",  [3 . 40000E+017 ,  -1.0,  129000]) 
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#  Reaction  3 

reaction (  "N2  +  C02  <=>  N  +  N  +  C02" ,  [7 . 00000E+018 ,  -1.6,  113200]) 

#  Reaction  3 

reaction (  "N2  +  CO  <=>  N  +  N  +  CO",  [7 . 00000E+018 ,  -1.6,  113200]) 

#  Reaction  3 

reaction (  "N2  +  N2  <=>  N  +  N  +  N2",  [7 . 00000E+018 ,  -1.6,  113200]) 

#  Reaction  3 

reaction (  "N2  +  02  <=>  N  +  N  +  02",  [7 . 00000E+018 ,  -1.6,  113200]) 

#  Reaction  3 

reaction (  "N2  +  NO  <=>  N  +  N  +  NO",  [7 . 00000E+018 ,  -1.6,  113200]) 

#  Reaction  3 

reaction (  "N2  +  C  <=>  N  +  N  +  C" ,  [3.00000E+019,  -1.6,  113200]) 

#  Reaction  3 

reaction (  "N2  +  N  <=>  N  +  N  +  N" ,  [3.00000E+019,  -1.6,  113200]) 

#  Reaction  3 

reaction (  "N2  +  0  <=>  N  +  N  +  0",  [3.00000E+019,  -1.6,  113200]) 

#  Reaction  4 

reaction (  "02  +  C02  <=>  0  +  0  +  C02" ,  [2 . 00000E+018 ,  -1.5,  59750]) 

#  Reaction  4 

reaction (  "02  +  CO  <=>  0  +  0  +  CO",  [2 . 00000E+018 ,  -1.5,  59750]) 

#  Reaction  4 

reaction (  "02  +  N2  <=>  0  +  0  +  N2",  [2 . 00000E+018 ,  -1.5,  59750]) 

#  Reaction  4 

reaction (  "02  +  02  <=>  0  +  0  +  02",  [2 . 00000E+018 ,  -1.5,  59750]) 

#  Reaction  4 

reaction (  "02  +  NO  <=>  0  +  0  +  NO",  [2 . 00000E+018 ,  -1.5,  59750]) 

#  Reaction  4 

reaction (  "02  +  C  <=>  0  +  0  +  C",  [1 . 00000E+019 ,  -1.5,  59750]) 

#  Reaction  4 

reaction (  "02  +  N  <=>  0  +  0  +  N" ,  [1 . 00000E+019 ,  -1.5,  59750]) 

#  Reaction  4 

reaction (  "02  +  0  <=>  0  +  0  +  0",  [1 . 00000E+019 ,  -1.5,  59750]) 

#  Reaction  5 

reaction (  "NO  +  C02  <=>  N  +  0  +  C02",  [1 . 10000E+014,  0.0,  75500]) 

#  Reaction  5 

reaction (  "NO  +  CO  <=>  N  +  0  +  CO",  [5 . 00000E+012 ,  0.0,  75500]) 

#  Reaction  5 

reaction (  "NO  +  N2  <=>  N  +  0  +  N2",  [5 . 00000E+012 ,  0.0,  75500]) 


#  Reaction  5 
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reaction (  "NO  +  02  <=>  N  +  0  +  02",  [5 . 00000E+012 ,  0.0,  75500]) 

#  Reaction  5 

reaction (  "NO  +  NO  <=>  N  +  0  +  NO",  [1 . 10000E+014,  0.0,  75500]) 

#  Reaction  5 

reaction (  "NO  +  C  <=>  N  +  0  +  C" ,  [1 . 10000E+014,  0.0,  75500]) 

#  Reaction  5 

reaction (  "NO  +  N  <=>  N  +  0  +  N" ,  [1 . 10000E+014,  0.0,  75500]) 

#  Reaction  5 

reaction (  "NO  +  0  <=>  N  +  0  +  0",  [1 . 10000E+014,  0.0,  75500]) 

#  Reaction  6 

reaction (  "NO  +  0  <=>  02  +  N" ,  [8 . 40000E+009 ,  0.0,  19450]) 

#  Reaction  7 

reaction (  "N2  +  0  <=>  NO  +  N" ,  [6 . 40000E+014 ,  -1.0,  38370]) 

#  Reaction  8 

reaction (  "CO  +  0  <=>  02  +  C" ,  [3 . 90000E+010 ,  -0.18,  69200]) 

#  Reaction  9 

reaction (  "C02  +  0  <=>  02  +  CO",  [2. 10000E+010,  0.0,  27800]) 

#  Reaction  10 

reaction (  "CO  +  N  <=>  NO  +  C" ,  [2 . 86000E+008 ,  0.5,  53630]) 

#  Reaction  11 

reaction (  "CO  +  CO  <=>  C02  +  C",  [2 . 33000E+006 ,  0.5,  65710]) 

#  Reaction  12 

reaction (  "NO  +  CO  <=>  C02  +  N" ,  [4 . 59000E+005 ,  0.5,  12070]) 

B.1.3  Viscosity  Input  File 

vise  Jnputs2440. m 

"/,  Computation  of  free  stream  viscosity  in  T5  shock  tunnel 

"/,  Free  Stream  temperature  (K) 

"/.  Tinf  =  1276.2 
Tinf  =  1867.5 

’/,  Constants  for  Viscosity  Curve  Fits  (T>=1000  K)  taken  from 

"/,  Nasa  Reference  1232  (The  constants  for  viscosity  of  argon  are  missing) 

‘/,  The  viscosity  of  argon  is  taken  from  M.  N.  Macrossan  and  C.  R.  Lilley, 

"/„  Physics  of  Fluids  v  15  (11)  ,  Nov  2003 
7,  Original  8  elements:  e+,  N2,  02,  Ar,  N,  0,  N0,N0+ 

"/,  Additional  elements:  C02,  CO,  C2,  C3,  CN,  C,  C+  added  by  I  Leyva  Jan  ’08 
"/,  Ref:  Olynick.  ,  D.  et  al.,  J  of  Spacecraft  and  Rockets, 
y,  Vol .  36,  No.  3,  May- June  1999 

chem. spec .MUFA(l) =  0; 
chem . spec .MUFA(2) =  0.0203; 
chem . spec .MUFA(3) =  0.0484; 
chem . spec .MUFA(4) =  0; 
chem . spec .MUFA(5) =  0.0120; 
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chem . spec ,MUFA(6)=  0.0205; 
chem . spec ,MUFA(7)=  0.0452; 
chem . spec ,MUFA(8)=  0; 
chem. spec. MUFA(9)=  -.019527387; 
chem. spec. MUFA( 10)=  -.019527394; 
chem . spec ,MUFA(11)=  -8.4311e-3; 
chem . spec ,MUFA( 12)=  -8.4312e-3; 
chem . spec ,MUFA( 13)=  -8.3811e-3; 
chem . spec ,MUFA( 14)=  -8.3285e-3; 
chem . spec ,MUFA( 15)=  -8.3285e-3; 


chem . spec .MUFB (1)=  2.5; 
chem . spec ,MUFB(2)=  0.4329; 
chem . spec .MUFB (3) =  -0.1455; 
chem . spec .MUFB (4) =  0; 
chem . spec .MUFB (5) =  0.5930; 
chem . spec .MUFB (6) =  0.4257; 
chem . spec .MUFB (7) =  -0.0609; 
chem . spec .MUFB (8) =  2.5; 
chem. spec. MUFB (9) =  1.047818; 
chem. spec. MUFB (10)=  1.013295; 
chem. spec. MUFB(11)=  0.7876060; 
chem. spec. MUFB (12)=  0.7876090; 
chem. spec. MUFB (13)=  0.7860330; 
chem. spec. MUFB (14)=  0.7703240; 
chem. spec. MUFB (15)=  0.7703240; 

"/.  additional  elements:  C02,  CO,  C2,  C3,  CN,  C,  C+ 


chem. spec. MUFC(1)=  -37.4475; 
chem. spec. MUFC(2)=  -11.8153; 
chem . spec ,MUFC(3)=  -8.9231; 
chem . spec ,MUFC(4)=  0; 
chem . spec ,MUFC(5)=  -12.3805; 
chem . spec ,MUFC(6)=  -11.5803; 
chem . spec ,MUFC(7)=  -9.4596; 
chem . spec ,MUFC(8)=  -32.0453; 
chem. spec. MUFC(9)=  -14.32212; 
chem. spec. MUFC( 10)=  -13.97873; 
chem. spec. MUFC(11)=  -13.02680; 
chem . spec ,MUFC( 12)=  -12.8240; 
chem . spec ,MUFC( 13)=  -12.9406; 
chem. spec. MUFC( 14)=  -12.7378; 
chem. spec. MUFC( 15)=  -12.7378; 

"/.Number  of  atoms  of  j~th  element  per  molecure  of  the  ith  species 
chem.ALPIJ  =[00010;  20000;  02000;  00100;  10000; 

0  1  0  0  0;  1  1  0  0  0;  1  1  0  -1  0;  0  2  0  0  1;  0  1  0  0  1;  ... 

0000  2;  0000  3;  1000  1;  0000  1;  000  -1  1]; 

‘/,  Element  molecular  weight 
chem .  elem .  CMW  ( 1 )  =  14.008;  '/.N 

chem.elem.CMW(2)  =  16;  "/.0 

chem.  elem.  CMW  (3)  =  39.944;  "/.Ar 

chem. elem. CMW (4)  =  0.00054847;  °/,e- 

chem .  elem .  CMW  (5)  =  12.00;  °/,C 
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%  Computation  of  the  molecular  weight  of  species 
chem . spec . CGI  =  chem. ALPIJ*chem.elem.CMW’ ; 

"/.  Species  concentration  (mole  fractions)  from  NENZF 


chem . spec . CAPX ( 1 ) 

o 

1 

CD 

t — i 

II 

"/.e- 

chem . spec . CAPX (2) 

=  0.753972; 

”/.N2 

chem . spec . CAPX (3) 

=  0.173973; 

7.02 

chem . spec . CAPX (4) 

=  l.e-10 

"/.AR 

chem . spec . CAPX (5) 

=  l.e-10 

"/.N 

chem . spe  c . CAPX ( 6 ) 

=  l.e-10 

7.0 

chem . spec . CAPX (7) 

=  0.072055; 

7.  NO 

chem . spec . CAPX (8) 

o 

1 

(1) 

T — 1 

II 

"/.N0+ 

chem . spec . CAPX(9) 

=  l.e-10 

"/.C02 

chem . spec . CAPX (10) 

=  l.e-10; 

”/.C0 

chem . spec . CAPX (11) 

=  l.e-8 

”/.C2 

chem . spec . CAPX (12) 

=  l.e-8 

”/.C3 

chem . spec . CAPX (13) 

=  l.e-8 

”/.CN 

chem . spec . CAPX ( 14) 

=  l.e-8 

”/.C 

chem . spec . CAPX (15) 

=  l.e-8 

”/.C+ 

"/.  mixture  molecular  weight 

chem.CMA  =  chem. spec . CAPX*chem . spec  .  CGI ; 

chem . spec . CMU  =  viscspecC02(chem,Tinf ) 

chem .MU  =  viscmixC02(chem. spec . CAPX’ , chem. spec . CGI , chem. spec . CMU’ ) 


94 


95 


Appendix  C 

Reaction  Mechanisms 


A  number  of  different  mechanisms  were  considered  for  the  analysis  in  this  work.  A  summary 
of  the  reactions  and  relevant  Arrhenius  parameters  for  each  mechanism  considered  follows 
in  this  appendix. 


C.l  Micheltree  94 

C.1.1  Mixed  Gases 


Table  C.l:  CO2-N2  Mechanism:  8  species,  12  reactions 


Third  Body 

Reactants 

Products 

Rxn  No. 

see  #la 

C02  +  M 

CO  +  0  +  M 

#1 

see  #lb 

C02  +  M 

CO +  0  +  M 

#1 

see  #2a 

CO  +  M 

C  +  O  +  M 

#2 

see  #2b 

CO  +  M 

C  +  O  +  M 

#2 

see  #3a 

n2  +  m 

N  +  N  +  M 

#3 

see  #3b 

n2  +  m 

N  +  N  +  M 

#3 

see  #4a 

02  +  m 

O  +  O  + M 

#4 

see  #4b 

o2  +  M 

O  +  O  +  M 

#4 

see  #5a 

NO  +  M 

N  +  O  +  M 

#5 

see  #5b 

NO  +  M 

N  +  O  +  M 

#5 

none 

NO +  0 

02  +  N 

#6 

none 

N2  +  0 

NO  +  N 

#7 

none 

co  +  o 

O2  +  C 

#8 

none 

002  +  0 

O2  +  CO 

#9 

none 

CO  +  N 

NO  +  C 

#10 

none 

CO  +  CO 

CO2  +  C 

#11 

none 

NO  +  CO 

CO2  +  N 

#12 
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Table  C.2:  Arrhenius  parameters  for  each  reaction. 


Third  Body 

Cf 

[m3/(mole  x  s)] 

Vf 

©/ 

[K\ 

Rxn  No. 

co2 , 

CO, 

n2,  no, 

O2 

6.900  xlO18 

-1.50 

6.3275 

xlO4 

#la 

c, 

O,  N 

1.400  xlO19 

-1.50 

6.3275 

xlO4 

#lb 

co2, 

CO, 

n2,  no, 

O2 

2.300  xlO17 

-1.0 

1.2900 

xlO5 

#2a 

c, 

O,  N 

3.400  xlO17 

-1.0 

1.2900 

xlO5 

#2b 

co2 , 

CO, 

n2,  no, 

O2 

7.000  xlO18 

-1.6 

1.1320 

xlO5 

#3a 

C, 

0 ,  N 

3.000  xlO19 

-1.6 

1.1320 

xlO5 

#3b 

co2 , 

CO, 

n2,  no, 

O2 

2.000  xlO18 

-1.5 

5.9750 

xlO4 

#4a 

C, 

O,  N 

1.000  xlO19 

-1.5 

5.9750 

xlO4 

#4b 

co2, 

CO, 

n2,  no, 

O2 

1.100  xlO14 

0.0 

7.5500 

xlO4 

#5a 

C, 

O,  N 

5.000  xlO12 

0.0 

7.5500 

xlO4 

#5b 

8.400  xlO9 

0.0 

1.9450 

xlO4 

#6 

6.400  xlO14 

-1.0 

3.8370 

xlO4 

#7 

3.900  xlO10 

-0.18 

6.9200 

xlO4 

#8 

2.100  xlO10 

0.0 

2.7800 

xlO4 

#9 

2.860  xlO8 

0.5 

5.3630 

xlO4 

#10 

2.330  xlO6 

0.5 

6.5710 

xlO4 

#11 

4.590  xlO5 

0.5 

1.2070 

xlO4 

#12 

C.1.2  Carbon  Dioxide 

The  simplified  Mitcheltree  mechanism  for  CO2  uses  5  species  and  6  reactions.  They  are 
the  same  reactions  as  for  the  Mixed  Gases  mechanism  (8  species,  12  reactions)  without  any 
species  containing  nitrogen. 


C.2  Various  Sources  -  Carbon  Dioxide 


Table  C.3:  C02  Mechanism:  5  species,  6  reactions 


Third  Body 

Reactants 

Products 

Rxn  No. 

all 

C02  +  M 

CO  +  O  +  M 

#1 

see  #2a 

CO  +  M 

C  +  O  +  M 

#2 

CO 

CO  +  M 

C  +  O  +  M 

#2 

see  #2b 

CO  +  M 

C  +  O  +  M 

#2 

see  #3a 

O2  +  M 

O  +  O  +  M 

#3 

see  #3b 

02  +  M 

O  +  O  +  M 

#3 

see  #3c 

02  +  M 

O  +  O  +  M 

#3 

none 

CO  +  CO 

CO2  +  C 

#4 

none 

co +  0 

02  +  C 

#5 

none 

O2  +  CO 

CO2  +  O 

#6 

97 


Table  C.4:  Arrhenius  parameters  for  each  reaction  from  fort. 8. new. 


Third  Body 

[m3/(mole  x  s)] 

Vf 

©/ 

[K\ 

Rxn  No. 

c, 

0,  co2 ,  CO,  02 

6.900  xlO18 

-1.50 

6.3275 

xlO4 

#1 

C02,  02 

2.300  xlO17 

-1.0 

1.2900 

xlO5 

#2a 

CO 

3.400  xlO17 

-1.0 

1.2900 

xlO5 

#2b 

c,  0 

3.400  xlO17 

-1.0 

1.2900 

xlO5 

#2c 

co2 

,  CO ,  n2,  no, 

02 

7.000  xlO18 

-1.6 

1.1320 

xlO5 

#3a 

C,  0,  N 

3.000  xlO19 

-1.6 

1.1320 

xlO5 

#3b 

CO-2 

,  CO,  n2,  no, 

O2 

2.000  xlO18 

-1.5 

5.9750 

xlO4 

#4a 

C,  0 ,  N 

1.000  xlO19 

-1.5 

5.9750 

xlO4 

#4b 

co2 

,  CO,  n2,  no, 

02 

1.100  xlO14 

0.0 

7.5500 

xlO4 

#5a 

C,  0,  N 

5.000  xlO12 

0.0 

7.5500 

xlO4 

#5b 

8.400  xlO9 

0.0 

1.9450 

xlO4 

#6 

6.400  xlO14 

-1.0 

3.8370 

xlO4 

#7 

3.900  xlO10 

-0.18 

6.9200 

xlO4 

#8 

2.100  xlO10 

0.0 

2.7800 

xlO4 

#9 

2.860  xlO8 

0.5 

5.3630 

xlO4 

#10 

2.330  xlO6 

0.5 

6.5710 

xlO4 

#11 

4.590  xlO5 

0.5 

1.2070 

xlO4 

#12 
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Appendix  D 

Thermofit  data 


This  appendix  contains  the  thermofit  data  for  each  species,  together  with  the  Arrhenius 
parameters  for  each  reaction,  used  by  NENZF  to  calculate  the  nozzle  flow. 
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Appendix  E 

Checklists 


This  appendix  contains  dummy  T5  checklists  for  both  the  high  enthalpy  (normal  operating 
mode)  and  low  enthalpy  shots.  All  facility  data  recorded  on  the  checklists  from  this  project 
are  also  included. 

E.l  High  Enthalpy  Shots 

T5  Checklist,  Version  21,  October  28,  2008  Signature:  _ 

Operators:  _  Date:  _  Time:  _  Shot#:  _ 

Researcher:  _  Name  of  Test:  _ 

Action  Check  Signal 

Piston  Mass:  _  kg 


Diaphragm  Thickness:  _ 

_ Expected  burst  pressure:  MPa  Value:  _  MPa 

Position  Tunnel  correctly  LVDT1:  setting:  _  value:  _  O 

_ LVDT2:  setting:  value:  O  A 

Turn  off  blower  O 

Open  SV-24  (DT  pressure  gage)  O 

SV-23  (DT  vacuum  gage)  O 

HS-20  (DT  vacuum  pump  isolation  valve)  O 

Close  HS-19  (DT  vacuum  pump  vent)  O 

HS-21  (DT  vent)  O 

Turn  on  cooling  water  for  DT  pump  O 

Turn  on  PO/3  (DT  vacuum  pump) _  O 
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Open  V-17  (ST  vacuum  valve)  O 

Close  V-18  (ST  vacuum  pump  vent  valve)  O 

Turn  on  PO/2  (ST  vacuum  pump)  O 

Open  SV-25A  (ST  pressure  gage)  O 

SV-25B  (ST  pressure  gage)  O 

SV-26  (ST  vacuum  gage)  O 

_ V-14  (ST  Pressure  Isolation  valve) _  O 

Close  2R  pressure  isolation  valve  (near  PI)  O 

Open  V-ll  (High  pressure  CT  isolation)  O 

V-8  (CT  vacuum)  O 

V-7  (CT/2R  vacuum  valve)  O 

V-12  (CT/2R/PS  vacuum  valve)  O 

Close  ZS-4  manual  and  remote  (emergency  2R  vent)  O 

ZS-13  (CT  vacuum  pump  vent  valve)  O 

Turn  on  PO/1  (2R,  Piston  Space  and  CT  vacuum  pump)  O 

Open  SV-29  (CT  vacuum  gage)  O 

SV-27  (2R  vacuum  gage)  O 

SV-31  (PS  vacuum  gage)  O  D 

Open  CT  gauge(s)  O 

Open  Flex  hose  isolation  valve _  O 

_ Wait  30  -  60  minutes  for  vacuum _  O 

Check  vacuum  levels  DT:  _  rnniHg  O 

ST:  _  rnniHg  O 

CT:  _  rnniHg  O 

2R:  _  rnniHg  O 

_ PS:  mmHg  O 

Disconnect  nozzle  support  trolley  O 

Isolate  ST  vacuum  valve  (V-17)  O 

Turn  off  ST  vacuum  pump  (PO/2)  O 

Open  ST  vacuum  pump  vent  valve  (V18)  O 

Isolate  ST  vacuum  gauge  (SV-26) _  O 


Fill  ST 


setting: 


kPa  abs 


(multiply  setting  by  20) 


thumbwheel: 


value: 


Press  hand  switch  HS-C3  to  confirm  ST  pressure 

Close  V-14  (ST  isolation) _ 

Isolate  DT  pressure  gauge,  HS-24 


Close  ZS-20  (DT  vacuum  isolation) 

Turn  off  HS-3  (dump  tank  vacuum  pump) 
Open  vent  HS-19  (DT  vacuum  pump  vent) 
Close  V-8  (CT  isolation) 

Close  (ST  gages):  SV-25A  & 


Close  SV-29  (CT  vacuum  gage) 
Fill  CT  gas:  _ 


Absolute  pressure  desired: 


SV-25B 


vacuum: 


setting: 

thumbwheel: 

value: 


kPa  gage 
kPa  abs 
kPa  gage 
kPa  gage 


O  E 
O  F 
O  G 


Close  V-ll  (CT  high  pressure  isolation  valve) 
Close  V-7  (CT/2R  vacuum) 

Close  SV-23  (DT  vacuum  gage) _ 

WARNING  LIGHTS  ON  (PI). 

Close  SV-27  (2R  vacuum  gage) 

2R  Temperature  before _ °C 

Fill  2R: 


setting: 

thumbwheel: 


value: 


Verify  door  is  closed 


setting: 

thumbwheel: 


kPa  abs 
kPa  gage 
kPa  gage 


PSI  gage 
PSI  gage 


O  H 

o 

O  I,J 


value: 


PSI  gage 


Check  for  changes  in  tunnel  position  CT: 


CT-2R: 
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Toggle  Dry  Run  Relay  box  to  “Experiment’' 

Arm  “T5  das” _ 

Isolate  flex  hose 

Confirm  2R  pressure,  Hand  switch  HS-C1 
Isolate  2R  gage 

Confirm  isolation,  Hand  switch  HS-28 _ 

Close  SV-31  (PS  vacuum  gage) _ 

Confirm  CT  pressure,  Hand  switch  HS-C2 
Isolate  CT  gage 

Confirm  isolation,  Hand  Switch  HS-30 
2R  Temperature  after  _  °C 

Close  V-12  (2R  vacuum) 

Open  V-13  to  vent  CT  vacuum  pump 
Turn  of!  CT  vacuum  pump  (PO/1) 

READY  TO  FIRE 

Press  and  hold  down  firing  buttons 


value: 


kPa  gage 


O  K,L 
O  M 
O  N 


Time: 


O  O 

o 

O  P 
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E.2  Low  Enthalpy  Shots 

T5  Checklist  -  Low  ho,  Vers. 21b,  October  28,  2008  Signature:  _ 

Operators:  _  Date:  _  Time:  _  Shot#: 


Researcher:  _  Name  of  Test: 


Action 

Check  Signal 

Piston  Mass: 

kg 

Diaphragm  Thickness:  _ 

_ Expected  burst  pressure:  MPa  Value:  _  MPa 

Position  Tunnel  correctly  LVDT1:  setting:  _  value:  _  O 

_ LVDT2:  setting:  value:  O  A 

Turn  off  blower  O 

Open  SV-24  (DT  pressure  gage)  O 

SV-23  (DT  vacuum  gage)  O 

HS-20  (DT  vacuum  pump  isolation  valve)  O 

Close  HS-19  (DT  vacuum  pump  vent)  O 

HS-21  (DT  vent)  O 

Turn  on  cooling  water  for  DT  pump  O 

Turn  on  PO/3  (DT  vacuum  pump) _  O 

Open  V-17  (ST  vacuum  valve)  O 

Close  V-18  (ST  vacuum  pump  vent  valve)  O 

Turn  on  PO/2  (ST  vacuum  pump)  O 

Open  SV-25A  (ST  pressure  gage)  O 

SV-25B  (ST  pressure  gage)  O 

SV-26  (ST  vacuum  gage)  O 

_ V-14  (ST  Pressure  Isolation  valve) _  O 

Close  2R  pressure  isolation  valve  (near  PI)  O 

Open  V-ll  (High  pressure  CT  isolation)  O 

V-8  (CT  vacuum)  O 

V-7  (CT/2R  vacuum  valve)  O 

V-12  (CT/2R/PS  vacuum  valve)  O 
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Close  ZS-4  manual  and  remote  (emergency  2R  vent) 
ZS-13  (CT  vacuum  pump  vent  valve) 

Turn  on  PO/1  (2R,  Piston  Space  and  CT  vacuum  pump) 
Open  SV-29  (CT  vacuum  gage) 

SV-27  (2R  vacuum  gage) 

SV-31  (PS  vacuum  gage) 

Open  CT  gauge(s) 

Open  Flex  hose  isolation  valve _ 

Wait  30  -  60  minutes  for  vacuum 


Check  vacuum  levels  DT:  _  mrnHg 

ST:  _  mrnHg 

CT:  _  mrnHg 

2R:  _  mrnHg 

_ PS:  mrnHg 


Disconnect  nozzle  support  trolley 
Isolate  ST  vacuum  valve  (V-17) 

Turn  off  ST  vacuum  pump  (PO/2) 

Open  ST  vacuum  pump  vent  valve  (V18) 

Isolate  ST  vacuum  gauge  (SV-26) 

Close  V-8  (CT  isolation) 

Close  SV-29  (CT  vacuum  gage) _ 

Fill  ST  &  CT  gas:  _  setting:  _  kPa  abs 

(multiply  setting  by  20)  thumbwheel:  _ 

value:  _ 

Press  hand  switch  HS-C3  to  confirm  ST  pressure 

Close  V-14  (ST  isolation) _ 

Isolate  DT  pressure  gauge,  HS-24 
Close  ZS-20  (DT  vacuum  isolation) 

Turn  off  HS-3  (dump  tank  vacuum  pump) 

Open  vent  HS-19  (DT  vacuum  pump  vent) 

Close  (ST  gages):  SV-25A  &  SV-25B 


O 

O 

O 

o 

o 

O  D 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O  E,F,G 
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Close  V-ll  (CT  high  pressure  isolation  valve) 
Close  V-7  (CT/2R  vacuum) 

Close  SV-23  (DT  vacuum  gage) _ 

WARNING  LIGHTS  ON  (PI). 

Close  SV-27  (2R  vacuum  gage) 

2R  Temperature  before _ °C 

Close  vent  HS-19  (DT  vacuum  pump  vent) 
Turn  on  HS-3  (dump  tank  vacuum  pump) 
Open  ZS-20  (DT  vacuum  isolation) 

Fill  2R: 


Close  ZS-20  (DT  vacuum  isolation) 

Turn  off  HS-3  (dump  tank  vacuum  pump) 
Open  vent  HS-19  (DT  vacuum  pump  vent) 

Check  for  changes  in  tunnel  position  CT: _ 

Toggle  Dry  Run  Relay  box  to  “Experiment’' 
Arm  “T5  das” 

Arm  oscilloscope _ 

Isolate  flex  hose 

Confirm  2R  pressure,  Hand  switch  HS-C1 
Isolate  2R  gage 

Confirm  isolation,  Hand  switch  HS-28 _ 

Close  SV-31  (PS  vacuum  gage) _ 

Confirm  CT  pressure,  Hand  switch  HS-C2 
Isolate  CT  gage 

Confirm  isolation,  Hand  Switch  HS-30 
2R  Temperature  after  _  °C 

Close  V-12  (2R  vacuum) 

Open  V-13  to  vent  CT  vacuum  pump 
Turn  off  CT  vacuum  pump  (PO/1) 


Verify  door  is  closed 


setting: 

thumbwheel: 


value: 


PSI  gage 
PSI  gage 
PSI  gage 


CT-2R: 


value: 


O  H 

o 

O  I,J 


_  o  K,L 

_  o  M 

kPa  gage  O  N 

o 

O  O 

_  o 

O  P 

o 
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READY  TO  FIRE 

Press  and  hold  down  firing  buttons 


Time: 


o  o 
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E.3  Checklist  data 

E.3.1  Series  1 


Table  E.l:  Series  1  Checklist  data 


Shot# 

Date 

[y/m/d] 

ref 

Shot# 

oper 

diaphm 

[10-3”] 

Pi 

[MPa] 

LVDTi 

lvdt2 

DT 

ST 

CT 

PS 

2331 

2006/02/28 

1999 

BA 

242/185 

95.4 

+52 

+505 

0.2 

0.25 

0.2 

0.2 

2332 

2006/03/02 

1999 

BA 

249/197 

107.7 

+48 

+503 

0.15 

0.22 

0.07 

0.06 

2333 

2006/03/06 

1999 

BA 

250/197 

104.5 

+39 

+504 

0.15 

0.15 

0.07 

0.055 

2334 

2006/03/07 

1999 

BA 

249/197 

102.6 

+45 

+498 

0.2 

0.2 

0.08 

0.06 

2335 

2006/03/09 

1993 

BA 

251/178 

77.5 

+190 

+705 

0.2 

0.18 

0.075 

0.065 

2336 

2006/05/17 

1993 

BA 

251/179 

78.0 

+193 

+719 

0.2 

0.3 

0.5 

0.55 

2337 

2006/05/18 

1999 

BA 

251/199 

110.2 

+53 

+487 

0.2 

0.7 

0.07 

0.06 

Table  E.2:  Series  1  Checklist  data  (cont’d) 

Shot# 

n2 

[%] 

Ppart'l 

[kPaa] 

co2 

[%] 

PST 

[kPaa] 

He 

[%] 

Ppart'l 

[kPaJ 

Ar 

[%] 

Ppart'l 

[kPaJ 

PCT 

[kPa9] 

T\ 

[K] 

P‘2R 

[psi] 

t2 

[K] 

2331 

100 

85 

85 

98.5 

15 

17.5 

+16.0 

23.9 

1100 

25.4 

2332 

95 

80.75 

5 

85 

85 

98.5 

15 

17.5 

+16.0 

24.8 

1100 

26.8 

2333 

100 

54.09 

85 

98.5 

15 

17.5 

+16.0 

24.3 

1100 

26.1 

2334 

50 

42.5 

50 

69.55 

85 

98.5 

15 

17.5 

+16.0 

24.1 

1100 

25.8 

2335 

100 

100 

82 

76.0 

18 

17.0 

-7.0 

23.8 

800 

24.8 

2336 

100 

100 

82 

76.0 

18 

17.0 

-7.0 

25.4 

800 

26.7 

2337 

95 

80.75 

5 

83.45 

85 

98.5 

15 

17.5 

+16.0 

26.6 

1100 

28.7 

E.3. 2 

Series  2 

Table  E.3:  Series  2  Checklist  data 

Shot# 

Date 

[y/m/d] 

ref 

Shot# 

oper 

diaphm 

[10-3”] 

Pi 

[MPa] 

LVDTi 

lvdt2 

DT 

ST 

CT 

PS 

2433 

2007/10/10 

2331 

BA 

247/197 

63.6 

+55 

+498 

0.2 

0.25 

0.2 

0.08 

2434 

2007/10/11 

2331 

BA 

247/197 

72.7 

+50 

+499 

0.15 

0.3 

0.3 

0.2 

2435 

2007/10/12 

2331 

EA 

246/204 

88.9 

+47 

+516 

0.18 

0.2 

0.18 

0.07 

2436 

2007/10/16 

2331 

BEA 

252/206 

105.8 

+53 

+496 

0.35 

0.2 

0.2 

0.09 

2437 

2007/10/17 

2331 

BIA 

252/207 

106.2 

+53 

+504 

0.25 

0.23 

0.18 

0.07 

2438 

2007/10/25 

2331 

BA 

248/204 

107.3 

+55 

+303 

0.2 

0.23 

0.17 

0.05 

2439 

2007/11/14 

2364 

BEA 

246/188 

98.6 

+56 

+402 

0.2 

0.3 

0.3 

0.2 
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Table  E.4:  Series  2  Checklist  data  (cont’d) 


Shot# 

air 

[%] 

Ppart'l 

[kPaa] 

co2 

[%] 

Pst 

[kPaa] 

He 

[%] 

Ppart'l 

[kPaa] 

Ar 

[%] 

Ppart'l 

[kPaJ 

PCT 

[kPa9] 

T-\ 

[K] 

P‘2R 

[psi] 

t2 

[K] 

2433 

100 

85 

85 

98.5 

15 

17.5 

+16 

22.2 

1100 

24.1 

2434 

100 

85 

85 

98.5 

15 

17.5 

+16 

22.6 

1100 

24.6 

2435 

100 

85 

85 

98.5 

15 

17.5 

+16 

22.1 

1100 

24.1 

2436 

100 

85 

85 

98.5 

15 

17.5 

+16 

21.0 

1100 

23.3 

2437 

100 

60 

85 

98.5 

15 

17.5 

+16 

21.9 

1100 

23.8 

2438 

100 

85 

85 

98.5 

15 

17.5 

+16 

24.1 

1100 

26.3 

2439 

100 

98 

82 

80.4 

18 

17.6 

-2.0 

22.4 

1150 

24.7 

E.3.3  Series  3 


Table  E.5:  Series  3  Checklist  data 


Shot# 

Date 

[y/m/d] 

ref 

Shot# 

oper 

diaphm 

[10-3”] 

Pi 

[MPa] 

LVDTi 

lvdt2 

DT 

ST 

CT 

PS 

2440 

2007/11/19 

2365 

BA 

246/185 

104.0 

-46 

+398 

0.5 

0.2 

0.2 

0.09 

2441 

2007/11/20 

2365 

BA 

246/189 

94.0 

-53 

+391 

0.3 

0.2 

0.2 

0.08 

2442 

2007/11/21 

2365 

BA 

245/188 

93.6 

-52 

+406 

0.3 

0.28 

0.18 

0.06 

2443 

2007/11/21 

2365 

BA 

246/189 

92.5 

-49 

+404 

0.3 

0.25 

0.18 

0.06 

2444 

2007/11/30 

2365 

EA 

245/188 

90.2 

-49 

+411 

0.3 

0.2 

0.09 

0.07 

2445 

2007/12/03 

2365 

BA 

246/185 

103.8 

-60 

+404 

0.2 

0.25 

0.6 

0.7 

Table  E.6:  Series  3  Checklist  data  (cont’d) 


Shot# 

air 

[%] 

Ppart'l 

[kPaa] 

co2 

[%] 

PST 

[kPaa] 

He 

[%] 

Ppart'l 

[kPaa] 

Ar 

[%] 

Ppart'l 

[kPaa] 

PCT 

[kPa9] 

Ti 

[K] 

P‘2R 

[psi] 

t2 

[K] 

2440 

100 

117 

74 

72.5 

26 

25.5 

-2 

21.0 

1150 

22.7 

2441 

80 

93.6 

20 

108.94 

74 

72.5 

26 

25.5 

-2 

21.6 

1150 

23.4 

2442 

60 

70.2 

40 

100.88 

74 

72.5 

26 

25.5 

-2 

22.2 

1150 

24.2 

2443 

30 

35.1 

70 

88.78 

74 

72.5 

26 

25.5 

-2 

21.8 

1150 

24.1 

2444 

100 

76.69 

74 

72.5 

26 

25.5 

-2 

21.9 

1150 

23.8 

2445 

100 

76.69 

74 

72.5 

26 

25.5 

-2 

21.9 

1150 

24.1 

Ill 


E.3.4  Series  4 


Table  E.7:  Series  4  Checklist  data 


Shot# 

Date 

[y/m/d] 

ref 

Shot# 

oper 

diaphm 

[10"3”] 

Pi 

[MPa] 

LVDTi 

lvdt2 

DT 

ST 

CT 

PS 

2446 

2008/01/11 

2365 

IA 

245/187 

94.2 

-47 

+391 

0.3 

0.25 

0.2 

0.07 

2447 

2008/01/16 

2365 

SA 

223/191 

90.0 

-47 

+399 

0.2 

0.15 

0.15 

0.06 

2448 

2008/01/17 

2365 

SA 

225/196 

92.1 

-50 

+395 

0.2 

0.2 

0.2 

0.06 

2449 

2008/01/18 

2365 

BISA 

223/201 

89.0 

-61 

+403 

0.2 

0.2 

0.25 

0.23 

Table  E.8:  Series  4  Checklist  data  (cont’d) 

Shot# 

air 

[%] 

Ppart'l 

[kPaa] 

co2 

[%] 

Pst 

[kPaa] 

He 

[%] 

Ppart'l 

[kPaJ 

Ar 

[%] 

Ppart'l 

[kPaJ 

PCT 

[kPa9] 

T\ 

[K] 

P2R 

[psi] 

t2 

[K] 

2446 

100 

117 

74 

72.5 

26 

25.5 

-2 

22.8 

1150 

24.4 

2447 

100 

117 

74 

72.5 

26 

25.5 

-2 

23.4 

1150 

25.3 

2448 

100 

117 

74 

72.5 

26 

25.5 

-2 

23.2 

1150 

24.9 

2449 

100 

117 

74 

72.5 

26 

25.5 

-2 

23.3 

1150 

24.8 

E.3.5 

Series  5 

Table  E.9:  Series  5  Checklist  data 

Shot# 

Date 

[y/m/d] 

ref 

Shot# 

oper 

diaphm 

[lO"3”] 

LVDTi 

[MPa] 

lvdt2 

DT 

ST 

CT 

PS 

2452 

2008/02/14 

1198 

BSA 

none 

+208 

+990 

0.15 

0.75 

0.35 

0.15 

2453 

2008/02/20 

1198 

SA 

none 

+199 

+983 

0.13 

0.3 

0.3 

0.08 

2454 

2008/02/22 

1198 

SA 

none 

+210 

+998 

0.3 

0.3 

0.18 

0.06 

2455 

2008/02/26 

1198 

SA 

none 

+206 

+990 

0.2 

0.3 

0.2 

0.1 

2456 

2008/02/27 

1198 

SA 

none 

+203 

+1014 

0.2 

0.4 

0.3 

0.3 

Table  E.10:  Series  5  Checklist  data  (cont’d) 


Shot# 

n2 

[%] 

Pst/ct 

[kPaa] 

Ti 

[K] 

P2R 

[psi] 

t2 

[K] 

Prto 

[psig] 

Prtf 

[psig] 

2452 

100 

122 

22.7 

170 

23.0 

360.9 

294.3 

2453 

100 

122 

22.7 

170 

22.6 

162.26 

85.2 

2454 

100 

122 

22.5 

170 

22.5 

50.24 

26.28 

2455 

100 

122 

24.0 

170 

24.2 

15.57 

6.01 

2456 

100 

122 

24.3 

170 

24.6 

112 


113 
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